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Abstract 


The determination of the radial and latitudinal temperature and wind profiles of 
the solar corona is of great importance in understanding the coronal heating mechanism 
and the dynamics of coronal expansion. Cram (1976) presented the theory for the 
formation of the K-coronal spectrum and identified two important observations. He 
observed the existence of temperature sensitive anti-nodes at certain wavelengths in the 
theoretical K-coronal spectra. The anti-nodes are separated by temperature-insensitive 
nodes. Remarkably, Cram showed that the wavelengths of the nodes and anti-nodes are 
almost independent of altitude above the solar limb. Because of these features, Cram 
suggested that the intensity ratios at two anti-nodes could be used as a diagnostic of the 
electron temperature in the K-corona. Based on this temperature diagnostic technique 
prescribed by Cram a slit-based spectroscopic study was performed by Ichimoto et.al 
(1996) on the solar corona in conjunction with the total solar eclipse of 3 November 1994 
in Putre, Chile to determine the temperature profile of the solar corona, hi this thesis 
Cram’s theory has been extended to incorporate the role of the solar wind in the 
formation of the K-corona, and we have identified both temperature and wind sensitive 
intensity ratios. The instalment, MACS, for Multi Aperture Coronal Spectrometer, a fiber 
optic based spectrograph, was designed for global and simultaneous measurement of the 
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thermal electron temperature and the solar wind velocity in the solar corona. The first 
ever experiment of this nature was conducted in conjunction with the total solar eclipse of 
11 August 1999 in Elazig, Turkey. In this instrument one end of each of twenty fiber 
optic tips were positioned in the focal plane of the telescope in such a way that we could 
observe conditions simultaneously at many different latitudes and two different radial 
distances in the solar corona. The other ends of the fibers were vertically aligned and 
placed at the primary focus of the collimating lens of the spectrograph to obtain 
simultaneous and global spectra on the solar corona. By isolating the K-coronal spectrum 
from the spectrum recorded by each fiber the temperature and the wind sensitive intensity 
ratios were calculated to obtain simultaneous and global measurements of the thermal 
electron temperature and the solar wind velocity. We were successful in obtaining 
reliable estimates of the coronal temperature at many positions in the corona. This is the 
first time that simultaneous measurements of coronal temperatures have been obtained at 
so many points. However due to instrumental scattering encountered during observations, 
reliable estimates of the wind velocity turned out to be impossible to obtain. Although 
remedial measures were taken prior to observation, this task proved to be difficult owing 
to the inability to replicate the conditions expected during an eclipse in the laboratory. 
The full extent of the instrumental scattering was apparent only when we analyzed the 
observational sequence. Nevertheless the experience obtained from this very first attempt 
to simultaneously and globally measure both the wind velocity and the temperature on 
the solar corona have provided valuable information to conduct any future observations 
successfully. 
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Chapter 1 


INTRODUCTION 


1.1 Introduction to the solar corona 

The solar corona is the region which begins at an altitude of -2000 km above the 
solar photosphere, and then extends out to many solar radii. This region becomes visible 
to the naked eye during a solar eclipse as seen in figure ( 1 . 1 ). 



Figure (1.1). The solar corona as seen during the 
solar eclipse of February 26, 1998. Picture: NASA 
eclipse library 
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Accounts of solar eclipses date back to at least the fourteenth century BC. 
However it is not clear that the solar corona was ever mentioned. The first clear mention 
of the corona was from Kepler in his Astronomiae pars Optica, although he seemed to 
have attributed the corona to be a property of the moon. With the advent of photographic 
methods and spectroscopy in the nineteenth century experimental studies were initiated 
on the nature of the sun during eclipses. By then eclipse path predictions were also 
possible. Edmond Halley was the first to predict the path of an eclipse in 1715 using the 
Newton’s law of gravity. 

At the 1860 solar eclipse, comparisons of photographs by De la Rue and Secchi 
from different locations confirmed that prominences were solar in nature. Prominences 
are cloud-like features with a reddish color, often seen during solar eclipses off the sun’s 
limb. This confirmation also added weight to the fact the corona too was solar in nature. 

Further discoveries followed with each subsequent solar eclipse. In the solar 
eclipse of 1871, Janssen discovered and identified absorption features in the solar 
spectrum. In 1814, Fraunhofer had been the first to label certain absorption features in the 
spectrum of the sun, and Janssen was able to identify some of Fraunhofer’s lines in the 
coronal spectrum. Most of the Fraunhofer features in the solar spectrum are atomic lines 
that are formed in the solar photosphere, although some are caused by molecules in the 
earth’s atmosphere. Janssen’s discovery of the solar absorption lines in the coronal 
spectrum confirmed that the corona was also solar in nature. 
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At the same eclipse of 1871, through spectroscopic observations Lockyer showed 
that the coronal gases extended ~500,000 km above the solar limb. This observation 
provided the first clue that the corona was also very hot. Later observations have showed 
that the corona extended to over thirty solar radii above the solar limb. Because of this 
great spatial extent, it can be inferred that the coronal gas is very hot. As to how hot the 
gas actually is, we look to the coronal spectrum for an answer: this is a major aim of the 
present work. However it was already known in 1957 (from work of E.N. Parker) that the 
coronal temperature is so high that it has important dynamical effects. To see this, we 
note that the extended nature of the solar corona is a result of a competition between 
gravity and gas pressure: the radial gradient of gas pressure tends to force material away 
from the sun, whereas gravity tends to hold the material back. In the inner corona, within 
one or two solar radii of the surface, these opposing tendencies balance each other almost 
exactly. However, beyond a few solar radii, Parker (1965) showed that gravity loses the 
competition, and the corona accelerates away from the sun to form the solar wind. 
Spacecraft measurements of the wind flux show that if the corona were not continually 
re-supplied with material, the wind would empty the corona in a time-scale of a few days. 

With the dawn of the spectroscopic age it was now possible to identify the 
elements present in hot gases through spectrum analyses. Hot gases emit photons at 
wavelengths, which are precisely characteristic of the elements present in the gases. 
However the coronal spectrum posed a great problem because most of the wavelengths 
could not be identified with known elements on Earth. Edlen (1937) identified some of 
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the emission lines as originating in forbidden transitions of the element iron in highly 
stripped stages of ionization, FeX and FeXIV. Forbidden transitions are those which are 
not allowed by electric dipole selection ndes: the transitions must be induced by 
magnetic dipole or electric quadmpole interactions. For iron to be in the ionization stages 
FeX or FeXIV requires that the electron temperature in the coronal gas be at least 1.0 
MK. For the forbidden transitions to occur, the ambient gas must have densities that are 
lower than a certain value. The upper limit of the density depends on the details of the 
transition. Golub and Pasachoff (1997) give a very detailed account of the history that led 
to the discovery of the solar corona and its properties. 

Figure (1.1) shows the solar corona during an eclipse. It is also evident from this 
picture that the coronal intensity peaks around the solar equator and gradually diminishes 
as it approaches the solar poles. Due to the reduced coronal intensity around the solar 
poles, these regions are called coronal holes. Although the coronal light intensity is small 
in these holes compared to the equatorial region, the holes are also the region for the fast 
solar wind. Figure (1.1) also shows solar plumes that radiate out from the north pole and 
south pole of the sun tracing its magnetic field. 

In the solar wind, material consisting of protons, electrons and a mixture of heavy 
elements flow out from the sun into interplanetary space. This ionized flow is so highly 
conducting that it drags along the magnetic field from the sun. The spatial extent of the 
solar wind is unknown: eventually, the wind runs into the interstellar medium and is 
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terminated at a shock. Spacecraft has not yet located the termination shock of the solar 
wind, but estimates put it at a distance of order 100 AU from the sun. 

In essence the solar corona is an extremely hot (> million Kelvin), highly ionized gas 
surrounding the sun. This is visible only during a total solar eclipse as a white light 
region extending to several solar radii, displaying streamers, plumes and loops. Its 
appearance changes during the solar cycle. At solar maximum it consists of many 
structures around the disk, but at the solar minimum it is dominated by large coronal 
holes at each pole and sheet-like structure near the equator. 

The brightness of the solar corona surrounding the solar disk is composed of three 
main components, namely; 

a. K- ( Kontinuierlich ) corona: The K-corona is due to the scattering of the 
photospheric light incident on the rapidly moving free electrons in the solar corona. 
The free electrons are a result of the electrons stripped off the coronal gas elements 
due to extremely hot temperature in the corona. Due to the high thermal velocity of 
these electrons they tend to broaden the photospheric spectrum, which consists of 
narrow lines, giving rise to a continuous spectrum. This scattering process is 
commonly known as Thomson scattering. 
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b. E- ( Emission ) corona: The E-corona is due to the emission from coronal ions, 


especially in highly ionized states. These emission lines are in many cases forbidden 
transitions resulting from atomic transitions from the highly ionized ions in the 
corona. Although these transitions are difficult to replicate in laboratories, however, 
they are not over ruled by selection rules for atomic transitions. 

c. F- (Fraunhofer) corona: The F-corona is due to photospheric light being scattered 
by dust particles in interplanetary space. Dust particles which are refractory enough to 
be able to survive at radial distances of a few solar radii also emit radiation in the near 
and mid-infrared wavelength regimes. The F-corona dominates the visible coronal 
brightness from about three solar radii distance from the center of the sun and has an 
increasing contribution to the total coronal brightness at longer wavelength. 

Figure ( 1 . 2 ) shows the radial brightness distribution of the F-corona, the K-corona 
and the E-corona. From figure ( 1 . 2 ) it is evident that the K-corona dominates in the 
coronal brightness to a height of ~ 2.0 solar radii above the limb and the F-corona begins 
to dominate beyond that height. 
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Figure (1.2). This plot shows the variation of the E, 
K and F components of the solar coronal brightness 
with height above the solar limb. 


As for the reason for the existence of the solar corona there does not exist a firm 
answer. Observations of other stars have revealed that many stars have coronae and many 
others do not. Empirically, the solar corona is densest and hottest in regions of closed 
magnetic fields (“loops”), and is coolest and fastest moving in regions of open fields 
(“coronal holes”). These results indicate that magnetic fields play an important role in the 
coronal heating process. 
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Figure ( 1 . 3 ) shows a cross section of the sun with some of its prominent features 
and the temperature scales between the core, photosphere and the corona. The core is the 
innermost part of the sun where energy is generated by nuclear reactions and its 
temperature is -15,000,000 K. The photosphere is the visible surface of the sun and also 
the layer which emits the light the human eye sees. Its temperature is - 6,000 K. The 
corona is the outer most layer in the solar atmosphere and consists of highly rarefied gas. 
This layer begins at altitudes of about - 2000 km above the photosphere and its 
temperature is -1,000,000 K. The corona is visible to the naked eye only during a solar 
eclipse. The region between the photosphere and the corona is the chromosphere. The 
region demarcating the chromosphere and the corona is the transition region that 
features a sharp temperature rise from -10 4 K to -10 6 K within a range of - 500 km. The 
coronal streamers are large-scale magnetic structures in the corona. 

From the point of view of physics, the most interesting aspect of the corona is the 
huge increase in temperature that occurs as we move from the photosphere (~ 6000 K) up 
into the corona (-1.0 MK). Now, coronal material must be replenished every few days 
(due to the solar wind): the only source of this material is the photospheric gas. 
Therefore, the creation of the corona is a continual process that involves heating the 6000 
K gas to several million K. It is impossible for a thermal process, i.e. one based on solely 
on exchange of internal energy, to do this. Therefore, a variety of non-thermal processes 
have been considered in the literature (Narain and Ulmschneider, 1996). The energy 
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Figure (1.3). A cross section of the sun to highlight the prominent 
parts of it and the temperature scales between the core, 
photosphere and the corona. 

Picture: SOHO pictures library 


The most likely source of non-thermal energy input to the solar corona is 
associated with disturbances generated in the solar convective layers, as shown in figure 
(1.3). These disturbances in the solar interior are manifested as supergranulations, 
granulations, magnetic flux tubes, and wave motions on the solar photosphere and 
acoustic energy above the photosphere. Now the question is to determine, which of these 
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energy forms can propagate upwards and then dissipate in the solar corona. However this 
process is further compounded by the need for the widely differing power requirements 
for the different coronal structures. The following are some of the popular method of 
non-thermal energy deposition in the solar corona. 

The most favored candidates are the Ohmic dissipation of the coronal electric 
currents and the viscous dissipation of magnetohydrodynamic (MHD) waves and 
turbulence (Narain and Ulmschneider 1990). The coronal electric currents are believed 
to be generated by the twisting of the magnetic flux tubes and, depending on the 
resistivity of the coronal plasma, dissipate through Joule heating in the coronal plasma. 
The MHD disturbances are believed to propagate along the magnetic fields protruding 
from the solar surface into the solar corona and dissipate their energy. These MHD waves 
can be visualized as a combination of the longitudinal sound waves and the transverse 
Alfven waves. Alfven waves are waves generated by the disturbance of the magnetic 
field lines and the sound waves could be produced by pressure perturbation. See Zirker 
(1993) for a survey of coronal heating theories, and conditions a successful theory must 
satisfy. 


The sun is also occasionally the site of transient releases of energy in the form of 
flares. These events emit high-energy particles and shock waves that may create aurorae 
on the earth. Although the reason for the solar wind is well understood to be driven by 
gas pressure due to the high temperature of the gas, however, the driving mechanism that 
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causes the solar wind flow velocity to change from subsonic velocities to supersonic 
velocities during its flow in the solar corona is yet to be established. The first analytical 
treatment of the solar wind model is due to Parker (1965). It was based on a thermally 
driven wind model. However the large velocities predicted at 1.0 AU do not seem to 
match the physical parameters at the coronal base. Here again the treatment of the solar 
wind as a purely thermally driven wind cannot produce the experimentally observed high 
velocities for the solar wind and some form of other non-thermal energy deposition 
mechanisms are needed to describe its transformation into supersonic velocities. Solar 
Wind Nine (1998) and the references therein give a summary of the present situation. 

Since the corona itself does not contain any heat sources that could possibly heat 
the corona to above million degrees or drive the solar wind flow to supersonic velocities 
the possibilities need to be linked to various phenomena on the solar surface itself. And, 
as a natural consequence to link various solar surface phenomena with the properties of 
the solar corona it is important to measure the coronal properties as accurately as 
possible. In this regard two of the most important physical parameters on the solar corona 
are the temperature and the solar wind velocity at different latitudes and radii in the solar 
corona. Therefore measurement of the coronal temperature and the solar wind velocity, 
simultaneously and globally on the solar corona, would help in gaining further insight 
into the state of the matter and its dynamics in the solar corona and its association with 
the solar surface phenomena. 
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The focus of this research is on the measurement of the coronal temperature and 
the solar wind velocity, simultaneously and globally, on the solar corona from the 
measurement of a single observable quantity. In section (1.2) some of the methods used 
by other researchers for the determination of the coronal temperature and the solar wind 
velocity are discussed. 

1.2 Other methods to measure the solar wind velocity and the coronal 
temperature 

The following are some of the techniques adopted in the measurement of the 
coronal temperature and the solar wind velocity. 

(a) Measurements of the radial intensity distribution of the white light corona (Saito 
1965, Koutchmy 1971, Guhathakurta et al. 1992) allow one to derive the scale-height 
of electron density. This can be converted to an equivalent electron temperature if one 
assumes hydrostatic equilibrium. 

(b) Emission line intensities of various lines have been used to determine the 
temperature of coronal electrons (Gabriel 1971, Nakada et al. 1975, Tsubaki 1975, 
Guhathakurta et al. 1992). These emission line theories rely on the ionization 
balances, various excitation mechanisms and atomic constants. Some of these 
quantities are subject to major changes owing to new discoveries in atomic physics. 
Ionization balance occurs when the rate at which an ion loses electrons per second (as 
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a result of collisional ionization by free electrons) is balanced by the rate at which the 
ion gains electrons per second (as a result of recombination). The excitation 
mechanism that is assumed to be predominant is collisional excitation. This model is 
a simple two-level atom consisting of the ground level m and an excited state n. The 
population of the upper state is set by a balance between collisional excitation from m 
to n due to electron-atom collisions, and by spontaneous radiative decay from n to m 
via an allowed atomic transition. For the coronal conditions the induced emission 
from n to m, the photoexcitation from m to n and the collisional de-excitation from 
impacts with electrons are considered negligible. The atomic constants that are 
inherent in this model are the Einstein coefficients for spontaneous emission, 
statistical weights of the various atomic levels, collisional excitation coefficient and 
the collision strength factor. The latter is a slowly varying function of the incident 
electron energy and involves difficult quantum mechanical calculations. In this regard 
we have preliminary evidence to suggest that Thomson scattering could be partly 
responsible for the emission lines in the EUV region using the Thomson scattering 
code developed as part of this dissertation work. 

(c) Hara et al. (1994) used the ratio of the soft X-ray intensities in different energy 
bands to determine the coronal ‘color temperature’. They used the X-ray analysis 
filters onboard Yohkoh to formulate a filter ratio method to determine the temperature 
of the quiet corona and have reported a value of 2.7 MK (Hara et al. 1992). Here 
photon energy ratios are determined for different filters. These are again functions of 
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the wavelength and temperature dependent emissivities and the response 
characteristics of the filters. Hara et al. (1994) and Withbroe (1988) contains reviews 
of coronal temperature measurements by various other groups. 

(d) Withbroe et al. (1985) measured the temperature of coronal ions using the 
Doppler widths of emission lines. The Doppler width measures the wavelength shift 
due to motion of the gas responsible for the observed emission line, which could be 
determined by comparison with a calibration lamp for the same emission line. Here 
the motion is attributed to the kinetic temperature, which includes both thermal and 
non-therntal motions. However to isolate the thermal effect on the motion the non- 
thermal effects on the motion, such as motions due to turbulence and waves, have to 
be eliminated. 

(e) Radio observational techniques have also been used to infer the temperature of 
coronal electrons (Zirin 1966). This method involves the radio observation of the 
quiet sun in determining the coronal temperature. From each level of the solar 
atmosphere only emission at frequencies higher than the plasma frequency may 
escape. Therefore, by tuning on to different frequencies the radio emission down to 
the region where the plasma frequency matches to the tuned frequency could be 
measured. This will include levels on the solar atmosphere with different 
temperatures too. If these temperatures are irregular the hotter regions are heavily 
weighted. For low optical depth the radio emission intensity is a product of the source 
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function and the optical depth where optical depth is a measure of how far the 
radiation will travel before being absorbed or scattered, which is also a function of the 
frequency. The radiation in the radio region is assumed to be due to free-free 
emission neglecting other non-thermal processes where the free-free emission is due 
to deceleration of an unbound electron during its passage close to an ionized atom. In 
addition the Planck function is assumed to represent the source function. Here again 
these are beset with uncertainties associated with the optical depth effects, refraction 
and the emission mechanisms. Furthermore, the presence of magnetic fields could 
produce other cut-off-frequencies associated with gyroresonance, which complicates 
the propagation. This reversal of the above procedure is also possible where radar 
signals are directed at the corona and are bounced off the layer where the plasma 
frequency is greater than the frequency of the incoming signal. 

(f) Doppler dimming technique has been used to measure solar wind velocity. In this 
technique (Strachan et al. 1993) the aim is to measure the intensity of a 
chromospheric line which has been resonantly scattered off coronal material. If the 
coronal material were to be at rest relative to the chromosphere, then the resonant 
scattering would be optimized, and the intensity of the scattered line would be 
maximized. However in the event that the coronal material has a non- zero velocity 
relative to the chromosphere, there is more or less significant reduction in the overlap 
of the line profiles in the chromosphere and corona. As a result of this mis-match 
between chromospheric and coronal line profiles, the efficiency of the resonant 
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scattering is reduced. The scattered line becomes dimmer, the larger the velocity 
difference between corona and chromosphere. Quantitative interpretation of an 
observed dimming in terms of wind velocity is subject to a number of uncertainties: 
among these are the emission mechanism of the line, and the associated atomic 
constants. 

(g) Solar wind velocities have also been measured using remote sensing techniques. 
Here a known radio signal is made to propagate through the solar corona and through 
the detection of fluctuations in various properties of the radio signal is utilized to 
determine the motions of density inhomogeneities in the solar corona. For these 
remote sensing techniques, a variety of methods are available, depending on whether 
the distant source is a natural source (broadband) or artificial (narrow-band). The 
fluctuations which the solar corona creates in intensity or in phase or in frequency or 
in line-width contain information on coronal density inhomogeneities on a broad 
variety of length scales, from a few kilometers up to tens of thousands of kilometers 
(Yakovlev and Mullan 1996). Observations can also be made from widely spaced 
ground stations or from single stations (Coles et al. 1986, Watanabe & Schwenn 
1989, Efimov 1994). Here the results are acknowledged to be substantially lower than 
the Doppler dimming techniques (Strachan et al. 1993). These too could be beset with 
problems owing to solar radio interference, carrier signal broadening and the 
sensitivity to variations in the density of the ionosphere, hi addition they also detect 
the motions of density inhomogeneities that include not only the bulk outflow of the 
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solar wind, but may also include wave motions if the latter have a compressive 
component. A review on this subject is found in Bird (1982). 

In all of the techniques summarized above, the aim is to measure the temperature of 
the coronal electrons or ions, or the solar wind speed. The techniques depend in different 
ways on a number of parameters, each of which is subject to uncertainty. However the 
theoretical idea utilized in this dissertation work is based on an well-understood simple 
physical principle, which is based on Thomson scattering of the photospheric radiation by 
the coronal free electrons, which is responsible for the formation of the K-corona. Here 
the quantity measured, which is the intensity, is dispersed over a wavelength region. The 
scattered intensity over this wavelength region is dependent on both the coronal 
temperature and the solar wind velocity with certain regions in the wavelength spread 
being sensitive to either the temperature or the wind velocity. This feature of selective 
sensitivity is exploited to formulate temperature and velocity diagnostics. In chapter-2 the 
physical overview for the formation of the K-corona is explained, together with a 
methodology to derive the thermal electron temperature and the solar wind velocity from 
the K-coronal spectrum. 


17 



Chapter 2 


FORMATION OF THE K-CORONA AND THE FOCUS OF THIS RESERACH 

2.1 The K-corona 

The K-coronal spectrum, as described in chapter- 1, is the coronal brightness 
attributed to the photospheric spectrum scattered off the free electrons in the solar corona. 
The existence of the free electrons is due to the high temperatures prevalent in the solar 
corona. This scattering phenomenon is commonly known as Thomson scattering. This 
brightness is only revealed during a total solar eclipse. 

2.2 The photospheric spectrum 

Figure (2.1) shows the solar spectrum in the visible region, as it would appear 
above the earth’s atmosphere. The shape of the photospheric spectrum is similar in 
nature, however, with a much greater flux. As evident from the solar spectrum the visible 
region is replete with absorption lines. Among the deep and conspicuous absorption 
lines, which were labeled by Fraunhofer, are the calcium H and K lines at 3968.5 
angstrom and 3933.7 angstrom, respectively. The dark Fraunhofer lines observed in the 
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solar spectrum are caused by selective process of line scattering or line absorption, 
accompanied by continuum absorption. 


extra terrestrial irradiance VS observed wavelength 
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Figure (2.1). This is a plot of the extraterrestrial solar 
irradiance spectrum. This was obtained by starting with 
the ground-based Fourier Transform Spectrometer at 
the Me Math/Pierce Solar Telescope at Kitt Peak, 
Arizona, and then correcting for wavelength-dependent 
absorption in the Earth’s atmosphere 
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2.3 The scattering source 

The scattering source for the formation of the K-corona is the free electrons in the 
solar corona. For a corona of one million Kelvin these free electrons have a mean thermal 
velocity of ~ 5500 km/sec and given by equation (2.1). 


v = 


2kT 


where 


V m 

k = Boltzmann constant = 1.38 x 10" J/K 

- 31 

m = Electron mass = 9.11x10 kg 
T = Coronal temperature ~ 10 


( 2 . 1 ) 


In an equilibrium state the electron velocity distribution can be described as Maxwelian 
in the solar corona where the electron number density in the vicinity of the solar limb is 
~10 8 cm 3 . This means that due to frequent collisions the velocity distribution is the same 
in all coordinates. However cataclysmic phenomena such as flares or coronal mass 
ejection can cause a deviation from a Maxwellian distribution locally. Furthermore, in the 
case of dilute plasma where collisions are infrequent the velocities may differ in different 
directions. This may be more pronounced in the interplanetary space where the particles 
stream from the sun in one direction. 
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2.4 The effects on the incident radiation by the scattering source 

In section (2.3) it was shown that a free coronal electrons move with an average 
thermal velocity of ~5500.0 km/sec in a corona of million degree Kelvin. Such motions 
are called Thermal Motions. This also means that a photon of wavelength (X ) incident 
in an electron will be Doppler broadened by equation (2.2). 


y 

AX = X x — where 
c 

v ~ 5500 km/sec for T = 1.0 MK 
c = 3 x 10' km/sec 


( 2 . 2 ) 


From equation (2.2) for an incident wavelength at 4000.0 angstrom the 
broadening is -70.0 angstroms. This phenomenon is called Thermal Doppler 
Broadening. This broadening increases with increasing thermal velocities for the 
electrons, which in turn depend on the coronal temperature. This phenomenon could be 
expected to cause all the narrow absorption lines, as evident in figure ( 2 . 1 ), to smooth 
into a continuum upon incidence with the coronal electrons. However weak depressions 
are still expected in the vicinity of the deep Fraunhofer lines thus smoothening the small- 
scale depressions while preserving the large-scale depressions. This near-total 
obliteration of the deep Fraunhofer lines by the thermal Doppler broadening served as 
one of the first clues for a hot corona. From equation (2.1) and equation (2.2) it is 
evident that the extent of smoothening (AX) is related to the coronal temperature (T) by 
equation (2.3). 

AX ~ Vt (2.3) 
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Once the radiation is incident on the electrons the scattering process is 
independent of wavelength and equal numbers of photons are scattered forwards and 
backwards. That is the scattering is coherent in the rest frame of the electrons. Since the 
scattering source itself is in thermal motion the radiation reaching the observer is 
Doppler shifted. In addition to the thermal velocity the coronal electrons may also 
acquire velocity due to turbulence. 

2.5 The K-coronal spectrum 

The photospheric radiation scattered off the free coronal electrons is known as the 
K-coronal radiation. Cram (1976) presented a physical model for the formation of the K- 
corona based on the Thomson scattering of the photospheric radiation by the free 
electrons in the solar corona. In figure (2.2) a schematic drawing of Cram’s model is 
shown. Here a known electron density distribution, determined by other independent 
means, is assumed for all positions along the line of sight. The normal to the line of sight 
from the center of the sun is at a distance ( p ). 

Each point on the line of sight has certain number of electrons and is subject to 
collisions with photons that make up the photospheric spectrum. The electrons at a 
particular point on the line of sight are exposed to the same photospheric spectrum from 
all the points on the solar surface that they see and are subjected to Thomson scattering. 
However the photospheric radiation incident on these electrons are also subjected to the 
limb-darkening effect. That is the brightness of the incident radiation is decreased as the 
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radiation from the edges travel through thicker and thicker levels on the solar atmosphere 
to reach the electrons. The limb darkening is wavelength dependent. The limb-darkening 
effect is most noticeable in the visible wavelength region. 

The thermal Doppler broadened scattered radiation observed by an observer on 
earth, due to Thomson scattering, is the accumulation of all the radiation that is scattered 
in the direction of the line of sight of the observer by all the electrons lying along the line 
of sight. 



Figure (2.2). The photospheric radiation is emitted from a point S 
on the surface of the sun. The electrons at point P in the corona 
scatter some of the incident radiation towards a terrestrial 
observer at E, which is known as K-coronal radiation. 
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Figure (2.3) shows a modeled spectrum for the scattered wavelength for an 
independently measured electron number density, given by equation (5.7), along the line 
of sight at 1.1 solar radii from the center of the sun. Here an isothermal corona was 
assumed and the plot shows models for different isothermal coronal temperatures. 


Modeled K— coronal intensity YS Wavelength. 



Figure (2.3). This plot shows the modeled K-coronal 
intensity spectra against wavelength. The composite plots 
are for different isothermal temperatures assumed for the 
solar corona for a given line of sight at 1.1 solar radii. 


As predicted in section (2.3) the modeled K-coronal intensity spectrum, as 
depicted in figure (2.3) for a line of sight at 1.1 solar radii, is very much smoother than 
the shape of the radiation incident on the coronal electrons. From the shape of the 
incident radiation spectmm, as is shown in figure ( 2 . 1 ), the vertical spread in the intensity 
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reflects the existence of millions of narrow absorption lines together with wide absorption 
lines. These millions of narrow lines in the incident radiation have been completely 
washed away, as revealed in figure (2.3), giving rise to a continuous spectrum. However 
the large-scale depressions are preserved. 


Figure (2.4) and figure (2.5) show the modeled K-coronal intensity spectra for 
lines of sight at 1.3 and 1.5 solar radii, respectively. 


Modeled K— coronal intensity VS Wavelength 



Figure (2.4). Modeled K-coronal intensity spectrum against 
wavelength for different isothermal coronal temperatures for 
a given line of sight at 1.3 solar radii. 
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Modeled K— coronal intensity VS Wavelength 



Figure (2.5). Modeled K-coronal intensity spectrum against 
wavelength for different isothermal coronal temperatures 
for a given line of sight at 1.5 solar radii. 


From these modeled K-coronal intensity spectra for lines of sight at different 
heights off the solar limb, as depicted in figure (2.3), figure (2.4) and figure (2.5), Cram 
(1976) made two very important observations. 


(a) The existence of temperature insensitive nodes in the K-coronal intensity spectra 
for various isothermal coronal temperatures. 


(b) The remarkable independence of the positions of these nodes on the wavelength 
scale for lines of sight at various heights above the solar limb. 
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From the above two properties for the K-coronal spectrum it is now possible to 
implement a temperature diagnostic scheme for the measurement of the coronal 
temperature by means of measuring the K-coronal intensity spectrum. For this purpose 
suitable intensity ratios have to be determined at wavelength positions where sufficient 
temperature dependent parameterization of the shape could be considered. In this regard, 
from figure (2.3)-(2.5), it is evident that the suitable wavelength positions occur at 3850.0 
angstrom and 4100.0 angstroms, respectively, which are also anti-nodes. Figure (2.6) is a 
plot of the K-coronal intensity ratios at 3850.0 angstrom to 4100.0 angstrom, from the K- 
coronal models for different isothermal coronal temperatures, against the temperatures. 


1.70 
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Figure (2.6). This is a plot of the K-coronal intensity 
ratios at 3850 .0 to 4100.0 angstroms against the 
assumed isothermal coronal temperatures. The ratios 
were calculated from the modeled K-coronal spectra for 
different temperatures shown in figure (2.3). These 
values pertain to a line of sight at 1.1 solar radii. 
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It is apparent from figure (2.6) that the coronal temperature could be determined by 
measuring the K-coronal intensity ratio at 4100.0 to 3850.0 angstroms from an observed 
K-coronal spectrum. Also from figure (2.6) it can be shown that if the above K-coronal 
intensity ratio could be determined within uncertainty of ± 1.0% then the coronal 
temperature could be inferred within an uncertainty of ± 0.2 MK. However it is 
unfortunate that the intensity measurement at 3850.0 angstrom lies in the blue end of the 
spectrum, which is subject to atmospheric absorption. This requires special chemical 
coated optics for signal enhancement. Below 3700.0 angstrom difficulties arise from the 
wavelength dependent atmospheric extinction while above 4700.0 angstrom the spectrum 
becomes relatively temperature insensitive. 

The feasibility of this method for the determination of the electron temperature in the 
solar corona, as first suggested by Cram (1976), was demonstrated by Ichimoto et al. 
(1996). Ichimoto et al. (1996) used a slit spectrograph for spectroscopic observation on 
the solar corona during the total solar eclipse of 3 November 1994 in Putre, Chile. By 
isolating the K-coronal spectra form the observed total coronal spectra temperatures of 
1.07 MK and 1.71 MK were determined for the base of a coronal hole and a streamer, 
respectively. They also reported cirrus covering of the sky during the totality, which had 
caused some difficulties in the analyses. 
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2.6 Further explanation for the formation of the temperature sensitive anti- 
nodes 

To adequately answer the reason for the formation of the temperature sensitive 
anti-nodes three simple intensity profiles are considered as input intensity spectra in order 
to understand the scattered intensity spectra produced by the theory. Figure (2.7) shows 
an absorption line centered at 4000.0 angstrom with a FWHM (Full width at half 
maximum) of 40.0 angstrom on an otherwise uniform continuum. In this example the 
intensity spectrum shown by figure (2.7) is assumed to be the sun spectrum as observed 
from earth. 


Input -spectrum to the scattering code 



3000 3500 400Q 4500 5000 5500 

Wavelength in angatroma 


Figure (2.7). An input intensity spectrum containing an 
absorption line centered at 4000.0 angstrom with a 
FWHM of 40.0 angstrom in an otherwise uniform 
continuum. 
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Figure (2.8) shows the theoretical scattered K-coronal intensity spectrum for an 
input spectrum given by figure (2.7). This intensity spectrum was obtained by allowing to 
be scattered off the coronal electrons. The theory assumes an isothermal corona of 
temperatures 0.5 MK. 1.0 MK, 1.5 MK and 2.0 MK. The essential feature observed in 
the scattered intensity spectrum is that the smearing of photons by coronal electrons fills 
in the deepest absorption thus making the absorption at the line center shallower. This 
smearing process increases with increasing temperatures, which in essence increases the 
Thermal Doppler Broadening as given by equation (2.2). The essential point is that the 
intensity of radiation at line center increases as the temperature of the scattering electrons 
increases. 


Modeled scattered intensity VS wavelength 
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Figure (2.8). The theoretical scattered intensity spectrum 
by the coronal electrons for an input spectrum given by 
figure (2.7). The theory assumes an isothermal corona. 
The plots show that the smearing of the absorption line in 
the scattered increases with increasing temperatures. 


30 




Ill figure (2.8), the gradient d(Intensity)/d(Temperature) is positive at line center. 


hi order to ascertain the behavior for an input intensity spectrum with two identical 
absorption features consider figure (2.9) where two strong absorption lines are centered at 
400.0 and 4500.0 angstroms in an otherwise uniform continuum. 


Input spectrum to the scattering code 
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Figure (2.9). An input intensity spectrum containing two 
absorption lines centered at 4000.0 and 4500.0 angstroms 
with a FWHM of 40.0 angstrom in an otherwise uniform 
continuum. 


Figure (2.10) shows the theoretical scattered K-coronal intensity spectrum for an 
input spectrum given by figure (2.9). Here we see that whereas the centers of absorption 
lines are filled in more and more as the temperature increases, the continuum between the 
line centers become progressively fainter as temperature increases. Therefore, at 
wavelength between the two lines (4100 - 4400 angstroms), the intensity of the 
continuum has a gradient d(Intensity)/d(Temperature) which is negative. This is in 
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striking contrast to the positive gradient d(Intensity)/d(Temperature) at line centers. It is 
this different behavior between line and non-line intensities that gives rise to anti-nodes, 
and provides the basis for Cram’s method. The reason for the difference in the depths of 
the absorption lines is due to the effect of the temperature dependent limb-darkening 
coefficient. 
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Figure (2.10). The theoretical scattered intensity spectrum by the 
coronal electrons for an input spectrum given by figure (2.9). The 
theory assumes an isothermal corona. The difference in the depths of 
the absorption lines in the scattered spectrum is due to the wavelength 
dependent limb-darkening coefficient. 
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Now consider the implications of introducing an emission line between the two 
absorption lines hi an otherwise uniform continuum as shown in figure (2.11). Here the 
absorption lines are centered at 4000.0 and 4500.0 angstroms while the emission line is 
centered at 4250.0 angstrom with both the absorption and the emission lines with a 
FWHM of 40.0 angstrom. Figure (2.12) shows the theoretical scattered K-coronal 
spectrum for an input spectrum given by figure ( 2 . 11 ). 
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Figure (2.11). An input intensity spectrum consisting of two 
absorption lines centered at 4000.0 and 4500.0 angstroms while 
an emission line centered at 4250.0 angstrom and with a 
FWHM of 40.0 angstrom for the absorption and emission lines 
in an otherwise uniform continuum. 
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Modeled scattered intensity VS wavelength 
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Figure (2.12). The theoretical scattered intensity spectrum by the 
coronal electrons for an input spectrum given by figure (2.11). The 
theory assumes an isothermal corona. The difference in the depths 
of the absorption lines in the scattered spectrum is due to the 
wavelength dependent limb-darkening coefficient. 


With the introduction of the emission line it is evident that in the scattered K- 
coronal intensity spectrum, as shown in figure (2.12), the absorption lines become 
shallower while the emission line peak becomes reduced in strength. This again gives rise 
to the intensity on one side of the node to decrease while the intensity on the other side 
increase. In analogy, when the input spectra of the examples shown in figure (2.7), figure 
(2.9) and figure (2.11) are substituted by the photospheric spectrum the dips in the 
photospheric spectrum are progressively filled in while the higher points are 
progressively reduced, giving rise to anti-nodes. This is tine of an integral over 
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wavelength of the product of the photospheric spectrum and a Gaussian kernel whose 
half-width lies in the range of 50 - 200 angstroms. 

hi the solar spectrum, the anti-nodes in the wavelength region 3900 - 4500 
angstroms owe their existence mainly to two wide and strong absorption features: (1) a 
pair of lines labeled H and K by Fraunhofer at 3968.5 and 3933.7 angstrom, respectively; 
(2) the combined effects of the Fraunhofer G band at 4300 - 4320 angstroms and the 
Fraunhofer feature H y at 4340.5 angstrom. The G-band is a collection of closely spaced 
absoiption lines predominantly due to the diatomic radical CH and iron. 

In essence, the uniqueness of the spectral region around 4000 angstrom for the 
formation of the anti-nodes arises from the following coincidences. On the one hand, the 
H and K lines happen to provide essentially a single absorption feature with a width of 
tens of angstroms. On the other hand, the only molecular band that contributes 
significantly to the visible solar spectrum (the G-band: almost 100 detectable lines within 
a range of 60 angstroms) happen to lie within a few tens of angstroms of a strong atomic 
line H ,. Both of these strong absorption features are readily identifiable in the low- 

resolution solar spectrum shown in figure ( 2 . 13 ): the figure is taken from Cram (1976) 
with his theoretical K-coronal spectra superposed. Note that in low-resolution, the solar 
photospheric spectrum between 4000 and 4300 angstroms has corrugations that give the 
impression of “quasi-emission lines”. For this reason, when we were modeling the 
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scattering of synthetic solar spectra above, as shown in figure ( 2 . 11 ), we included an 
emission feature between the two absorption lines. 

Are there any other parts of the solar spectrum that might be useful to test Cram’s 
idea? It seems unlikely. If we go shortward in wavelength from 3800 angstrom, the solar 
flux becomes rapidly so low that it is difficult to obtain adequate signal to noise ratio. 
And if we go to longer wavelengths we do not find pairs of strong absorption features 
that are wide enough to remain detectable when we smooth them by many tens of 
angstroms. 

The difference in the intensity of the photospheric spectra given by figure ( 2 . 1 ) 
and figure ( 2 . 13 ) is due to the units of measurement. In figure ( 2 . 1 ) the photospheric 
intensity in given in units of ergs/sec. cnr. angstrom while in figure ( 2 . 13 ) the intensity is 
given in units of ergs/sec. cnr.angstrom.steradians. To reconcile the photospheric 
intensities between figure ( 2 . 13 ) and figure ( 2 . 1 ) the latter needs to be divided by the 
solid angle of the sun’s disk, which is 6.8 x 10 ' . The differences in the intensities of the 
K-coronal spectra, given by figure ( 2 . 3 ) and ( 2 . 13 ), owe to the different electron number 
density profiles and the line of sight chosen for the analyses. 
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Figure (2.13). This figure shows the theoretical K-coronal 
intensity spectra for four isothermal coronal 
temperatures whose electron density is given by equation 
(5.8). The absolute intensity of the disk center 
photospheric spectrum is also shown. The vertical bars at 
the bottom of the figure indicate the relative strengths of 
the three coronal lines “c”, and nine flash spectrum lines 
“f”. 

Reproduced from Cram (1976). 
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2.7 The influence of the solar wind on the K-coronal spectrum 

The solar wind is an inherent property of the hot solar corona, a constant outflow 
of gas into interplanetary space. There the coronal electrons in the vicinity of the solar 
limb can be assumed to have a bulk outflow velocity equivalent to the solar wind 
velocity. In the vicinity of the sun this flow is radial and begins to take a spiral pattern at 
large distances, which is attributed to the rotation of the sun. However these flow 
velocities are latitude dependent with the fastest outflows occurring in the coronal holes. 

Figure (2.14) is a schematic diagram to explain the redshift induced on the K- 
coronal spectrum by the radially outflowing coronal electrons due to the solar wind. For 
the illustration of the above consider three electrons located at positions 1 ,2 and 3 on the 
line of sight. All three electrons are embedded in a radially directed out flow: the 
magnitude of the outflow velocity is the same at three electrons, but the velocity vectors 
point in different directions. Electrons 1 and 3 are symmetrically located on either sides 
of electron 2 while electron 2 is located at the closest point on the line of sight to the 
center of the sun. Since they are all moving away from the sun the photospheric radiation 
incident on these electrons are redshifted in the frame of reference of the electrons. 
Furthermore, for Thomson scattering, the photospheric radiation scattered off these 
electrons is wavelength independent in the reference frame of the electrons and they 
scatter a fraction of this incident radiation along the line of sight of the observer. In 
figure (2.14), the electrons 1 and 3 have common line of sight velocity magnitudes, away 
and towards the observer, respectively. This allows for the scattered radiation off 
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electrons 1 and 3 to be red and bine shifted, respectively, in the frame of reference of the 
observer, thus rendering a canceling effect on the observed radiation from these two 
electrons. However the line of sight velocity of electron 2 is zero, which simply scatters 
the redshifted radiation incident on it without causing anymore wavelength shifts. Thus 
the Doppler shifting of the radiation scattered by these three electrons is dominated by the 
effect due to electron 2, which is an overall redshifting of the scattered radiation along the 
line of sight. 



Figure (2.14). This is a schematic drawing to illustrate the redshifting 
of the photospheric radiation scattered off the coronal electrons, 
which are radially flowing away from the sun at a uniform velocity 
(w). Electron 2 is located at the closest point on the line of sight to the 
center of the sun. Electrons 1 and 3 are located symmetrically on 
either sides of electron 2. 
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The magnitude of this redshift on the K-coronal spectrum for an incident radiation 
of 4000.0 angstrom on the coronal electrons with an outflow velocity of 300.0 km/sec in 
a million degree Kelvin Corona is, from equation (2.2), ~4.0 angstrom. Figure (2.15) is a 
reproduction of figure ( 2 . 1 ) with horizontal bars to emphasize the magnitude of the 
thermal Doppler broadening and the redshift associated with the solar wind velocity. For 
this illustration these quantities are calculated for an incident wavelength of 4000.0 
angstrom on the coronal electrons, which are radially outflowing at a velocity of 300.0 
km/sec in a million degree Kelvin corona. 


m 



250 


200 


150 


extra_terrestrial_Lrra.diance VS observed wavelength 


Intense Fraunhofer Line: 

l 38383 A ! 



100 



Redshift due to solar wind velocity 
as indicated by the red dot within the circle 




Thermal Doppler Broadening T ipi 


. r 

phi 


W!' 1 


iillfij 

II 

iiifi' 


_i I i i i L_i I i i i 


3000 


3500 


4000 4500 

angatroma 


5000 


5500 


Figure (2.15). This figure shows the scale of the thermal 
Doppler broadening and the redshift associated with the 
solar wind velocity in the wavelength scale in comparison to 
the features of the photospheric spectrum. The scale of the 
redshift is shown by a red dot, which is highlighted by a 
circle. This comparison is made for an incident wavelength 
of 4000.0 angstrom on the coronal electrons, which are 
radially outflowing at a velocity of 300.0 km/sec in a million 
degree Kelvin corona. 
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The effect of the solar wind velocity in the model for the formation of the K- 
corona was considered to be negligible by Cram (1976). However an interesting 
additional property was revealed by the inclusion of a constant radial solar wind velocity 
to Cram’s model. That is the temperature insensitive nodes, which were observed by 
Cram (1976), were now observed to shift in wavelength positions with increasing solar 
wind velocity while maintaining the remarkable independence of the nodal positions with 
height above the solar limb. Figure (2.16) and figure (2.17) are K-coronal intensity 
models for different isothermal coronal temperatures at a line of sight at 1 . 1 solar radii 
with assumed radial solar wind velocities of 400.0 km/sec and 800.0 km/sec, 
respectively. 


Modeled K— coronal intensity VS Wavelength 



Figure (2.16). Modeled K-coronal intensity 
against wavelength for different isothermal 
coronal temperatures, for line of sight at 1.1 
solar radii and a wind velocity of 400.0 km/sec. 
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Figure (2.17). Modeled K-coronal intensity 
against wavelength for different isothermal 
coronal temperatures, for line of sight at 1.1 
solar radii and a wind velocity of 800.0 km/sec. 


Comparing figure (2.16) and figure (2.17) with figure (2.3), which is a K-coronal 
model for zero wind velocity, it is evident that the temperature insensitive nodes shift in 
wavelength positions with increasing solar wind velocities. This newfound physical 
property pertaining to the shifting in the wavelength positions of the temperature 
insensitive nodes with increasing solar wind velocities can be utilized to construct wind 
sensitive diagnostics. 
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Figure (2.18) is a plot of the K-coronal intensity model for different wind 
velocities ranging from 0.0 km/sec to 900.0 km/sec in intervals of 100.0 km/sec for an 
isothermal coronal temperature of 1.0 MK with the line of sight at 1.1 solar radii. It is 
evident from figure (2.18) on the existence of a wind insensitive node at 4233.0 angstrom 
with the maximum spread in intensity (anti-node) at 3987.0 angstrom. This property 
holds for other isothermal coronal temperatures and heights above the solar limb too. 


Modeled K— coronal intensity VS Wavelength 



Figure (2.18). Modeled K-coronal intensity against 
wavelength for solar wind velocities of 0.0 kni/sec to 900.0 
km/sec in intervals of 100.0 km/sec, for line of sight at 1.1 
solar radii and an isothermal coronal temperature of 1.0 MK. 
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Tlius the summary of the theoretical observations on the modified version of the Cram’s 
model, with the inclusion of the solar wind velocity is, 

(a) The existence of temperature sensitive anti-nodes at certain wavelengths in the K- 
coronal intensity spectrum provides a meaningful temperature diagnostics. 

(b) The wavelengths of the anti-nodes are remarkably independent of heights above 
the solar limb. 

(c) The simultaneous existence of nodes in the spectrum, where temperature effects 
are negligible, enables us to develop a diagnostic of wind speed that is more or less 
independent of temperature. 

2.8 Simultaneous determination of the coronal temperature and the solar wind 
velocity 

From the theoretical properties observed from the modified Cram’s model it is 
evident on the existence of both temperature and wind diagnostics. Figure (2.19) is a plot 
of the intensity ratios at the anti-nodes 3850.0 to 4100.0 angstroms against isothermal 
coronal temperatures for a range of solar wind velocities from 0.0 km/sec to 900.0 km/sec 
and with the line of sight at 1.1 solar radii. It too reveals the feasibility of the coronal 
temperature measurement within an uncertainty of ± 0.2 MK if the above intensity ratio 
at the anti-nodes 3850.0 to 4100.0 angstroms could be experimentally determined within 
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an uncertainty of ± 1.0 %. The width at the low end of the temperature range is ~0.02 


MK over a wind velocity range of 0.0 km/sec to 900.0 km/sec. 


1(41 00 )/I( 3850) VS Temperature 



Figure (2.19). The intensity ratio at 3850.0 to 4100.0 
angstroms against isothermal coronal temperatures for a 
range of solar wind velocities with the line of sight at 1.1 
solar radii. The horizontal spread in the low end of the 
temperature scale is ~0.02 MK over the velocity range of 
0.0 km/sec to 900.0 km/sec. 


Figure (2.20) is a plot of the intensity ratios of the wind sensitive anti-node at 
3987.0 to the wind insensitive node at 4233.0 angstroms against solar wind velocities for 
a range of isothermal coronal temperatures with the line of sight at 1.1 solar radii. The 
solar wind velocity can be determined by first determining the temperature from figure 
(2.19) and then by plotting in figure (2.20) the intensity ratios of the wind sensitive anti- 
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node at 3987.0 to the wind insensitive node at 4233.0 angstroms against the wind 
velocities at this temperature. Here the solar wind velocity measurements are highly 
sensitive to the wind- sensitive intensity ratio measurements. An uncertainty of ± 1.0% 
in the wind- sensitive intensity ratio measurements can cause an uncertainty of ± 200.0 
km/sec at 0.50 MK. This uncertainty increases with increasing coronal temperatures. 
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Figure (2.20). The intensity ratio at 3987.0 to 4233.0 
angstroms against solar wind velocities for a range of 
isothermal coronal temperatures with the line of sight at 
1.1 solar radii. To determine the wind velocity, first the 
temperature has to be determined from figure (2.19). 
Then the K-coronal models need to be created for a 
range of wind velocities at this temperature and plotted 
in figure (2.20). 
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2.9 Global determination of the coronal temperature and the solar wind velocity 

In section (2.7) a methodology was developed to simultaneously determine both 
the coronal temperature and the solar wind velocity from the same wavelength dispersed 
K-coronal intensity spectrum. That is, by measuring the intensity ratios at 3850.0 
angstrom to 4100.0 angstrom and 3987.0 angstrom to 4233.0 angstrom, from the 
observed K-coronal spectrum. The former is termed the, ‘temperature-sensitive’ 
intensity ratio, which allows for the temperamre measurements from figure (2.19). The 
latter is termed the, ‘wind-sensitive’ intensity ratio, which allows for the wind 
measurements from figure (2.20). However for the wind measurement, further theoretical 
K-coronal models need to be created for various wind velocities using the measured 
temperature as the isothermal coronal temperature. From these models the wind-sensitive 
intensity ratios are plotted against the wind velocities assumed for those models in figure 
(2.20). Now by calculating the wind-sensitive intensity ratio from the observed K- 
coronal spectrum its corresponding wind velocity can be obtained. 

In order to extend this measurement globally on the solar corona the wavelength 
dispersed K-coronal intensities need to be measured simultaneously and globally on the 
solar corona. Chapter-3 explains an instrument designed for simultaneous and global 
measurement of the wavelength dispersed coronal intensity spectra from multiple 
locations on the solar corona. Chapter-4 discusses the results from the very first attempt 
to measure both the coronal temperature and the solar wind velocity, simultaneously and 
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globally, on the solar corona. This experiment was performed during the total solar 
eclipse of 11 August 1999 hi Elazig, Turkey. 

2.10 The essence of this dissertation 

In this dissertation the influence of the solar wind velocity was added to the 
Cram’s model to search for both coronal temperature and solar wind velocity diagnostics. 
Based on the feasibility for the measurement of the coronal temperature and solar wind 
velocity diagnostics, form a single spectroscopic observation of the solar corona, an 
instrument was designed that would globally and simultaneously measures both these 
quantities at twenty different locations on the solar corona (Reginald & Davila, Paper I, 
2000). The very first experimental effort was carried out in conjunction with the total 
solar eclipse of 1 1 August 1999 in Elazig, Turkey and the results were analyzed to obtain 
both the coronal temperatures and the solar wind velocities at multiple locations 
(Reginald & Davila, Paper II, 2000). 
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Chapter 3 


THE INSTRUMENT 
3.1 Overview of the instrument 

Basically, the instrument designed was a spectrograph to obtain coronal spectra of 
the solar corona during the total solar eclipse of August 11, 1999. 

However this design incorporated a new feature in the way the light was collected 
and fed to the spectrograph in obtaining spectra. The exception from the usual slit-based 
spectrograph was that light was collected at the focal plane of the telescope via fiber 
optics. The fiber optic tips at the focal plane of the telescope were positioned at various 
latitudes and radii on the image of the sun formed on the focal plane of the telescope 
during the solar eclipse. This method served two very important purposes. 

1. This way the instalment was able to obtain simultaneous coronal spectra from 
different latitudes and radii on the solar corona. 
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2 . Also as a complimentary the telescope was spared the task of carrying the weight of 
the spectrograph. This was because they were detached systems connected via fiber 
optic cables. 

In essence the instrument was composed of three vital components. 

1. The telescope imaged the sun during the eclipse on the focal plane of the telescope. 

2 . The fiber optic tips, positioned at various latitudes and radii on the sun’s image on the 
focal plane of the telescope, carried the coronal light to the spectrograph. 

3. The spectrograph produced simultaneous spectra from the light received by the 
individual fibers. 

The telescope used in this experiment was a Schnhdt-Cassegrain telescope with a 
primary mirror diameter of 12-inches and a focal ratio of 10.0. A focal reducer was 
attached to reduce the focal ratio to 6.3. Schnhdt-Cassegrain telescopes by design are free 
from coma, spherical aberration and astigmatism. The telescope also contained auto- 
tracking capability. Figure ( 3 . 1 ) is a picture of the telescope used in this experiment and 
produced by the Meade Corporation. Figure ( 3 . 2 ) is a picture of the focal reducer 
attached to the back end of the telescope that increased the field of view by 56.0 % at the 
expense of lowering the magnifying power by 36.0 %. 
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At the focal plane of the telescope was a glass plate embedded with twenty-one fiber 
optic tips with ten each spread hi equal angles in a circular loop corresponding to 1.1 and 
1.5 solar radii on the sun’s image and one at the center of the frame. The one in the center 
of the frame was used to record the background signal using the lunar shadow during the 
eclipse. Another four fibers were attached to a mercury calibration lamp. 



Figure (3.1). Meade 12-inch, F/10.0 Schniidt- 
Cassegrain telescope with auto tracking 
capabilities. 


51 







1 6.3 focal reducer 



Figure (3.3) shows a schematic of the glass plate, as it would appear when the sun is 
in focus. Figure (3.4) shows the picture of the fiber optic tips embedded on the glass 
plate. The front end of the glass plate, where the fiber ends were exposed to the coronal 
light, was roughened with the back end smoothened. This allowed for the focusing of the 
sun prior to the eclipse on the front end of the glass while observing the image going in 
and out of focus from the back end. The fiber in the center was to record the background 
signal during the eclipse. 


Figure (3.2). F/6.3 focal reducer reduces the focal length by a 
factor of 0.63. The device shown in figure (3.4) was placed at 
the focal plane behind the F/6.3 focal reducer. 
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Figure (3.3) A drawing depicting the location of the fiber tips on 
the image plane of the telescope. The inner and outer rings are at 
1.1 and 1.5 solar radii, respectively. During the eclipse the lunar 
shadow replaces the sun. 


In figure (3.4) is the glass plate imbedded with twenty-one fibers that was placed 
at the focal plane of the telescope. The coronal light entering the fibers, which were 
located at the focal plane of the telescope and shown in figure (3.4), were all vertically 
aligned at the other end. The other end substitutes the position of the slit in a slit-based 
spectrograph. This slit, made up of multiple fibers at regular spacing between neighbors, 
was placed at the prime focus of the collimating lens of the spectrograph, hi figures (3.5) 
and (3.6) are images of the spectrograph. 
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Figure (3.4). The glass plate imbedded with twenty-one fibers. The inner 
and the outer circles are located at 1.1 and 1.5 solar radii, respectively, on 
the image plane of the sun. The fiber at the center was to record the 
background signal during the eclipse. 


From figure (3.5) it is apparent that, apart from the capability of this setup to 
simultaneously produce spectra, this methodology also allowed for the telescope and the 
spectrograph to be two independent units connected via light fiber optic cables. This is 
not possible in a slit-based spectrograph. In such a case the spectrograph needs to be 
attached to the back end of the telescope with the slit at the focal points of both the 
telescope and the collimating lens of the spectrograph. The purpose of the collimating 
lens is to transmit collimated light to the diffraction grating for dispersion. 
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Fiber optic link to the spectrograph 


Fiber optic link to the telescope 


1 Fibers to the mercury calibration lamp 


Figure (3.5). This is a picture of the spectrograph. The coronal light was fed 
into the spectrograph by twenty-one fibers that were vertically arranged at a 
regular spacing. Its location in the spectrograph was at the focal point of the 
collimating lens. 


Figure (3.6) shows the inside of the spectrograph between the grating and the 
camera lens that focuses the dispersed light to the camera in the backend. The 
spectrograph body contained many features that made it light tight from stray light that 
could have found its way into the spectrograph. Also the setscrews ensured that the many 
parts of the spectrograph fitted in place during assembly as a single unit. 


55 










Camera aligned for first order 


Figure (3.6). This is a picture of the inside of the spectrograph between the 
grating and the camera lens. Baffling were used to prevent all orders but the 
first order from entering into the camera lens. 


In figure (3.7) is a picture of the CCD camera that was used in this experiment to 
record the light dispersed in the first order in the wavelength region of 3500.0 to 4500.0 
angstroms. A thermoelectric cooler cooled the CCD device. Baffling were placed to 
restrict only the dispersed light in the first order from entering the camera. All the other 
orders were subjected to multiple reflections within the saw-tooth shaped grooves in the 
baffles for absorption by the black-coated baffling material. 
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Figure (3.7). This is a picture of the thermoelectrically cooled CCD 
camera used in the experiment. The front end housed the camera lens 
that focused the dispersed light from the grating on to the camera. 

As a general overview of the instrument the following could be said. The 
telescope imaged the eclipsed sun on to a glass plate containing twenty-one fibers, which 
was located at the focal plane of the telescope. The fibers picked up the light from 
various locations on the sun’s corona and the moon’s center. This light was fed to a 
spectrograph via the fiber optic cables. Then the grating simultaneously dispersed lights 
from all the fibers. The camera located at the backend of the spectrograph recorded the 
dispersed light in the first order in the wavelength region of 3500.0 to 4500.0 angstroms. 
In view of this instruments ability to record simultaneous spectra it was deemed 
appropriate to name this system the Multi Aperture Coronal Spectrometer with the 
acronym (MACS). 
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3.2 Description of the optical elements in detail 

In the following subsections all the optical elements of MACS will be discussed 
in details. 

3.2.1 The telescope 

The telescope used in MACS was a 12-inch F/10.0 Schmidt-Cassegrain 
manufactured by the Meade Corporation. The f-number was reduced to F/6.3 using a 
focal reducer in order to increase the plate scale. The LX200 model used in MACS also 
featured auto tracking and a super wedge for polar alignment. By design the Schmidt- 
Cassegrain type telescopes eliminate image distortion due to coma, spherical aberration 
and astigmatism, which are explained below. These optical systems are also called 
anastigmats and consist of at least three optical elements. A Schmidt-Cassegrain 
telescope consists of a concave and convex mirrors as the primary and the secondary, 
respectively, and a thin aspheric corrector plate in the front. 

Spherical aberration: Parallel rays are not brought to a focus at a point, but 
along a line. Therefore off-axis rays are brought to a focus closer to the lens than are on- 
axis rays. 

Coma: Off-axis rays do not quite converge at the focal plane, thus rendering off- 
axis points with tails, reminiscent of comets, hence the name. 
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Astigmatism: Rays of different orientations having different focal lengths 
resulting in the distortion of the image. Rays of light from horizontal and vertical lines in 
a plane on the object are not focused to the same plane on the edges of the image. Off- 
axis points are blurred in their radial or tangential direction, and focusing can reduce one 
at the expense of the other. 

The plate scale (P) of the telescope is defined as the spatial area seen in the sky, 
which is measured in arcseconds, corresponding to the spatial scale in the focal plane of 
the telescope, which is measured in micrometers and shown in figure (3.8). 



Figure (3.8). Schematic drawing on the image formation by the 
primary mirror of the telescope. 
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From figure (3.8) the expression for angle A is given by equation (3.1). 


Tan(A) =— f31) 

F 

Since (A) is a small angle, expressing (A) in radians, where (Tan (A) ~ A), and then 
converting radians into arcseconds, where (1 radian = 206265) arcseconds, gives the 
following expression for (A), as shown in equation (3.2). 

A = 206265 xD ( 3 . 2 ) 

F 


From the definition of the plate scale (P) it could be written as shown in equation (3.3). 


D 

206265XD/F 

D 

206265 

F~ 

Dx — 

D 

206265 

arcseconds/pm 

DxF/# 


(3.3) 
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Substituting for the diameter of the primary D in pm (12-inch = 12x2.54 cm = 
12x2.54xl0 4 p m) and the f-number of the telescope (F/# = 6.3) in equation (3.3) gives 
the following value for the plate scale (P) as shown in equation (3.4). 


P = 0.107 arcseconds/pm 


(3.4) 


3.2.2 The fibers 

The fibers linking the focal plane of the telescope and the spectrograph 
were twenty-one fibers with four more connected to a mercury calibration lamp. Each of 
these fibers was 1.0 m in length with a core diameter of 200.0 p m. Figure (3.9) shows the 
transmission efficiency of the fibers, marked as UV/VIS, for the region of wavelength 
interest which was from 3500.0 - 4500.0 angstroms. These fibers were purchased from 
Oriel Corporation. 
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Figure (3.9) Transmission efficiency of the fibers used in 
MACS and marked as UV/YIS. The region of wavelength 
interest in MACS was from 3500.0 -4500.0 angstroms. 


All optical fibers transmit light signals by total internal reflection. That is, if a ray 
of light in a medium of refractive index ni strikes the interface with another medium of 
refractive index n 2 (n 2 < nl) at an angle 0 , and if 0 is greater than 0 C , the ray is totally 

reflected back in to the first medium where 0 C is given by equation (3.5). 
0 c =sin- 1 (n i /n 2 ) (3.5) 
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Acceptance cone 
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Figure (3.10). Schematic diagram depicting the light ray path at 
the focal plane of the telescope into a fiber optic cable. 


Figure (3.10) is a schematic drawing depicting the light ray path at the focal plane 
of the telescope into a fiber optic cable. At point (A) applying the Snell’s law gives the 
relationship between angles (a and a ) and refractive indices (n 0 and ni) by equation 

(3.6). 

sin(a)xn 0 =sin(a)xn 1 (3.6) 


Now from triangle ABC the expression for angle 0 is given by equation (3.7). 

0 = 90 - a 
.\sin(0) = cos(a) 



(3.7) 
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At the critical angle (0 = 0 C ) using equations (3.5), (3.6) and (3.7) the expression for the 
acceptance angle (a) in terms of the refractive indices is given by equation (3.8). 


sin (a) = 


V n i 2 “ n 2 


n f 


andn 0 ~ 1 


(3.8) 


The sine of the acceptance angle (a) is defined as the numerical aperture (NA). Also from 
figure (3.10), the expression for sin(a) is given by equation (3.9). 

777 D/2 

sin(a) = — j— ( 3 . 9 ) 

1 

V 

D 

1 

~ 2(F/#) 

The fibers used in MACS have a high purity UV grade silica core and silica glass 
cladding with a numerical aperture of 0.27. This corresponds to an acceptance cone angle 
(2a) of 3 1.0 degrees and an f-number (F/#) of 1.9. 

All the light rays entering the fiber within the acceptance cone is totally reflected 
at the core-cladding boundary inside the fiber and the signal is transmitted to the other 
end. That is, all light rays within the acceptance cone entering the fiber at an angle less 
than the acceptance angle (a) strike the core-cladding interface at an angle greater than 
the critical angle 0 C . This fulfils the requirement for total internal reflection at the core- 
cladding boundary. 
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The f-number of the telescope used in MACS was F/6.3, which corresponds to an 
acceptance cone of 9.1 degrees. However this value was much less than the acceptance 
cone of 31.0 degrees for the fibers. Therefore all the light that was focused onto each of 
the fibers at the focal plane of the telescope was expected to be transmitted. 

The emergent cone form the other end of the fiber is dependent upon the input 
illumination, the fiber properties and the layout of the fiber bundle. For long fibers, the 
fiber properties dominate, while for short fibers (< 1.0 to 2.0 m) the input conditions 
dominate. Since the fibers in MACS were 1.0 m in length the emergent cone was 
expected to closely resemble the input conditions. 

3.2.3 The lenses 

The spectrograph contained two lenses of diameter 8.0 cm each with f- 
numbers F/3.9 and F/2.0. The former functioned as a collimator to collimate the light 
given out by the fibers. This enabled the grating to receive collimated light for dispersion. 
The latter acted as a focusing lens to focus the dispersed light on to the camera. Spindler 
& Hoyer manufactured these achromatic lenses. An anti UV-reflective coating was 
applied on the lens surfaces in order to enhance the transmission below 4000.0 angstrom. 
Figure (3.11) shows the reflection properties for the anti UV-reflective coating. Very low 
reflection ensures greater transmission. This was achieved with ARB2 UV chemical 
coating. 
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Figure (3.11). Reflection properties for the anti UV-reflective 
coating on the collimating and the camera lenses used in the 
spectrograph. 
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3.2.4 The grating 


The dispersing element used in MACS was a square transmission grating 
of dimensions 6.0X6.0 cm 2 with 6000.0 lines per cm and anti UV-reflective coating for 
enhanced transmission. Figure (3.12) shows the transmission properties for the anti UV- 
reflective coating. This grating was manufactured by American Holographic. 


Efficiency 



Figure (3.12). Transmission properties at normal incidence for the 
anti UV-reflection coating on the transmission grating. Here the 
back surface reflection is not considered. The substrate material is 
fused silica. The black squares pertain to the grating used in 
MACS. 
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3.2.5 The camera 


The camera used in MACS was a CCD camera with a back thinned SITE 
512.0X512.0, 24.0pm square chip and controlled with PMIS camera-control software. 
PMIS contains a powerful command interpreter designed for automating repetitive or 
complex actions. The camera can also be thermoelectrically cooled from 50.0° - 60.0° C 
below ambient and was manufactured by Apogee Instruments, Inc. Figure (3.13) is a plot 
of the quantum efficiency (QE) curve for the for the CCD chip in the camera. 



Figure (3.13). Quantum efficiency (QE) curve for the back thinned SITE 
512.0X512.0, 24.0 pm square chip used in MACS. This chip also featured 
a full well depth in excess of 300,000.0 electrons, a dynamic range of ~ 90.0 
dB and a readout time of ~ 5.0 seconds with an ISA card. 


68 








3.3 Analytical expressions for the integration time, the spectral resolution and 


the spatial resolution 



Figure (3.14). Schematic diagrams of the optical components that make 
up MACS and the parameters that determine the spatial and the 
spectral resolutions. 


Figure (3.14) is a schematic diagram of the optical components that contributed to 
the operation of MACS and the physical parameters of the different components that 
determined the spectral and the spatial resolutions, hi MACS the fibers substituted for a 
slit in a slit based spectrograph. 
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Consider a single fiber of diameter (df) and approximate the circular end of the 
fiber to a square of side (df). The projected fiber width (w) and height (h) on the camera 
focal plane are, respectively, given by equation (3.10) and equation (3.11). 


w 


V 


v 1 y 


COS(P) 


cos(ot) 


(3.10) 


h = d. 


V 

v f l. 


In equation (3.10) ( (3 ) is the angle of 


incidence, which is the case for MACS. 


(3.11) 


incidence on the grating and is zero for normal 
The direction of the dispersion is along (w). 


For normal incidence the grating equation is given by equation (3.12). 

dsin(oc) = (3.12) 

In equation (3.12) (d) is the ruling density, (a ) the diffraction angle, (m) the diffraction 
order and (X ) the wavelength. 
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Figure (3.15). Schematic diagram showing the projected width and 
height of a single fiber on the camera plane. 


Figure (3.15) is a schematic diagram showing the projected width and height of a single 
fiber on the camera plane. 


For collimated light incident on the grating the angular dispersion is given by 
(da/dX), with the dispersion parallel to the rulings on the grating. Then the linear 
dispersion is given by ^ x (da/ d/y) ) where (f2) is the focal length of the camera lens. 
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Using equation (3.12) the angular dispersion (AD) and the linear dispersion (LD) 


are given, respectively, by equation (3.13) and equation (3.14). 


AD 


da 

d)i 


m 

dcos(a) 



(3.13) 


(3.14) 



Figure (3.16). Schematic diagram of the dispersing element showing 
the angular and the linear dispersions. 
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Figure (3.17). Schematic layout of the slit spectrometer. In this diagram 
collimated light is incident on the grating. 


From figure (3.17) the height and width of the entrance slit (single fiber) subtend angles 
co and cp on the sky and the collimating lens, respectively, and are given by equation 
(3.15) and equation (3.16). 


co = 



(3.15) 



(3.16) 
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If for the CCD detector the pixel width is (p) then it needs to be decided on how 
many pixels are to be allowed to match the image width (w). The usual convention is to 
consider a proper match to be one where two pixels cover the width (w). However in 
situations constrained by, time, amount of light expected from the object to be studied, 
transmission efficiency of the optics and resources available in obtaining the desired 
optical components, the decision rests on the compromise between resolution and 
allowing for sufficient transmission of light to the detector. 


In the absence of such constraints, say (n) number of pixels are to be matched to 
the image width (w), then from equations (3.10), (3.11) and (3.15) the expression for the 
pixel width (w) is given by equation (3.17). 



In equation (3.17) the relationship (fi/di = F/D p ) is used together with the physical 
parameter for the beamwidth (di) at the collimating lens, as shown in figure (3.17). Here 
it is assumed that the grating diameter is at least (di). 
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Consider the spectrometer slit of width (dr) illuminated with light of wavelengths 
(X ) and (A. + dX ). The slit image in each of these wavelengths will have a width (w) on 
the detector. From equation (3.14) the separation between the centers of these two images 
is given by equation (3.18). 

f-in 

dl = — ^-dA, (3.18) 

dcos(oc) 

Now defining spectral purity ( 8A. ) as the wavelength difference for which 
( A1 = w ), ensuring the images are just resolved. Substituting this condition in equation 
(3.18) and also using equation (3.10) the spectral purity is given by equation (3.19). 


dcos(oc) 
oA = w 

f 2 m 

dj.dcos(|3) 


infj 


(3.19) 


where ( P = 0 ) for collimated light incident on the grating. However it needs to be noted 
that the expression for the spectral purity given by equation (3.19) cannot be smaller than 
the spectral purity limit (5?i 0 ) set by diffraction. The resolving power (R) of a grating is 
the ratio between the smallest change of wavelength that the grating can resolve and the 
wavelength at which it is operating. The resolving power is given by equation (3.20). 

R=— (3.20) 
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The spectral purity ( 8X q ) derived in equation (5.19) does not take into account 

the limit on the image size set by diffraction. What follows here is the expression for 
spectral resolution (S^ Q ) and resolving power (R 0 ) limited by diffraction. The spectral 

purity (8X) cannot be smaller than (SX q ) while the resolving power (R) cannot be 

greater than the diffraction limited resolving power (R 0 ). 



Figure (3.18). Schematic diagrams showing the diffraction limited spectral 
resolution and the resolving power. 
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From figure (3.18) the path difference (DO) between contributions from the center and 
the edge of the grating is given by equation (3.21). 


i ' 

DO = — sin(Soc) 


(3.21) 


Then the associated phase difference (£ ) is given by equation (3.22). 



(3.22) 


The intensity distribution (Iz) at the focal plane, due to diffraction, is given by equation 

(3.23). 

si » 2 (5) (3.23) 


I. =1, 




The first minimum of equation (3.23) is for (£ ) given by equation (3.24). 


2 tc 

£ = — x 
s X 

L'Z 0 

2 f 2 

= K 



(3.24) 


From equation (3.14) two wavelengths that are differing by (bX^), with each 

producing a diffraction pattern at the focal plane, as shown in figure (3.18), will be 
displaced by ( 8l 0 ) and given by equation (3.25). 


2. f 2 m SI 

8l« — / \8X 0 


(3.25) 


dcos(oc) 
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The closest the two wavelength patterns could be without merging into each other 
is given by the Rayleigh’s criterion. As depicted in figure (3.18) the two wavelength 
patterns differing by ( 8X q ) could just be resolved if the maximum of one pattern 

coincides with the first minimum of the other. In this case (81^ = Z^). Substituting (Z 0 ) 

form equation (3.24) in equation (3.25) and from figure (3.18) using ( 1 / = Lcos(a)) 
gives the following expression for (87.^) given by equation (3.26). 


SX 0 = — 

Lm 

Nm 
= ^L 


where (N) is the total number of grooves in the grating. 

The theoretical resolving power R is then given by equation (3.27). 

R 0 = Nm 


(3.26) 


(3.27) 
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The values for the physical parameters of the optical components in MACS are given in 


equation (3.28) and table (3.1). 


m = 1.0 

p = o° 

N = 3.6 xlO 4 
fq = 31.0 cm 
f 2 = 16.0 cm 
d = 1.67 x 10“ 3 mm 
d f = 200.0 pm 
pixel = 24.0 pm 
R 0 = 3.6 xlO 4 


(3.28) 


Table (3.1). A table showing the values for the angle of diffraction, the linear 
dispersion, the line width and the resolvable theoretical wavelength range for 
three different wavelengths using the physical parameters of the optical 
components used in MACS. In column 3 the reciprocal linear dispersion is 
calculated per pixel, which is the dimension of a pixel in the CCD camera 
used in MACS. 


wavelength 

o 

A 

diffraction 

angle 

a degree 

reciprocal 

linear 

dispersion 

1 ° 

A / pixel 
LD 

spectral 

purity 

8A. A 

resolvable 

theoretical 

wavelength 

range 

5X„ A 

3600.0 

12.47 

2.44 

10.77 

10.00 

4000.0 

13.86 

2.43 

10.77 

9.00 

4500.0 

15.63 

2.41 

10.77 

8.00 
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The exposure time allowed to record the spectra is dependent on the source 
brightness, the light collecting area, the transmission efficiency of the telescope- 
spectrograph system, the slit dimension and the sensitivity of the detector. Consider a 
source of brightness given by equation (3.29). 

I'; Joules/sec/cm'^/angs/steradians 13.29) 

The flux received by the Schmidt-Cassegrain telescope is then given by equation (3.30). 


F x k \ 


D 


P 


I) 


2 ) 


Joules/sec/angs/steradians 


(3.30) 


where (D p ) and (D s ) are, respectively, the diameters of the primary and the secondary 
mirrors of the telescope. 


The flux through a single fiber at the focal plane of the telescope is then given by 
equation (3.31). 


f^XJTl 




I) 


2 y 


d f d f 
F F 


Joules/sec/angs 


(3.31) 


In equation (3.31), (d^ / F ) is the angle subtended by the height and width of a single 
fiber on the sky and F the effective focal length of the telescope. 
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Let (T k ) be the wavelength dependent transmission efficiency of the telescope- 


spectrometer combination. 


The spectral flux received by a single pixel (E) is then given by equation (3.32). 


El =f\ xd 


Dp -D 


2 \ 


— — 8X — 

F F wh 


XT. Joules/sec 


(3.32) 


V J 

In equation (3.32) (w) and (h) are the projected fiber width and height on the detector and 
given by equations (3.10) and (3.11), respectively, (8X) the spectral purity and (p) the 
dimensions of a square pixel. From equation (3.32) the number photons of wavelength 
(X) and energy (h c/X) incident on a pixel per unit time (Nphotons) is then given by 

p 

equation (3.33) where (h p ) and (c) are the Planck constant and the speed of light, 
respectively. 


N 


photons 


X_ 

h c/X 

p 


sec 


-1 


(3.33) 


The basic detection mechanism of a CCD is related to the photoelectric effect 
where the light incident on the semiconductor produces electron -hole pairs. These 
electrons are then trapped in potential wells produced by numerous small electrodes. The 
maximum number of electrons that could be contained in each a pixel is called the well 
depth (WD). The quantum efficiency (QE) is the ratio of the actual number of photons 
detected to the number of photons incident at a given wavelength {X ). If the time taken 


81 



for a pixel to reach the well depth (WD) is (T) then the wavelength dependent integration 


time (T) is given by equation (3.34). 
h c/X 

T = x(WD) seconds 

E x x(QE), 


(3.34) 


The physical parameters that are used in determining the integration time (T) are given in 
equation (3.35). 


D p = 30.5 cm, 
d f =200.0 pm, 
p = 24.0 pm, 
f 2 = 16.0 cm, 

d = 1.67 x HE 3 mm, 


D s = 10.2 cm 
F = 192.0 cm 
= 31.0 cm 
m = 1.0 


a = sin 


inX 


p = o°, 

w = d f 


( e \ 




COS(P) 


cos(a) 


6X = 10.77 A 
h = d f 


V 

v r ‘ y 


c = 3.0xl0 8 ms' 1 
h = 6.63 x lO 14 Joules/sec 

p 

WD = 300,000 


QE(?i = 3600 A) = 0.32 


QE(?i = 4000 A) = 0.60 
o 

QE(?i = 4500 A) = 0.70 


(3.35) 


The transmission effciency (T x ) is a wavelength dependent parameter that 

depends on the transmission efficiency of the collimating lens, the fibers, the camera lens 
and the grating used in MACS and its wavelength dependent values are listed in table 
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(3.2). The number of optical surfaces that contribute to the transmission efficiency are 


listed in table (3.3). Table (3.4) gives effective transmission efficiency at different 
wavelength by considering the wavelength dependent transmission efficiencies given in 
table (3.2) and the number of optical surfaces through which the the light has to pass 
before entering the detector and given in table (3.3). Table (3.5) gives the K-coronal 
brightness at 1.1 and 1.6 solar radii at different wavelengths. 


Table (3.2), The transmission efficiencies of the fibers, the lenses and the 
grating for three different wavelengths. This information is from figure (3.9), 
figure (3.11) and figure (3.12). 


Wavelength 

Angstrom 

Fibers 

Lenses 

Grating 

3600.0 

0.70 

1.0 

0.33 

4000.0 

0.90 

0.99 

0.32 

4500.0 

0.92 

0.98 

0.29 


Table (3.3). The list of lenses and the associated number of surfaces through 
which the incident light had to pass before entering the detector. 


Optical element 

Number of surfaces 

Corrector plate 

2 

Focal reducer 

4 

Collimating lens 

4 

Camera lens 

4 

CCD glass protector 

2 
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Table (3.4) The effective transmission effciency ( if ) of MACS for the 
wavelengths 3600.0, 4000.0 and 4500.0 angstroms. The values are based on the 
information from tables (3.2) and (3.3). 



Table (3.5) The K-coronal brightness at 1.1 and 1.6 solar radii for three 
different wavelengths during the solar maximum. The values were obtained 
from Allen (1973) from pages 172 and 176. The units are in 
ergs/sec/cm 2 /angs/steradians. (1 Joule=10 7 ergs) 










Substituting equation (3.32) in equation (3.34) the expression for the integration time (T) 


is given by equation (3.36). 


(h cA)x(WD) 

p 


L X Til 


V 


D 


2 \ 


seconds 


^-S^xt' , 'x(QE ) 1 

F F wh x 


(3.36) 


Using the values for the physical parameters from equation (3.35) and tables (3.4) and 
(3.5) gives the following chart, as shown in table (3.6), for the lower limit of the 
integration time (T) for the K-corona at 1.1 and 1.6 solar radii. 

Table (3.6). Lower limit of the integration time for the K-coronal spectrum for 
MACS based on equation (3.36). 


Wavelength 

T in seconds 

T in seconds 

angstrom 

at 1.1 solar radii 

at 1.6 solar radii 

3600.0 

122.0 

177.0 

4000.0 

38.0 

55.0 

4500.0 

27.0 

39.0 
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However the actual integration time will depend on the amount of light absorbed 
by the various optical elements, the amount of light escaping the light collecting area of 
the optical elements and further contributions from the F-corona. 

As for the spatial resolution of MACS the width of a single fiber of 200.0 pm and 
a plate scale of 0.107 arcseconds/pm corresponds to a spatial scale of 21.4" (0.022 solar 
radii). Taking the height of a single fiber as 200.0 pm and from equation (3.11) the 
projected height of the fiber on the spatial direction of the detector is given by equation 
(3.37). 


h = 200pm x 


16.0 cm 


31.0 cm 


= 103.23pm 


(3.37) 


This corresponds to (103.23pm /24.0pm) 4.30 pixels in the spatial direction. Therefore 
form the spatial scale per fiber and the projected height of the fiber given in equation 
(3.37) the spatial resolution is given by equation (3.38). 


Spatial resolution = : 

4.30 pixels 

= 4.98"/ pixel 
1 

~ 5.19 x 10'" solar radii / pixel 


(3.38) 
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3.4 Cost of the instrument 


Table (3.7). Cost associated with the construction of MACS. 


Item 

$000 

Manufacturer 

Telescope 

6.79 

Meade Corporation 

Solar filters 

1.02 

Thousand Oak Optical 

Fibers 

4.05 

Oriel Instruments 

Collimating lens 

1.39 

Spindler & Hoyer 

Camera lens 

1.58 

Spindler & Hoyer 

Grating 

2.35 

American Holographic 

CCD camera 

10.5 

Apogee Instruments 

Camera software 

0.60 

PMIS 

Computer 

2.90 

DPS 9000 

GPS 

0.41 

Eagle Map Guide Pro 

Glass base 

0.45 

United Glass Company 

Spectrograph body 

7.25 

Lorr Company 

Generator 

0.77 

Honda 

Total 

40.31 
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Chapter 4 


OBSERVATIONAL RESULTS AND ANALYSES 
4.1 Observational site 

The coronal intensity measurements on the total solar eclipse of 11 August 1999 
were conducted by Lake Hazar in Elazig, Turkey. The totality at this location 
(38°41 , N,039°14 , E) lasted for 124 seconds. The second contact was at ~ 11.36 GMT. 
Also no clouds could be reported for the entire duration of the totality. 

Figure (4.1) is a photograph of the total solar eclipse of 11 August 1999 taken 
from the Black Sea area. From figure (4.1) it is evident that the coronal brightness was 
almost uniformly distributed all around the sun. This was expected of the sun that was 
approaching maximum phase in 2000. This photograph also reveals a detached 
prominence in the southwest limb during the totality. Figure (4.2) depicts the path of the 
eclipse where the maximum duration of the totality occurred in Romania. Figure (4.3) 
shows the layout of MACS. The spectrograph lies on the metal box. The telescope, 
CCD-camera, and the computer were powered from a portable generator. 
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Figure (4.1). A photograph of the 
August 1999 eclipse from the Black 
Sea area. The photograph also reveals 
a detached prominence eruption above 

the southwest limb. 


Total Solar Bclipsc of 1999 August 1 1 

FluiiRK 12: Till EcUPSK P.VI1I TiiROlOH Tl'KKEV. SVRt\ wo IRAQ 



Figure (4.2). A map showing the path of the 
eclipse over Turkey. 

Courtesy: NASA Reference Publication 1398 
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4.2 Spectra observed by MACS 

Figure (4.4) shows the spectra due to the twenty-one fibers being exposed to the 
sky and the remaining four fibers being exposed to the mercury calibration lamp in 
Elazig, Turkey. The horizontal lines indicate the prominent absorption lines in the shy 
spectrum and the four bright dots on the two ends correspond to the emission line at 
4046.7 angstrom in the mercury calibration lamp. This information was used to determine 
the wavelength scale and was determined to have a linear dispersion of 2.39 
angstrom/pixel. The positions numbered from 1 to 4 are the locations corresponding to 


Figure (4.3). A picture showing the layout of MACS and the telescope. 
The telescope is on a polar mount and MACS is laid out on the metal box. 
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the temperature and the wind sensitive wavelength positions, which is explained in 


chapter-2. 





0 100 £00 500 400 

Pr<*i (ttimtver - Spat® director; (a) 


Figure (4.4). The spectra from twenty-one fibers exposed to the sky 
while four other fibers exposed to the mercury calibration lamp in 
Elazig, Turkey. 


MACS was exposed to the coronal light during totality for duration of 4.0, 18.0 
and 83.0-seconds. Figures (4.5 a), (4.5 b) and (4.5 c) show the spectra recorded by 
MACS for the above duration, respectively. Each vertical line corresponds to the 
spectrum from a single fiber with a total of twenty-one vertical lines corresponding to the 
twenty-one fibers located in the focal plane of the telescope. 
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Figure (4.5 a, b, and c). The coronal spectra recorded by 
MACS for the exposure times 4.0, 8.0 and 83.0 seconds, 
respectively. Each vertical line corresponds to a spectrum 
from a single fiber. It is also clear that five of the fibers 
show emission lines. 


4.3 Identifying the individual fiber location on the corona 

Prior to the eclipse the sun was focused to the position identified by the moon 
during eclipse, as shown in figure (4.6). A solar filter attached at the front end of the 
telescope facilitated this process. The telescope was then allowed to track. Watching the 
direction of movement of the sun’s image upon slewing the telescope due north and south 
with the auto controls the fiber #19 was oriented at the position perceived to be the 
earth’s South Pole. Its relation to the heliographic coordinates was determined from the 
value listed for the P-angle for that day in, ‘The Astronomical Almanac 1999’. 
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Figure (4.6) shows the envisaged location of the twenty-one fibers in the focal 


plane of the telescope during the eclipse. However the envisaged location could have 
drifted due to imperfect tracking or minute movements of the whole telescope structure 
for not being anchored on hard ground. Time constraints and the remoteness of the 
location compounded these problems. The term “envisaged” is emphasized because of 
the incapability of MACS at the present state to simultaneously record the image. This 
capability would have verified the exact location of the fiber rings during the eclipse. 



Figure (4.6). This picture shows the envisaged location of the twenty-one 
fibers in the focal plane during the eclipse. The inner and the outer rings 
correspond to 1.1 and 1.5 solar radii, respectively. The picture is not the 
actual eclipse on 11 August 1999. The fiber in the center was expected to 
record the background signal. 
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Ill order to identify the fiber orientation to a reasonable accuracy, some other 
information shown in the coronal spectra was utilized. This task was further aided by the 
very distinctive emission lines recorded in the spectra of five of the fibers as revealed in 
figure (4.5 a), figure (4.5 b), and figure (4.5 c). These emission lines are most pertinent to 
be from fibers exposed to prominence eruptions. Using imaging data at totality from 
other space and ground based solar instruments that are dedicated to observe specific 
solar features, it was now possible to infer the solar features that would have contributed 
to the emission lines in those fibers and accordingly orient the fibers. 

The fibers that show emission lines, going from left to right in figure (4.5 a), 
figure (4.5 b), and figure (4.5 c), are the fibers numbered 22, 20, 19, 15 and 14 in figure 
(4.6). All the emission lines of fiber #19 correspond to the signatures of a prominence. 
Table (4.1) gives a list of prominence lines reported by Athay and Orrall (1957) from a 
prominence in the west limb of the total solar eclipse of 1952. Figure (4.7) shows the 
spectrum recorded by fiber #19. In this paper an emission line of Sr II listed at 3077.7 
angstrom needs to be corrected as 4077.7angstrom. Comparison between table (4.1) and 
figure (4.7) reveals that the fiber #19 had recorded all the prominence emission lines 
reported by Athay and Orrall (1957), that falls within the wavelength range of the 
detector in MACS. 
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Table (4.1) List of the emission lines reported by Athay and Orrall 
(1957) of a prominence from the total solar eclipse of 1952. In this 
paper the Sr II line at 4077.7 has been erroneously recorded as 
3077.7 angstrom. In the following table only the emission lines that 
could lie within the wavelength range of MACS is enumerated. 


Emission line 

Wavelength (angstrom) 

Sr II 

4215.5 

H 6 (2-6) 

4101.7 

Sr II 

4077.7 

He I 

4026.3 

Ca II H 

3968.5 

Ca II K 

3933.7 

H 8 (2-8) 

3889.1 

H 9 (2-9) 

3835.4 


Figure (4.7) shows the spectrum recorded by the fiber #19 in the 4.0 seconds 
exposure. This is after averaging over three pixels about its center and subtracting the 
dark signal. However the spectrum is not yet corrected for the wavelength sensitivity of 
the telescope-spectrometer system. 
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Pixel number - spectral direction 


Figure (4.7). The spectrum recorded by the fiber #19 in the 4.0 
seconds exposure. The emission lines and their corresponding 
wavelengths and pixel positions are identified. 


From figure (4.7) and table (4.1) it is apparent that the fiber #19 had recorded all 
the prominence lines that fall within the wavelength range of the detector. The yellow 
curves superimposed on the spectrum are the gaussian curve fits. The wavelengths 
correspond to the peaks of the gaussian curves. However it needs to be yet ascertained 
whether the emission lines were correctly identified. This was verified by constructing a 
regression fit between the emission peaks and the peak positions. Figure (4.8) shows the 
result of the regression fit and the wavelength scale between different emission lines. 
These values give an average linear dispersion per pixel of 2.39 angstrom and matches 
with the linear dispersion per pixel calculated from the optical parameters of MACS. 
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Figure (4.8). The regression fit between the emission line peaks 
and the peak positions of the fiber #19 in the 4.0 seconds 
exposure. The table incorporated shows the wavelength scale 
between any two-emission lines. 


Furthermore, the spectrum recorded by the fiber #20 also shows all the emission 
lines from a prominence as listed in table (4.1). Figure (4.9) is the spectrum recorded by 
the fiber #20. 
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Figure (4.9). The spectrum recorded by the fiber #20 in the 4.0 
seconds exposure. This too matches with all the emission lines 
recorded by the fiber #19 as shown in figure (4.7). 


In table (4.2) is a list of emission lines, categorized as very strong and strong, 
from the flare of 18 September 1957 and reported by Jefferies, Smith and Smith (1959). 
Figure (4.10) shows the spectrum recorded by the fiber #22 hi the 4.0 seconds exposure. 
This is after averaging over three pixels about its center and subtracting the dark signal. 
However the spectrum is not yet corrected for the wavelength sensitivity of the telescope- 
spectrometer system. 
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Table (4.2). List of the emission lines reported by Jefferies, Smith 
and Smith (1959) from the flare of 18 September 1957. This table 
only includes those emission lines categorized as very strong and 
strong. Also in this following table only the emission lines that 
could lie within the wavelength range of MACS is enumerated. 


Emission line 

Wavelength 

angstrom 

Category 

h 8 

3889.1 

strong 

Ca II- K 

3933.7 

very strong 

Ca II-H 

3968.5 

very strong 

H £ 

3970.1 

strong 

H S 

4101.7 

very strong 

H y 

4340.5 

very strong 

Hel 

4471.7 

strong 


Figure (4.10) showing the spectrum recorded by the fiber #20 does indicate that 
all but the H g line at 3970.1 angstrom listed in table (4.1) to be present. Since the linear 

dispersion is ~2.40 angstrom/pixel and, with the H e line at 3970.1 angstrom being only 

2.0 angstrom away from the very strong Ca II-H line at 3968.5 angstrom, this may have 
caused this line to blend with the stronger line. 
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Figure (4.10). The spectrum recorded by the fiber #22 in the 4.0 
seconds exposure. The emission lines and their corresponding 
wavelengths and pixel positions are identified. 


Figure ( 4 . 11 ) shows the result of the regression fit and the wavelength scale 
between different emission lines of figure ( 4 . 10 ). These values give an average linear 
dispersion per pixel of 2.40 and matches with the linear dispersion per pixel calculated 
from the optical parameters for MACS, as listed in table ( 3 . 1 ). 
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Figure (4.11). The regression fit between the emission line peaks 
and the peak positions of the fiber #19 in the 4.0 seconds 
exposure. The table incorporated shows the wavelength scale 
between any two-emission lines. 


Based on the comparison on the spectra recorded by the fiber#19, fiber#20 and 
fiber #22 and the prominent emission lines expected of a prominence or a flare as listed 
in tables (4.1) and (4.2), respectively, it is reasonable to assume that these fibers were 
exposed to a prominence or a flare or to both of these features. 
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Ill order to identify the individual fiber location on the corona as accurately as 

possible the following analyses were performed. 

1. Identifying the fibers exposed to prominence spectrum. The fibers #19, 20 and 22 
were positively identified to have recorded the typical emission lines due to a 
prominence. 

2 . Orienting an eclipse photograph showing prominent eruptions with EIT images from 
SOHO on the day of the eclipse and matching these features. This matching was 
performed using the conspicuous prominence eruption identified at ~ 68° southwest 
in the EIT picture in 304 angstrom by SOHO at 1 1:26 UT with an eclipse photograph 
taken at -11:16 UT from the Black Sea. 

3 . Drawing the polar axis of the earth, which is p-angle from the polar axis of the sun. 
The p-angle, which is the orientation of the polar axis of the sun with respect to the 
polar axis of the earth, was obtained from, ‘The Astronomical Almanac 1999’. 

4 . Placing fiber #19 at the intersection point of the South Pole with the orbital circle at 
1.5 solar radii. Prior to the eclipse the sun was focused to the position identified by 
the moon during eclipse, as shown in figure ( 4 . 6 ). A solar filter attached at the front 
end of the telescope facilitated this process. The telescope was then allowed to track. 
Watching the direction of movement of the sun’s image upon slewing the telescope 
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due north and south with the auto controls the fiber #19 was oriented at the position 


perceived to be the earth’s South Pole. The locations of the other fibers could be then 
determined based on their positions with respect to fiber #19. 


Figure (4.12) shows the envisaged positions of the fibers on the corona based upon the 
above procedure. 



Figure (4.12). This is a comparison between solar image by SOHO in 
304.0 angstrom at 11:14 UT and an eclipse photograph taken from the 
Black Sea at ~ 11:16UT. The orientation was made by matching the 
prominence eruption ~ 68° southwest in the SOHO image. The orientation 
of the polar axis of the sun with respect to the polar axis of the earth was 
obtained from the p-angle data given in, ‘The Astronomical Almanac 
1999’ for 11 August 1999. 
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Figure (4.13), figure (4.14) and figure (4.15) are pictures of the sun taken by 


LASCO instrument on SOHO, EIT instrument on SOHO and radio image from Nancy 
and x-ray image by Yohkoh, respectively, on 1 1 August 1999. 



A 1 '• 


1999/03/11 12:08 


lasco C2 12:08 


lasco C2 12:40 


lasco C3 12:15 


lasco C3 21:42 


Figure (4.13). Images of the sun taken by the LASCO instrument 
on SOHO on the 11 August 1999. The numbers denote the time in 
UT. C2 and C3 are coronal images of the sun from 2.0 - 6.0 solar 
radii and 3.7 - 32.0 solar radii, respectively, using an externally 
occulted disk. 
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Figure (4.14). Images of the sun taken by the EIT instrument 
on SOHO on the 11 August 1999. The numbers denote the 
time in UT. 
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YOHKOH 14:42 UT 164 MHz 10:45 UT 327 MHz 10:45 UT 


Figure (4.15). Images of the sun taken by the radiograph in Nancy, France 
and the x-ray telescope on Yohkoh on the 11 August 1999. 


In figure (4.13) the coronal images of the sun were obtained from the C2 and C3 
coronagraphs, of the LASCO (Large Angle and Spectrometric Coronagraph) instrument 
on SOHO, providing coronal images of the sun from 2.0 - 6.0 solar radii and 3.7 - 32.0 
solar radii, respectively. These pictures reveal large-scale activities in the form of 
streamers in the southeast limb as compared to the rest of the limb. From figure (4.12) the 
fibers #19, 20, 21 and 22 were located in this region of high streamer activity. 

In figure (4.14) are EIT (Extreme ultraviolet Imaging Telescope) on SOHO are 
images of the sun through various filters centered at different wavelengths as indicated on 

o 

the images. These are images, through selected bandpasses in Fe I X/X at 171.0 A, Fe 

XII at 195.0 A , Fe XV at 284.0 A and He II at 304.0 A , reflect temperatures of ~1.0 
MK, ~2.0 MK, ~3.0 MK and ~4.0 MK, respectively. 
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Ill figure (4.15) are radio images from Nancy, France and a x-ray image from 
Yohkoh. The radio image in the center reveals that the east limb was either hotter or 
lesser in electron number density than the west limb. The bright spot on the image on the 

right shows the subflare #8662 reported at southl4 east04 on 1 1 August 1999. 

4.4 The telescope-spectrometer sensitivity curve 

In section 4.4 the wavelength correction factor to convert terrestrial intensity 
observation to extraterrestrial intensity is discussed. The intensity recorded by MACS 
differs from the extraterrestrial intensity primarily due to the following two reasons. 

1. The wavelength dependency of the physical parameters such as transmission 
efficiency, focal length, dispersion and quantum efficiency, etc of the optical 
components of MACS. 

2. The wavelength dependency of the measured coronal intensity on the transmission 
efficiency and the scattering in the earth’s atmosphere. 

In order to find a correction term to convert terrestrial intensity observation to 
extraterrestrial intensity for each of the fibers the following steps were followed. This 
follows the procedure adopted by Ichimoto et al. (1996). 
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1. The dark-current subtracted terrestrial sky spectrum observed by each fiber was 
corrected for Rayleigh scattering ( X” 4 ). The sky represents the F-component of the 
solar spectrum and also closely resembles the solar spectrum itself. However the 
Rayleigh scattering affects the overall distribution of the sky spectrum. Thus the 
Rayleigh corrected sky spectrum was assumed to closely represent the terrestrial solar 
spectrum at the observation site. Figure (4.16) shows the sky spectrum recorded by 
fiber #06 before and after correcting for Rayleigh scattering. 


Sky spectrum recorded by fiber $06 in Elazig corrected for Rayleigh scattering 



Figure (4.16). Plots of the sky spectrum recorded by fiber #06 in Elazig, 
Turkey before and after correcting for Rayleigh scattering. The plots 
have been normalized at 3968.5 Angstrom. 
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2. The terrestrial photospheric spectrum (Rayleigh corrected sky spectrum) recorded by 
MACS was calibrated with a wavelength scale using prominent absorption lines Ca II 
K & H at 3933.7 & 3968.5 angstrom, respectively and H 3 at 4101.7 angstrom. The 

positions of these absorption lines were matched with the corresponding positions of 
the same absoiption lines in the extraterrestrial photospheric spectrum. These 
intensities were then divided by their respective intensities at 3968.5 angstrom. Figure 
(4.17) shows the terrestrial photospheric spectrum (Rayleigh corrected sky spectrum) 

recorded by fiber #06 and the extraterrestrial photospheric spectrum. 


Elaziq sky spectrum compared with the extraterrestrial photospheric spectrum 



Figure (4.17). Plots of the extraterrestrial photospheric 
intensity spectrum and the terrestrial photospheric 
intensity spectrum recorded by fiber #06. The 
extraterrestrial photospheric spectrum is a smoothed 
version of the extraterrestrial photospheric spectrum 
shown in figure (2.1). 
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3. The extraterrestrial photospheric spectrum was smoothed until it matched the 
resolution of the terrestrial photospheric spectrum recorded by MACS. Figure (4.18) 
shows the terrestrial photospheric spectrum recorded by fiber #06 and the 
extraterrestrial photospheric spectrum smoothed in resolution to match the resolution 
of the terrestrial photospheric spectrum. 


Elazig sky spectrum compared with the smoothed extraterrestrial photospheric spectrum 



Figure (4.18). Plots of the smoothed extraterrestrial 
photospheric intensity spectrum and the terrestrial 
photospheric intensity spectrum recorded by fiber #06. 
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4. 


The ratios between the extraterrestrial and the terrestrial photospheric spectra 


were calculated against wavelength. Figure (4.19) shows plots of the smoothed 
extraterrestrial photospheric intensity spectrum, terrestrial photospheric intensity 
spectrum recorded by fiber #06 and the ratio between them. 


Smoothed extraterrestrial photospheric spectrum and the Rayleigh corrected sky spectrum 



Figure (4.19). Plots of the smoothed extraterrestrial 
photospheric intensity spectrum, terrestrial photospheric 
intensity spectrum recorded by fiber #06 and the ratio between 
these two against wavelength. 
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5 . 


A polynomial fit was made to the ratios between the extraterrestrial and the 


terrestrial photospheric spectra that were calculated against wavelength. Figure 
( 4 . 20 ) shows the polynomial fit made to the ratio between the smoothed 
extraterrestrial photospheric intensity spectrum, and the terrestrial photospheric 
intensity spectrum recorded by fiber # 06 . 


Wavelength dependent intensity correction factor 
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Figure (4.20). Polynomial fit made to the ratio 
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Figure (4.21) shows the application of the polynomial fit to the terrestrial 
photospheric intensity spectrum to obtain the extraterrestrial photospheric intensity. This 
indicates a reasonably good fit. Also plotted in figure (4.12), for comparison purpose, is 
the smoothed extraterrestrial photospheric intensity spectrum as shown in figure (4.18). 


Fitting Elazig sky spectrum to the extraterrestrial photospheric spectrum 



Figure (4.21). Application of the polynomial fit shown in figure 
(4.20) to the terrestrial photospheric intensity spectrum to 
obtain the extraterrestrial photospheric intensity. Also plotted 
for comparison purpose is the smoothed extraterrestrial 
photospheric intensity spectrum as shown in figure (4.18). 
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4.5 Application of the intensity correction factor 


The extraterrestrial photospheric intensity spectrum differs from the terrestrial 
photospheric intensity spectrum measured by MACS primarily due to the atmosphere and 
the optical elements in MACS. Therefore the ratio between the extraterrestrial and the 
terrestrial photospheric intensity at a given wavelength could be written as shown in 
equation (4.1). 


IpT (^V IpT M 

I EX PO/IeX M 


= f(A,x) = f(X.) 


x(X) = transmission efficiency of MACS 
A(A,) = atmospheric attenuation 


(4.1) 


where I PT and I EX are the photospheric intensity measured by MACS from ground and 

the extraterrestrial photospheric intensity, respectively. Here both the extraterrestrial and 
the terrestrial intensities are divided by their respective intensities at a reference 
wavelength X 0 . In producing figure (4.18) the reference wavelength X 0 was selected to 

be the Ca-II H line at 3968.5 angstrom. The same is true for the measurement of any 
other extraterrestrial intensity spectrum including the measurement of the solar coronal 
intensity measurement during an eclipse. 
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Using equation (4.1) the terrestrial coronal intensity at a given wavelength^ 
measured by MACS during the eclipse would correspond to the extraterrestrial coronal 
intensity given by equation (4.2). 




Iex(^o) 

f(*T igM 


xigW 


q (X) = terrestrial coronal intensity 
ex (X) = extraterrestrial coronal intensity 


(4.2) 


From equation (4.2) the extraterrestrial coronal intensity ratio at wavelengths A., and X 2 
is given by equation (4.3). 


I EX (^l) ^2) v (^1 ) ( . 

ig(M ( ,3) 

4.6 Determination of the Thermal electron temperature and the Solar wind 
velocity for each fiber location 

The coronal spectrum recorded by each fiber contains the K, F and E coronal 
spectra. In order to determine the thermal electron temperature and the solar wind 
velocity, as per the methodology discussed in chapter-2, the K-coronal spectrum has to be 
isolated. Figure (4.22) shows the coronal spectrum recorded by the twenty-one fibers in 
the 18.0 seconds exposure. 
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Figure (4.22). These plots show the coronal spectra recorded 
by the twenty-one fibers located at various latitudes and 
heights on the corona as depicted in figure (4.12). Fiber #23 
was centered on the moon. The above spectra were recorded in 
the 18.0 seconds exposure. The x-coordinate is in pixel units. 


117 










Comparison between figure ( 4 . 22 ) and ( 4 . 12 ) it is evident that the fibers below the 
equator prominently display emission lines. From fiber #23, which was located at the 
center of the moon, was expected to record the background counts in addition to the 
coronal light reflected off the clouds and re-reflected off the moon. However the 
observational site in Elazig, Turkey reportedly being free of clouds the spectrum from 
fiber #23 only reveals the extent of the instrumental stray light. This stray light is 
believed, largely, to have originated in the telescope and could pose a significant error in 
the wind measurement. This is because the scattering is significantly large in the bright 
end of the spectrum in MACS. This in turn could give a larger count number for the 
intensity at 4233.0 angstrom, which is used in determining the wind sensitive intensity 
ratio. This is evident in the spectrum of fiber #23. In this regard it a prudent to eliminate 
at least l-% of the counts in the bright end of the spectrum. The fibers recording emission 
lines can be attributed to exposure to prominence spectrum as discussed in section ( 4 . 3 ). 

Figure ( 4 . 23 ) is the coronal spectrum recorded by fiber #06. As per figure ( 4 . 12 ) 
this fiber was located midway between the solar equator and the solar North Pole at 1.1 
solar radii. Also LASCO images in figure ( 4 . 13 ) confirm this region free of streamer 
activities. The following steps were carried out in order to isolate the K-coronal spectrum 
in fiber #06 and closely resemble the methodology described by Ichimoto et.al (1996). 
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1 . 


It was assumed that the Rayleigh corrected sky spectrum to closely resemble the 


F-component of the solar spectrum. On this basis a fraction of the sky component 
corresponding to about l-% of the counts at 4250.0 angstrom in the coronal 
spectrum was subtracted. 


Coronal spectrum recorded by fiber #06 



Figure (4.23). The coronal spectrum recorded by fiber #06. 


2. A third order polynomial fit was made to the difference between the coronal spectrum 
substracted by the F-coronal spectrum. This was possible in the above case since the 
coronal spectrum of fiber #06 is devoid of any significant emission lines, hi the cases 
with strong emission lines gaussian fits would be made and subtracted with the 
background added to the continuum. In order to account for the stray light a constant 
value was subtracted from the coronal spectrum in order to make the overall variation 
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with wavelength (the ratio of the average intensity in 3700-3800 angstrom and the 
intensity around 4250 angstrom) to be identical with that of the brightest spectrum 
recorded by a fiber. For fiber #06 this correction amounted to 55 digital units in a 16- 
bit dynamic range. Figure (4.24) shows the reduced K-coronal spectrum from fiber 
#06. 


Terrestrial K— coronal spectrum recorded by fiber $06 



Figure (4.24). The terrestrial K-coronal spectrum reduced from 
the coronal spectrum recorded by fiber #06. 


3. The terrestrial K-coronal spectrum was adjusted by the wavelength sensitivity 
curve in figure (4.20). This would transform the terrestrial K-coronal spectrum to 
the extraterrestrial K-coronal spectrum. This curve was used to determine the 
wind and the temperature sensitive ratios using equation (4.2). Figure (4.25) 
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shows the extraterrestrial K-coronal intensity derived for fiber #06 together with 
the wind and the temperature sensitive intensity ratios. 


Extraterrestrial K— coronal spectrum calculated for fiber $06 



Figure (4.25). This plot shows the extraterrestrial K- 
coroual spectrum derived for fiber #06. The wind and the 
temperature sensitive intensity ratios are 1.245 and 
1.488, respectively. 

4. The theoretical temperature and the wind sensitive plots for the desired line of sight 
were then used to predict the thermal electron temperature and the solar wind 
velocity. The figure (4.26) and figure (4.27) show the locations of the temperature 
and the wind sensitive intensity ratios, derived from figure (4.25), on the theoretical 
thermal electron temperature and the solar wind velocity plots at 1.1 solar radii, 
respectively. These plots associate a thermal electron temperature of -1.73 MK and a 
solar wind velocity of -467.0 km/sec for fiber #06’s location on the solar corona. The 
error bars will be discussed in section (4.7). 


121 





1.70 


1(41 00}/i(3SBO) VS Temperature 



Temperature Mk 


rho in solar radius 1.1000000 



0 200 400 600 800 1000 


'Mind km/sec 


Figure (4.27). The wind sensitive plot. This assigns 
a solar wind velocity of ~467.0 kni/sec for fiber 


#06. 
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Similarly figure (4.28) and figure (4.29) show the locations of the temperature and 
the wind sensitive intensity ratios on the theoretical thermal electron temperature and the 
solar wind velocity plots at 1 . 1 solar radii, respectively. These plots associate a thermal 
electron temperature of ~1.29 MK and a solar wind velocity of ~300.0 knt/sec for fiber 
#06 ’s location on the solar corona. Figure (4.30) shows the extraterrestrial K-coronal 
spectrum derived for fiber #10. Figure (4.31) and figure (4.32) give a thermal electron 
temperature of ~1.34 MK and a solar wind velocity of ~571.0 knt/sec for fiber #10’s 
location on the solar corona. 


Extraterrestrial K— coronal spectrum calculated for fiber #10 



Figure (4.30). This plot shows the extraterrestrial K- 
coronal spectrum for fiber #10. The wind and the 
temperature sensitive intensity ratios are 1.248 and 
1.524, respectively. 
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1(4-1 GO)/i( 3850) VS Temperature 




Figure (4.32). The wind sensitive plot. This assigns 
a solar wind velocity of '571.0 kni/sec for fiber 
# 10 . 
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Figure (4.33) shows the coronal spectrum recorded by fiber #14 that includes 
prominent emission lines in calsium and hydrogen. The continuum was obtained by 
fitting a gaussian function of the form given by equation (4.4) to the emission lines and 
fitting equation (4.5) to the continuum. 


f(?i) = AqC " + A ^ + A jA, + 


where 


Z = 


f*-A, 


\ 



A ^ + A^X + A 


(4.4) 


(4.5) 


Coronal spectrum recorded by Fiber $14- in the 4— seconds exposure 



Figure (4.33). This plot shows the terrestrial coronal 
spectrum recorded by fiber #14. The continuum is a 
measure of the terrestrial F+K corona. The E-corona is 
removed by making Gaussian fits of the form given in 
equation (4.1) to the individual emission lines and fitting 
the background given by equation (4.2) to the continuum. 
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Following the steps described above figure (4.34) shows the extraterrestrial K- 


coronal spectrum reduced for fiber #14. 


Extraterrestrial K— coronal spectrum calculated for fiber $14 



Figure (4.34). This plot shows the extraterrestrial K- 
coronal spectrum for fiber #14. This is following the 
elimination of the E-corona from figure (4.33). 


From figure (4.35) and figure (4.36) the thermal electron temperature and the 
solar wind velocity derived for fiber #14 are ~ 0.73 MK and ~ 819.0 km/sec, respectively. 
Flowever from figure (4.12) and figure (4.13) it is evident of a streamer and some bright 
structures in the field of view of fiber #14. This fairly low temperature could result from 
one of two possibilities. Firstly the fiber may have seen the corona very close to the limb 
based on the possibilities of errors arising from improper alignment of the fiber optic 
plate in the focal plane of the telescope or improper tracking. 
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Secondly, it is apparent that local density enhancement may have resulted at 
points of intersection of the streamer with the line of sight of fiber #14. Although the 
actual geometry nor the magnitude of the density enhancement are known, nevertheless, 
using a trial geometry and an assumed magnitude for density enhancement it can be 
shown that the temperature determined for fiber #14 may differ widely from the true 
value. This phenomemon is described in Appendix-C. In the following example the 
geometry described in Appendix-C was used along with an assumed density 
enhancement by a factor of 50.0. 


I(+100)/I(3850) VS Streamer angle 



angle in degrees 


Figure (4.37). Plot showing the variation in the 
temperature-sensitive intensity ratio, from the 
theoretical models without streamers, for a streamer of 
geometry described in Appendix-C and density 
enhanced by a factor of 50.0 together with various 
streamer inclinations in front and behind the solar limb. 
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Ill figure (4.37) the horozontal lines show the temperature-sensitive intensity ratio 
for models without streamers. The curves with the matching colors show the variation of 
this ratio with streamers of the above geometry and density enhancement located at 
various angles in front and behind the limb. The two arrows indicate that a corona with 
temperature 1.0 MK and wind velocity 819.0 km/sec along with a streamer enhanced in 
density by a factor 50.0 and inclined at ~ 18.0° or ~58.0° degrees in front of the limb 
could have given the same temperature- sensitive intensity ratio as the theoretical model 
without streamers for an isothermal corona of temperature 0.73 MK and solar wind 
velocity 819.0 km/sec. This example shows the sensitivity of the streamers crossing the 
lines of sight of the fibers in determining the temeprature and the wind velocity. However 
this cannot be considered a detrimental factor in the broader science of determining the 
density enhancement associated with streamers, their temperatures and velocities. 

My belief for the low temperature recorded by fiber #14 has the following 
interpretation. From figure (4.1) and (4.14) it is apparent that the fiber #14 was closely 
placed to the extended prominence. It is very likely in this this scenario that the line of 
sight integration for fiber #14 contained cool material associated with a prominence and 
is a signature of its likely temperature. 

From similar analyses the schematic diagram in figure (4.38) shows the values 
determined for the thermal electron temperature and the solar wind velocity for fibers 
#04, 06, 08, 10 ,12, 14 & 22 located at 1.1 solar radii and fiber #15 & 17 located at 1.5 
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solar radii. Unfortunately it has not been successful in determining the wind velocities for 


fibers #15 and 22. The spectrum recorded by these fibers show a gradual decline in 
counts above 4200.0 angstrom for reasons, which are not obvious. All the temperature 
and wind values determined for the above fibers are based on the theoretical models 
without consideration for streamers. 



Figure (4.38). The thermal electron temperature and the solar wind 
velocity determined for solar coronal locations observed by fibers #04, 06, 
08, 10, 12, 14 & 22 that are located at 1.1 solar radii and fiber #15 & 17 
located at 1.5 solar radii. These values are based on the theoretical 
models. The temperature and the wind velocity are measured in MK and 
km/sec, respectively. However the wind measurements seem too far high. 
In this regard it is acknowledged that the experiment was not successful in 
measuring the wind velocities and the absolute values should be 
disregarded. 
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4.7 Sources of error and error bars 


The sources of systematic errors could have resulted from improper alignment of 
the fiber optic plate in the focal plane of the telescope and tracking inadequecies. These 
practical reasons for these problems in the present status of MACS are highlighted in 
chapter-7 with suggested remedial measures. 


Another disappointment in this experiment reflects from the spectrum recorded by 
finer #23, which was centered in the center of the lunar shadow to record the background 
counts. The spectrum recorded by fiber #23, which was intended to measure the 
background counts, is shown in figure (4.39). 





Figure (4.39). The spectrum recorded 
by fiber #23, which was centered on 
the lunar shadow during the eclipse. 
Its purpose was to measure the 
background counts. 
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Figure ( 4 . 39 ) highlights the scale of the instrumental scattering specially in the high end 
of the wavelength region that is sensitive to the wind-sensitive intensity ratio giving a 
higher count at 4233.0 angstrom. An attempt to quantify its influence on the spectra 
recorded by the other fibers is given in chapter-6. This may be the reason for obtaining 
very large measurements for the wind velocities. The wind velocities expected at 1.1 
solar radii are of the order of ~100.0 km/sec. In this regard the absolute values 
determined for the wind velocities should be disregarded. However the ratios of the wind 
velocities in the equatorial regions to the coronal hole regions do make sense where the 
wind velocities at the poles are of the order ~3.0 times the wind velocities in the 
equatorial regions. Instrumental scattering could also be responsible for causing the 
coronal light recorded by a number of fibers unusable for analyses. 

However scattering is an inevitable consequence arising from the financial 
constraints associated with obtaining superior optics. Also in respect to the time 
constraint associated with the eclipse duration it was not possible to do repeat 
measurements for random statistics. However this effort will be feasible with MACS 
attached to a coronagraph. 

As for the temperature-sensitive intensity ratio one of the concern was the 
assocition of the wavelength positions to be in close proximity to strong emission lines. 
For example the proximity of the 4100 angstrom wavelength position to the emission line 
of Hg at 4101.7 angstrom. 
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In addition there is also an error associated with the derivation of the instrumental 


wavelength sensitivity curve. Figure (4.21) shows the application of the wavelength 
sensitivity curve to the Rayleigh corrected sky spectrum to obtain the extraterrestrial 
photospheric spectrum (red curve). Matching this extraterrestrial photospheric spectrum 
(red curve) with what is considered as the true extraterrestrail photospheric spectrum 
(green curve), as depicted in figure (4.21), the overall deficiency is 0 . 44 %. This is the 
difference between the area under these two curves. However this does not mean that 
point to point differences are las low as 0 . 44 %. At the temperature-sensitive wavelength 
locations at 3850.0 and 4100.0 angstroms, these errors are - 4.1 % and + 1.9 %, 
respectively. This can cause the temperature to be under estimated by 5.7 %, which 
translates to 1.14 MK. The ideal remedy to this problem is an absolute wavelength 
calibration of MACS with sources of known strengths. However this does not still 
account for the instrumental scattering. 

Although no clouds were reported during observation, yet the haze due to 
convection that is associated with the very dry terrain in Elazig, Turkey is another hidden 
consequence on the accurate coronal brightness measurements. 

The influence of the numerical procedure and other physical parameters, such as 
the electron number density profile, used in creating the shape of the K-coronal intensity 
spectrum is discussed in chapter-5. 
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4.8 Summary 

It has been demonstrated that the concept of MACS could be used as an 
instrument for global and simultaneous measurements of the thermal electron 
temperature and the solar wind velocity in the solar corona. The inaugural mission of 
MACS was beset with problems associated with instrumental scattering, lacking an 
imaging device to image the eclipse to ascertain tracking reliability and the possible 
atmospheric dust pollution due to the experimental location in a dry terrain devoid of 
vegetation. However the remedial measures are very simple in nature and tied to adequete 
financial resources. Using MACS with a coronagraph will also eliminate the time and 
terrain constraints associated with eclipses. Using MACS in a space platform that could 
take it very close to the solar surface and look into the corona will also eliminate the need 
for a coronagraph and the problems associated with the F-corona. This would be very 
promising in light of the scattering and the F-corona concerns associated with terrestrial 
observations. The remedial measures and the uses of other platforms for MACS are 
detailed in chapter-7. 
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Chapter 5 


MODIFIED CRAM’S THEORY AND ITS DEPENDENCE ON VARIOUS 

PARAMETERS 


5.1 Modified Cram’s theory 

Lawrence Cram (1976) formulated the theory for the formation of the K-corona. 
The modified Cram’s theory includes the influence of the solar wind velocity in the 
formation of the K-corona, which is a continuous spectrum. The reasons for the 
formation of this continuous spectrum and the influence of the solar wind velocity are 
discussed in chapter-2. Figure (5.1) is a schematic diagram showing the mathematical 
description of the scattering phenomenon where the photospheric radiation from the sun 
is scattered by the free electrons lying along the line of sight of an observer. This 
scattering phenomenon is popularly known as Thomson scattering while the scattered 
radiation contributes to the formation of the K-corona. A detailed treatment of the 
scattered intensity, which includes the influence of the solar wind velocity, is provided in 
Appendix-A. Equation (5.1) gives the intensity of the radiation scattered at an observed 
wavelength (X ) by the free electron density distribution along the line of sight at height 
( p ) above the solar limb. 
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Figure (5.1). This is a schematic diagram showing the 
mathematical description of the scattering phenomenon where 
the photospheric radiation from the sun is scattered by the free 
electrons lying along the line of sight of an observer. The line of 
sight is at a distance ( p ) from the center of the sun. 
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( 5 . 2 ) 
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The distances r, p and x are expressed in solar radii as follows. 


r — » rR 


solar 


P P R solar 


x — » xR 


solar 


(5.3) 


The physical parameters that are inherent in equation (5.1) are given in equation (5.4). 


Uj(X r ) = limb darkening coefficient 

f(X') = extraterrestrial solar irradiance 
N e (rR solar ) = electron density model 
T(rR SO | ar ) = coronal temperature model 
W (rR SO | ar ) = solar wind model 


The detailed derivation of equation (5.1) (=Appendix-A.98) is shown in 
Appendix-A. The computer code written in IDL to solve equation (5.1) is given in 
Appendix-B. From the physical parameters that are inherent in equation (5.1), as shown 
in equation (5.4), it is evident that the shape of the theoretical K-coronal intensity 
spectrum is dependent upon the models used for the coronal temperature, the solar wind 
velocity and the electron density distribution structures, hi addition the numerical 
outcome of equation (5.1) can also be affected by the numerical methods employed in 
solving the various integrals. However the extraterrestrial solar irradiance spectrum and 
the limb darkening coefficients could be considered as well-established quantities. As 
such, it is prudent to analyze the dependence of the shape of the theoretical K-coronal 
intensity spectrum on the coronal temperature, the solar wind velocity and the electron 
density distribution structures and the numerical methods used in solving the integrals. 
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5.2 Dependence on the electron density distribution function 

The bright corona seen during a solar eclipse is due to the scattering of the 
photospheric light by the free coronal electrons (F-corona), by dust (F-corona) and the 
emission lines (E-corona). The F-corona can be easily eliminated with the use of 
polarizes while the E-corona could be identified with the prominent peaks in an otherwise 
continuous spectrum. These peaks could be easily removed, as demonstrated in figure 
(4.33), to isolate the K-corona. The point that needs to be made here is that the most 
straightforward way of measuring the coronal electron density distribution is by 
measuring the K-coronal intensity distribution during a solar eclipse that is filtered 
through polarizes orientated at three different angles to eliminate the K-corona. This 
information can then be used to determine the electron density distribution that would 
have yielded the measured intensity distribution due the scattering of the photospheric 
spectrum by the free coronal electrons. This method is most suitable for electron density 
distribution during the maximum phase of the sun, which allows for the corona to assume 
azimuth symmetry. 

One of the pioneering efforts in this regard was due to Baumbach (1937). 
Baumbach, from an analysis of the photometric material of ten eclipses from 1905 to 
1929, deduced the expression given in equation (5.5) for the distribution of brightness in 
the solar corona. This calculation was based upon averaging for all the observations and 
for all values of the position angle with the brightness of the center of the solar disk taken 
as 10 6 . 
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I(r) = 


( 5 . 5 ) 


0.0532 1.425 2.565 

„ 2.5 + 7.0 + 17.0 

r r r 

From equation ( 5 . 5 ) it is apparent from the last term that the brightness falls off 
most rapidly in the innermost region of the corona while at large distances the brightness 
is proportional to the first term, which decreases fairly slowly. This formula gives only a 
smoothed, average value for the distribution of the brightness. On a brighter note 
November and Koutchmy (1996) report of a best ever achieved white-light coronal 
observations. This unique opportunity was afforded by the path of total solar eclipse of 
11 July 1991 passing above the 3.6 m Canada-France-Hawaii Telescope (CFHT) on 
Mauna Kea, Hawaii, which allowed for high spatial resolution white-light coronal 
observations with the most modern detectors that was denied to Baumbach (1937). The 
sun was at its maximum phase during this eclipse. Like Baumbach (1937), November 
and Koutchmy (1996) use the sum of three power -law terms to best approximate three 
regions of the curve, namely, r < 1.1 SR, 1.1 SR < r < 1.5 SR, and r > 1.5 SR (SR = 
solar radii), to obtain the brightness distribution given by equation ( 5 . 6 ). Here again I(r ) 
is in units of 10 6 of the solar disk center and (r ) in units of solar radii. 


I(r) = 


0.0551 1.939 3.670 

„ 2.5 + 7.8 + 18.0 

r r r 


( 5 . 6 ) 
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Figure (5.2) shows the comparison between the brightness distributions given by 
equation (5.5) and equation (5.6). November and Koutchmy (1996) conclude that the 
small difference in the slope to be smaller than the uncertainties in the measurements. 


Log (Intensity) VS Solar radii 

Log (I) 



Figure (5.2). A plot showing the solar radial log-intensity 
variation as a function of solar radii (r) from the center of the 
sun. The solid curve represents the measurements by 
November and Koutchmy (1996) while the dashed curve 
represents the measurements by Baumbach (1937) based on 
several eclipses observations. 


The coronal medium is optically thin. For Thomson scattering in optically thin 
plasma, the measured K-coronal intensity is proportional to the integral through the line 
of sight of the electron density times a local dilution factor, which is the solar intensity 
illuminating the scattering region of the corona. One consequence of this property is that 
the average radial intensity variation determines the mean radial variation of the electron 
density. Figure (5.3) shows the electron density distribution derived by Baumbach 
(1937) and November and Koutchmy (1996). Here an isotropic scattering function is 
assumed. This is not a bad choice for the inner corona since the illum in ation dilution 
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factor becomes equivalent to the nearly isotropic scattering form because it represents a 


superposition of intensity contributions from a wide range of solar disk angles. 


Log(N e ) cm 


Log (Electron number density) VS Solar radii 



U 15 1.5 

(Solar radii) 


Figure (5.3). A plot showing the log of the electron density variation 
as a function of the solar radii (r) from the center of the sun. The 
solid curve represents the calculations by November and Koutchmy 
(1996) while the dashed curve represents the calculations by 
Baumbach (1937) based on several eclipses observations. 


The electron number density derived by Baumbach ( 1937 ) is given by equation (5.7). A 
detailed derivation of equation (5.7) is also shown in Shklovskii ( 1965 ). 


N (r) = 10 8 x ( 


0.036 1.55 2.99 


.1.5 


- + - 


.6.0 


- + - 


.16.0 


) cm 


-3 


(5.7) 


It is true that the total brightness of the corona changes from one eclipse to another. 
Generally speaking the coefficients of equation (5.7) too should change from eclipse to 
eclipse. At sunspot maximum, that is when the sun is at its maximum phase, the corona 
has approximately a circular form and uniformly bright as shown in figure (5.4). At 
sunspot minimum, that is when the sun is at its minimum phase, the corona is elongated 
in the equatorial regions and brighter than the polar regions as shown in figure (5.5). 
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Figure (5.4). A photograph of the solar 
corona during the total solar eclipse of 
11 August 1999 where the sun was 
approaching the maximum phase. 
Here the corona is uniformly bright 
and circular in shape. 



Figure (5.5). A photograph of the solar 
corona during the total solar eclipse of 
4 November 1994 where the sun was 
approaching the minimum phase. Here 
the corona is elongated in the 
equatorial regions and brighter than 
the polar regions. 
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For several reasons as listed below, the electron number density distribution given in 
equation (5.7) was selected as input parameter in equation (5.1) for the computation of 
the theoretical K-coronal intensity. 

(a) This derivation of the electron number density given in equation (5.7) was based 
upon ten eclipse observations. Although the observations were performed almost 60 - 
80 years ago the shape of the density distribution was reconfirmed by November and 
Koutchmy ( 1996 ) using the most modem telescope for the same purpose. Since the 
reconfirmation was in conjunction with the sun in maximum phase in July 1991 , this 
density distribution could be considered suitable for the total solar eclipse of August 
1999 , which again coincided with the maximum phase of the sun. The comparison 
between the electron number densities determined by Bauntbach ( 1937 ) and 
November and Koutchmy ( 1996 ) is shown in figure (5.3). 

(b) During the maximum phase of the sun the coronal brightness is approximately 
uniform all around the sun, which is not true for the sun at the minimum phase where 
the coronal brightness is concentrated in the equatorial regions. For a sun at the 
minimum phase a single density distribution function may not hold for all latitudes. 
Most of the other density models are specific regions and in particular hold for the 
minimum phase of the sun, which are detailed below. Therefore in conjunction with 
the maximum phase of the sun where the coronal brightness is uniform all around the 
sun the density model given by equation (5.7) seem to be most appropriate. 
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(c) The method employed in determining the electron density given by equation (5.7) 
is an average over so many other coronal features. In determining the theoretical K- 
coronal intensity from equation (5.1) it is difficult to integrate over the different 
coronal features that may cross the line of sight. For such a calculation to be possible 
individual coronal features have to be identified along with their electron density 
contrast. In light of the difficulties in such a process it is again most appropriate to 
use an electron density model averaged over all such coronal features that may occur 
in conjunction with the maximum phase of the sun. 

(d) The very definition of the white-light corona is the brightness produced by the 
scattering of the photospheric light by the free coronal electrons, which is identified 
as the K-corona. Although this may be contaminated by the scattering by dust 
particles, which is identified as F-corona, this contamination can easily quantified by 
the use of polarizes and eliminated. Also in the region of interest of the experiment 
described in this dissertation, which is from the limb to 1.5 solar radii, the 
contribution by F-corona is negligible. Therefore it is prudent to use an electron 
number density calculation determined by the measurement of the brightness of the 
white-light corona, which is the case for the electron number density determined by 
Baumbach (1937) and November and Koutchmy (1996). 
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The electron density model used Cram (1976), in formulating the properties of the K- 
coronal intensity spectrum for the determination of the thermal electron temperature in 
the corona, is given in equation (5.8). 

N e (r) = 1.67xlO (4 - 0+4 - #4/r) cm 3 (5 - 8) 

Ichimoto et al. (1996) used the same electron density model given in equation (5.8) in 
determining the coronal temperatures. Their eclipse observation of February 1994 
coincided with the sun at its minimum phase where the coronal brightness is prominently 
displayed in an elongated shape in the equatorial regions of the sun. 

According to Cram (1976) the model given in equation (5.8) agrees with the Van de 
Hulst (1950) minimum equator model within 2% for line of sight between 1.5 and 2.0 
solar radii. This model is appropriate for the temperature determination by Ichimoto et al. 
(1996) because their observations were confined to regions between 1.5 and 2.0 solar 
radii and also their observations on the total solar eclipse of February 1994 coincided 
with the minimum phase of the sun. 

Another model for the electron number density was due to Newkirk (1961). This 
model was based on radio and optical observations. Here the effects of refraction were 
considered while neglecting the effects of magnetic fields. The optical depth effects 
inherent in any radio observations are discussed in section (1.2). Separate values for the 
electron number densities were obtained for the average corona, the polar caps and the 
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active regions. The approximate formula for the electron number density presented by 


Newkirk (1961) is shown in equation (5.9). 


N e (r) = 4.2xlO <4 ‘ 0+4 ‘ 32/r> cm' 3 


(5.9) 



Figure (5.6). A series of plots showing the radial dependence 
of the mean electron density in coronal holes computed 
from various sets of polarization brightness data. This plot 
was obtained from Cranmer et al. (1999). 
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Various other research groups have computed radial dependence of the mean electron 
density specific to coronal holes. These calculations were based on various sets K- 
coronal polarization brightness. Figure (5.6) is a plot from Cranmer et al. (1999) showing 
the electron number density profiles calculated by various research groups. 

In figure (5.6) the heavy solid line is due to Cranmer et al. (1999) using the 
UVCS/WLC aboard the SOHO satellite. And the electron number density profile 
obtained by Cranmer et al. (1999) is given by equation (5.10). The observed data in this 
case is restricted to regions greater than 1.5 solar radii where (r) is in solar radii. 

N,(r) = 10 s x (3.890 xr' 1#i + 0.00869 xr' 157 ) cm' 3 (5.10) 

The filled circles in figure (5.6) is due to Strachen et al. (1993) using rocket based 
coronograph data. Here they considered the corona to be unstructured and spherically 
symmetric. However this assumption is reasonable for considerations restricted to the 
coronal hole. Here again electron density profile is restricted to the coronal hole region at 
heights above 1.5 solar radii. The electron number density profile obtained by Strachen et 
al. (1993) is given by equation (5.11) where (r) is in solar radii. 

N e (r) = 10 8 x(0.152xr“ 6 ' 71 + 0.0193 xr' 2 ' 40 ) cm 3 (5.11) 
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Using data obtained by Spartan 201-01, Guhathakurta and Holzer (1994) have 
presented an electron density profile for the polar coronal hole region specific to the 
minimum phase of the solar cycle. This again is valid for heights above 2.0 solar radii 
from the center of the sun. The electron density model by Guhathakurta and Holzer 
(1994) is given by equation (5.12) where (r) is in solar radii. 

N e (r) = 10 8 x (0.81 xr“ 16 87 +0.80 xr‘ 8 45 + 0.014xr 28 ) cm 3 ^ 5 - 12 ) 

The electron density measurements by Fisher and Guhathakurta (1995) were based on 
the use of both the Mk-III K-coronameter from Mauna Loa and White Light Coronagraph 
of the space-bome Spartan 201. This study coincided with the maximum phase of the sun 
in 1990. However the observations were restricted to the northern and southern coronal 
holes. Their measurements extended from 1.16 solar radii to 5.5 solar radii. The 
uncertainties reported for the electron density measurements at 1.16 - 1.3 and 1.4 - 1.8 
solar radii are -15% and (30 % - 35 %), respectively. This latter region coincided with 
data set where Mk-III and Spartan 201 data were cross calibrated. Figure (5.7) is a plot of 
the electron density measurements by Fisher and Guhathakurta (1995). 
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Radial distance r in units of Ro 


Figure (5.7). A plot showing the electron 
densities for the north (solid line) and south 
polar coronal holes (dashed line) and the north 
(dotted line) and south (dot-dashed line) polar 
coronal rays as a function of height. Reproduced 
from Fisher and Guhathakurta (1995). 


Figure (5.8) is a comparison of the electron number density functions based on the 
models by Baumbach (1937, equation (5.7)), Ichimoto et al. (1996, equation (5.8)), 
Newkirk (1961, equation (5.9)), Cranmer et al. (1999, equation (5.10)) and Guhathakurta 
and Flolzer (1994, equation (5.12)). Figure (5.9) is a comparison of the shapes of the K- 
coronal intensity spectra based on the Baumbach (1937, equation (5.7)), Ichimoto et al. 
(1996, equation (5.8)), Newkirk (1961, equation (5.9)), Cranmer et al. (1999, equation 
(5.10)) and Guhathakurta and Holzer (1994, equation (5.12)) models for the electron 
number density. These plots are normalized to intensity at 4000.0 angstrom and modeled 
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for line of sight at 2.0 solar radii above the center of the sun. This plot is simply for the 
comparison of the shapes of the theoretical K-coronal intensity spectra based on different 
electron density models. The theory by Cram (1976) for the determination of the thermal 
electron temperature of the solar corona is based upon the shape of the K-coronal 
intensity spectrum and not its absolute values. It is evident from figure (5.9) that no 
discernible differences could be seen in the shapes of the theoretical K-coronal spectra 
for five different electron density models. 


Electron number density VS radius 
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Figure (5.8). A comparison of the electron number density 
functions based on the models by Baumbach (1937, equation 
(5.7)), Ichimoto et al. (1996, equation (5.8)), Newkirk (1961, 
equation (5.9)), Cranmer et al. (1999, equation (5.10)) and 
Guhathakurta and Holzer (1994, equation (5.12)). 
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Modeled observed intensity for different electron number density profiles VS wavelength 
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Figure (5.9). Plots of the theoretical K-coronal intensity spectra 
based on the electron number density functions given by 
Baumbach (1937, equation (5.7)), Ichimoto et al. (1996, equation 
(5.8)), Newkirk (1961, equation (5.9)), Cranmer et al. (1999, 
equation (5.10)) and Guhathakurta and Holzer (1994, equation 
(5.12)) and the line of sight at 2.0 solar radii from the center of 
the sun. The curves have been normalized at 4000.0 angstrom. 


Figure (5.10) is another plot the theoretical K-coronal intensity based on different 
density models for the line of sight at 1.1 solar radii. Here again no discerning differences 
could be observed in the shapes of the K-coronal intensity distribution based on the 
different electron density models. However the absolute values would differ for the 
different electron density models although it is irrelevant to this experiment. 
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Modeled observed intensity for different temperature profiles VS wavelength 



Figure (5.10). Plots of the theoretical K-coronal intensity spectra 
based on the electron number density functions given by 
Baumbach (1937, equation (5.7)), Ichimoto et al. (1996, equation 
(5.8)), Newkirk (1961, equation (5.9)), Cranmer et al. (1999, 
equation (5.10)) and Guhathakurta and Holzer (1994, equation 
(5.12)) and the line of sight at 1.1 solar radii from the center of 
the sun. The curves have been normalized at 4000.0 angstrom. 


In summary the electron density measurements given by equations (5.8) to (5.12) are 
specific to minimum phase of the sun, or specific to the coronal holes or specific to 
heights above the region of interest of MACS. The region of interest of MACS extends 
from 1.1 solar radii to 1.5 solar radii and the experiment coincided with the maximum 
phase of the sun. In this regard it was decided to select the electron number density given 
by Baumbach (1937), as given by equation (5.7), as input parameter for the theoretical 
calculation of the K-coronal intensity spectrum. Here the calculations were based on ten 
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eclipse observations, which avoids the scattering inherent with coronagraphs, and the 
region of interest extended from 1.0 solar radius to ~3.0 solar radii. Although these 
measurements were conducted about 60 - 80 years ago, this distribution function was 
reconfirmed by November and Koutchmy (1996). They had a unique opportunity of 
performing the same experiment using the most modem 3.6-m CFHT telescope in Hawaii 
in conjunction with the total solar eclipse of July 1991, where the sun was at its 
maximum phase. 


5.3 Dependence on the temperature profile 

Here what need to be investigated are the differences in the theoretical K-coronal 
intensity profiles based on the assumption of an isothermal corona and temperature 
profiles with a radial dependence. Usually the radial dependence of the coronal 
temperature profiles is based on other experimental results. Since the purpose of this 
dissertation is to determine the temperature and the wind profiles at several coronal 
heights and latitudes, it will not be helpful in substituting radial temperature profiles 
determined experimentally. However an important theoretical consideration for the 
temperature profile of the solar corona is due to Chapman (1957). He assumed that above 
a certain height in the corona losses by emission to be insignificant compared to 
conductive losses and derived the following radial dependence for the temperature profile 
as given by equation (5.13) where (r) is in solar radii. 


T(r) = T # xr” 2/7 


(5.13) 
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For comparison between an isothermal corona and a temperature profile given by 
equation (5.13) assume the line of sight at 1.3 solar radii with the electron number density 
profile given by equation (5.7). Also assume that for the case of the temperature profile 
given by equation (5.13) that the temperature at the intersection between the line of sight 
at 1.3 solar radii and the plane of the solar limb to be 1.0 MK and then to fall off as r 2/7 . 
For such a case the temperature profile is given by equation (5.14) and its variation along 
the line of sight is given by figure (5.11). Here (r) is the distance to points along the line 
of sight from the center of the sun. The isothermal corona assumes a temperature of 1.0 
MK all along the line of sight. 


Temperature VS radius 



Figure (5.11). This is a plot of the temperature profile along the 
line of sight at 1.3 solar radii based on the assumption of a 
conductive corona Chapman (1957). The temperature at the 
point of intersection of the line of sight and the plane of the solar 
limb is 1.0 MK. 
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Figure (5.12) is a plot of the theoretical K-coronal intensity spectra due to an isothermal 
corona of temperature 1.0 MK and due to temperature profiles given by 
T(r) = 1.0x(1.3/r) 2/7 and T(r) = 1.0+ (1.0 + 9.0* (r -1.75) 2 )' 1 . The latter is a 
hypothetical temperature profile assumed by Cram (1976). From figure (5.12) no 
discernible differences could be seen between the spectra due to an isothermal corona and 
the model with a temperature profile given by T(r) = l.Ox (1.3/r) 2/7 . 
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Figure (5.12). This is a plot of the modeled K-coronal intensity 
spectra for an isothermal corona of 1.0 MK and a model with 
the temperature profiles given by T(r) = l.Ox (1.3/r) 2 7 and 
T(r) = 1.0 + (1.0+9.0*(r-1.75) 2 ) -1 . Here (r) is the distance 
between the center and points along the line of sight 1.3 solar 
radii. The temperature at the point of intersection of the line of 
sight and the plane of the solar limb is 1.0 MK. 
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5.4 Dependence on the wind profile 

Here again what need to be investigated are the differences in the theoretical K- 
coronal intensity profiles based on the assumption of an isothermal corona and wind 
profiles with a radial dependence. Usually the radial dependence of the coronal wind 
profiles is based on other experimental results. However one important wind profile 
based purely on theoretical consideration is due to Parker (1958). The wind structures 
due to isothermal and non-isothermal coronal conditions, driven by sound waves, driven 
by dust, driven by lines, driven by the combined effect of the magnetic fields and rotation 
and driven by Alfven waves are detailed in Laniers & Cassinelli (1999). 

Consider an isothermal wind in which the gas is subject to two forces: the inward 
directed gravity and the outward directed gradient of the gas pressure. It has been shown 
that the momentum equation has many solutions, depending on the boundary conditions 
but only one of them, the critical solution starts subsonic at the lower boundary of the 
wind and reaches supersonic velocities at large distances. The solution to the wind 
structure in such a case is given by equation (5.14). 


v x exp( — — ) = a 




a 


v r j 


xexp 


r 2 


where 


r_ = 


GM 


solar 


2a 2 


a = 




v ^ j 


(5.14) 
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This phenomenon reveals that a wind is generated as a natural consequence of an 
isothermal corona, hi equation (5.14), G, T, p , R, a and r c are the universal 

gravitational constant, isothermal temperature, mean atomic weight of the particles 
expressed in units of the mass of the proton, universal gas constant, isothermal sound 
speed and the critical radius, respectively. For solar composition p is 0.602. For an 
isothermal temperature of 1 .0 MK the isothermal sound speed (a) and the critical radius 
(r c ) are 117.5 km/sec and 6.9 solar radii, respectively. Figure (5.13) shows the profile of 
the wind for the above case. 


Solar wind velocity (based on Porker model) VS solar radius 



Figure (5.13). Coronal wind profile based on the Parker 
model for an isothermal corona. In the ideal solution 
the velocity is equal to the isothermal sound speed at 
the critical radius and starts with subsonic velocities at 
the lower boundary and reaches supersonic velocities at 
large distances. 
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The Parker wind model shown in figure (5.13) is a natural consequence of an 
isothermal corona. Therefore, even if the solar wind is neglected in the determination for 
the theoretical K-coronal intensity spectrum for an isothermal corona, it is worthwhile to 
investigate the effect on this spectrum by the bulk flow velocity introduced by the Parker 
model naturally. In figure (5.14) the comparison is made between a theoretical K-coronal 
spectrum for an isothermal corona of 1.0 MK, zero wind velocity and line of sight at 1.5 
solar radii with the spectrum due to the bulk flow velocity that naturally arises with the 
introduction of the Parker wind model given in equation (5.14). 



w 


Modeled observed intensity for different wind profiles VS wavelength 
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Figure (5.14). Comparison of the theoretical K-coronal 
intensity spectrum for an isothermal corona at 1.0 MK, 
zero wind velocity and line of sight at 1.5 solar radii with 
the intensity spectrum due to the bulk flow velocity that 
naturally arises with the Parker wind model for an 
isothermal corona as given in equation (5.14). 
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Ill figure (5.14) no discernible differences could be observed between the two K- 


coronal intensity spectra. The reasons for this are twofold. First, from the Parker wind 
profile shown in figure (5.13) it is evident that the solar wind velocity at 1.5 solar radii is 
very small. Second, although the wind profile shown in figure (5.13) indicates that the 
wind values are significant at large distances, nevertheless, this effect is negligible due to 
the drop in the electron number density at large distances. Figure (5.15) shows the 
theoretical K-coronal intensity spectra due to isothermal coronal temperatures of 1 .0 and 
1.5 MK with the inclusion of the Parker wind model for an isothermal corona given in 
equation (5.14). As expected the nodes begin to appear mainly due to the temperature 
differences. 


Modeled observed intensity for different wind profiles VS wavelength 



Figure (5.15). Comparison between two theoretical K- 
coronal intensity spectra due to isothermal coronal 
temperatures of 1.0 and 1.5 MK and with the inclusion of 
the Parker wind model for an isothermal corona. 
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5.5 Dependence on the numerical method 

The theoretical model for the K-corona, as shown in equation (5.1), consists of 
four coupled integrations. In the absence of an analytical solution this equation could be 
solved only through a four-dimensional numerical quadrature. While acknowledging the 
existence of multitude of numerical methods, the numerical methods used in solving 
equation (5.1) were based on procedures presented in the Handbook of Mathematical 
Functions (Ed. M. Abromowitz and L.A. Stegun). Two of the integrals were calculated 
using the Trapezoidal Composite Rule while the other two were solved using a Gaussian 
Quadrature. There is reason to believe in the results since they can reproduce the 
theoretical results obtained by Cram (1976) and Ichimoto et al. (1996). This numerical 
procedure was also sufficient to determine the polarization components to confirm that 
the degree of polarization is almost wavelength independent and that it increases with 
heights above the solar limb. The degree of polarization (p) is defined by equation (5.15). 


p(^p) = 




where 


^ a (X, p) = Intensity parallel to the radial plane 

(X,p) = Intensity perpendicular to the radial plane | 


(5.15) 


Figure (5.16) shows a plot of the degree of polarization, for a given isothermal coronal 
temperature of 1.0 MK, for lines of sight at 1.1, 1.3 and 1.5 solar radii. It is evident from 
figure (5.16) that the polarization components were calculated with sufficient accuracy to 
confirm the behavior expected of the polarization. That is, they are wavelength 
independent and increase with heights above the solar limb. 
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Figure (5.16). Plot of the polarization from theoretical 
modeling for the formation of the K-corona, at a given 
isothermal coronal temperature of 1.0 MK, for lines of 
sight at 1.1, 1.3 and 1.5 solar radii. The behavior of the 
polarization conforms to the observations. That is, the 
polarization is almost wavelength independent and 
increases with coronal heights. 
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The numerical code to solve equation (5.1) was written in IDL (Interactive Data 
Language) and took -100 minutes to calculate the intensity spectrum from 3700 - 4700 
angstroms in a Pentium-II 233 MHz PC. Of the four integrals two were solved using the 
Trapezoidal Composite Rule while the other two were solved using Gaussian Quadrature. 
In the situation where the Trapezoidal Composite Rule was used, the quadrature points 
were divided into higher and higher number of points to check for differences in the 
results. The final number of quadrature points was selected based on the point where the 
overall results for the intensity converges to within five decimal places. As for the other 
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two integrals, these were solved using Gaussian quadratures. Their integration limits 
allowed for the use of Hermite polynomials and its associated weights. That is, these 
integrals had the form given by equation (5.16). 


+'*’ , n 

Jf(x)e -X dx ~ where 

i=0 

Iq = zeros of Hermite polynomials 
W; = weights 


(5.16) 


However the line of sight integration was of the form given by equation (5.17), which 
was converted to the form given by equation (5.16), which follows the procedure given in 
the Handbook of Mathematical Functions (Ed. M. Abromowitz and L.A. Stegun, pg. 
924). 


+oo +oo n 

Jf(x)dx = J{f(x)e + * 2 }e~ x dx * £{f(h, > h ' }w, where 

i=0 

hj = zeros of Hermite polynomials 
Wj = weights 


(5.17) 


In equation (5.17) the function f(x) is of the following form given by equation (5.18), 
which is a convergent function without any singularities in the region of interest. 



(5.18) 
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For example assume p =1 .0 and p=l.l in equation (5.18). Then the equation 


(5.18) reduces to equation (5.19). 
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(5.19) 


Solving equation (5.19) by using the trapezoidal composite rule, where the integration 
limit is restricted to x = ±5000 and dividing this region into 100000 intervals of width 
0.1, the result yields a value of 1.1954237. The restriction in the x-range is perfectly 
valid because the integrand rapidly falls to zero about x=0. Performing this operation 
over a range of values for p is very time consuming. However by using the Gaussian 
quadrature over a Hermite polynomial roots of just twenty points, as specified in equation 
(5.17), the result yields a value of 1.1900812. This value is accurate to within 1.0 % of 
the result yielded by the trapezoidal composite nde and is also computed very much 
faster. This is also acceptable in light of the limits imposed by the instrument. An 
analytical solution also exists for equation (5.19) as shown in equation (5.20) (Standard 
Mathematical Tables and Formulae). 
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(5.20) 


165 





Writing equation ( 5 . 19 ), as shown in equation ( 5 . 21 ), and using equation ( 5 . 20 ) yields a 
value of 1.1954237, which matches with the value obtained by using the trapezoidal 
composite rule. 



The detailed derivation of the theory and the code written in IDL is provided in 
Appendix - A and B, respectively. Throughout the code the accuracy of the calculations 
were checked using the built in CHECK_MATH function, which returns the accumulated 
math error status, and manually. There were no errors reported. 
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Chapter 6 


COMPARISON BETWEEN OBSERVATION AND THEORY 

6.1 An unexpected phenomenon noticed in the observations 

As shown in figure (3.4) the fiber optic plate in the focal plane of the telescope 
contained a single fiber in the center of the frame. During the observations this fiber was 
positioned at the center of the lunar shadow. The purpose of this fiber was to measure the 
background signal and was to be subtracted from the signals measured by the other fibers 
that were exposed to the coronal light. 

However the shape of the signal measured by the fiber located in the center of the 
lunar shadow only provided a clue that the observations were affected by scattering. 
Figure (6.1) shows the signal recorded by the fiber at the center of the lunar shadow. 
This prominently shows the peaks of the Calcium H and K lines. The counts recorded by 
this fiber are not negligible compared to the counts recorded by fibers exposed to the 
coronal light. This indicates that we have a problem with scattering light. The extent of 
scattering could not be envisaged in advance for the simple reason that the true 
observational condition could not be replicated in the laboratory. However Cram (1976) 
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pointed out that the bright sources of scattered light at the time of a solar eclipse are 
prominences and the chromosphere, and these have a spectrum quite unlike that of either 
the F or K corona. The only remedial measure that was taken in order to minimize 
scattering was to maintain clean optical surfaces of the optical components in MACS. 
Although the extent of scattering in the different wavelength regions is very difficult to 
quantify, nevertheless, it is a hard lesson to be learnt for any future observations. In this 
regard some of the remedial measures for any future observations are highlighted in 
section (6.3). 


Fiber in the center of the lunar shadow 



Figure (6.1). The spectrum recorded by the fiber located in 
the center of the lunar shadow during the eclipse. The x- 
coordinate corresponds to wavelength scale. The two 
prominent peaks are associated with the Calcium K & H 
lines corresponding to 3933.7 and 3968.5 angstroms, 
respectively. 
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6.2 Comparison between theory and observation 

In chapter-2 temperature and wind sensitive intensity ratios were identified at 
particular wavelengths. These intensity ratios were at wavelength positions 
4100.0/3850.0 and 4233.0/3987.0, respectively, for temperature and wind measurements. 
In chapter-2, we also showed how to determine the temperature and the wind velocity 
using the above intensity ratios calculated from the K-coronal intensity spectrum. In 
contrast to our ratio technique, Ichimoto et al. (1996) chose to determine the coronal 
temperature by fitting theoretical isothermal K-coronal models to the observed K-coronal 
spectrum, as shown in figure (6.2). 



Figure (6.2). Theoretical isothermal K-coronal models fitted 
to the observed K-coronal intensity spectrum to determine 
the temperature by Ichimoto et al. (1996). The lower curves 
show the differences between the observational spectrum and 
the theoretical models. The model that showed the least 
difference determined the temperature. 
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The same temperature measurement, as shown in figure ( 6 . 2 ), could also be 
determined from the measurement of the temperature sensitive intensity ratio. Such ratio 
calculations were performed using filters centered at the temperature sensitive 
wavelength positions by J.M. Pasachoff and his team from Williams College, 
Massachusetts during the total solar eclipse of 11 August 1999 in Romania. N.L. 
Reginald and J.M. Davila performed the theoretical calculations for predictions of the 
temperature for the intensity ratios for J.M. Pasachoff and his team. No results have yet 
been published by Pasachoff et al. 

In chapter-4 the intensity ratio technique was employed to determine the 
temperature and the wind velocity for several fiber locations. The wind measurements 
were considered a total failure in this very first attempt to simultaneously and globally 
determine both temperature and wind velocity at twenty different locations on the solar 
corona. The reason for this is that the wind measurement is very sensitive to the intensity 
at 4233.0 angstrom, which is at the high end of the wavelength measurement by MACS. 
As revealed in figure ( 6 . 1 ) the instrumental scattering was high at the high end of the 
wavelength region thus predicting very high values for the wind velocity measurements. 
As for the temperature it too suffered from the instrumental scattering. The reason for this 
is that the temperature is sensitive to the intensity measurement at 3850.0 angstrom and 
4100.0 angstrom. As revealed in figure ( 6 . 1 ) the strong Calcium Hand K lines in the 
region of 3900-4000 angstrom seem to have contributed strongly to instrumental 
scattering. Predicting the contribution due to scattering on the coronal intensity measured 
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by individual fibers is further compounded by its dependency on wavelength. At this 
juncture the knowledge of instrumental scattering only provides awareness of the 
problem of instrumental scattering and remedial measures that need to be considered 
before any future observations. The sole handicap in such observations is the inability to 
replicate the experiment in the laboratory in the absence of a coronagraph. 

Figure ( 6 . 3 ) shows the comparison between experimentally obtained K-coronal 
spectrum for fiber # 04 with theoretical models for various isothermal temperatures. 
From the intensity ratio method the wind velocity and the temperature deduced for fiber # 
04 are 300.0 km/sec and 1.29 MK, respectively. The theoretical K-coronal models plotted 
in figure ( 6 . 3 ) correspond to a wind velocity of 300.0 knt/sec and isothermal coronal 
temperatures of 0.5MK, 1.0 MK, 1.5 MK and 2.0 MK. In such a scenario the 
experimentally obtained K-coronal spectrum may be expected to lie between the 
theoretical models for 1.0 MK and 1.5 MK. However it is obvious that the experimental 
K-coronal spectrum deviates from the predictions. This deviation can be divided into four 
distinct regions on the wavelength scale. 

1. 3800 - below: In this region the experimental curve falls below the predictions. 

2 . 3800 - 4000: In this region the experimental curve rises above the predictions. 

3 . 4000 - 4200: In this region the experimental curve falls below the predictions. 

4 . 4200 - above: In this region the experimental curve rises above the predictions. 
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Figure (6.3). This plot shows the comparison between 
experimentally obtained K-coronal spectrum for fiber # 04 with 
theoretical models for various isothermal temperatures. From the 
intensity ratio method the wind velocity and the temperature 
deduced for fiber # 04 are, 300.0 km/sec and 1.29 MK, 
respectively. The theoretical K-coronal models plotted in figure 
(6.3) correspond to a wind velocity of 300.0 km/sec and 
isothermal coronal temperatures of 0.5MK, 1.0 MK, 1.5 MK and 
2.0 MK. In such a scenario the experimentally obtained K- 
coronal spectrum may be expected to take a position between the 
theoretical models for 1.0 MK and 1.5 MK. (Expt: = 
Experimental, Theo: Theoretical) 


The comparison between observation and theory for the K-coronal spectrum recorded 
by fiber # 04, as shown in figure (6.3), needs to be considered in conjunction with the 
extent of instrumental scattering that is depicted in figure (6.1). Figure (6.1) too can be 
divided into four distinct regions. Assuming that the signal expected of this fiber located 
in the center of the lunar shadow to produce a straight line with a positive gradient 
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signaling a higher background count at higher wavelengths. Based on this assertion, and 
assuming the signal recorded in the 4000 - 4200 angstrom to represent the correct true 
background signal, figure (6.2) can be divided into four distinct regions. 

1. 3800 - below: Some optical component in MACS seems to possess a very low 
sensitivity to the blue region of the spectrum. 

2. 3800 - 4000: The Calcium II K & H lines contribute a large signal in this region. 

3. 4000 - 4200: Assumed to represent the true background count. 

4 . 4200 - above: Scattering is again significant in the high end of the spectrum. 

The above properties pertaining to the four distinct regions in the scattered spectrum 
is similar to the pattern of deviation of the observed K-coronal spectrum of fiber # 04 
from the theoretical models, as shown in figure ( 6 . 3 ). That is, the very low blue 
sensitivity in the wavelength region 3800 - below and the higher counts from the 
scattering of the Calcium II K & H lines in the region 3800 - 4000 angstrom may have 
caused the observed K-coronal spectrum to dip and rise, respectively, in these regions. 
The scattering associated with the high end of the wavelength region, i.e. 4200 - above, 
may have caused the observed K-coronal spectrum to again rise above the predicted 
level. Figure ( 6 . 4 ) is a plot of the superposition of figure ( 6 . 1 ) and figure ( 6 . 3 ). This 
shows that there is a certain amount of similarity between, on the one hand, the shape of 
the scattered spectrum, and on the other hand, the deviation between the predicted and 
observed K-coronal spectra. In figure ( 6 . 4 ), all spectra are normalized to unity at 4000.0 
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angstrom. It is reasonable to infer from figure (6.4) that some proportion of the scattered 
light spectrum recorded by the fiber located in the center of the lunar shadow has 
contaminated spectrum recorded by fiber #04. This could explain the deviations between 
the observed spectrum and the predictions. 


Comparison of shapes between the scattered light spectrum and the spectrum of fiber #04- 



3900 4000 4100 4200 

angstroms 

Figure (6.4). A plot of the superposition of figure (6.1) and figure 
(6.3). The shape of the spectrum recorded by the fiber located in 
the center of the lunar shadow corresponds with the shape of the 
observed K-coronal spectrum by fiber #04. The expected location 
of the observed spectrum is between the theoretical isothermal K- 
coronal spectra for 1.0 MK and 1.5 MK. 


Figure (6.5) shows the superposition of the observed K-coronal spectrum for fiber 
#04 and the theoretical K-coronal spectrum for a wind velocity and isothermal 
temperature of 300.0 km/sec and 1.29 MK, respectively, and the difference between the 
two. The wind velocity and isothermal temperature of 300.0 km/sec and 1.29 MK, 
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respectively, for fiber #04 were determined from the intensity ratio method. From figure 
(6.5) it is apparent that the deviations between the observational and the theoretical 
spectra are larger at wavelengths of 3900 - 3950 angstroms (where the scattered light 
spectrum peaks), and are smaller at wavelengths around 4100 angstrom In this regard, the 
deviations between observed and predicted spectra exhibit certain gross similarities to the 
shape of the scattered light spectrum. Even a small percentage of the scattered light 
spectrum being contaminated with fiber #04 could have caused its spectrum to deviate 
from the theoretical prediction. 


Comparison between theory and observation for fiber #04 



Figure (6.5). A plot showing the observed K-coronal spectrum 
for fiber #04 along with the theoretical K-coronal spectrum for 
the wind velocity and temperature derived from fiber #04. The 
plot also includes the difference between the observational and 
the theoretical plots along with the scattered light spectrum 
observed by the fiber in the center of the lunar shadow. 
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Iii figure ( 6 . 5 ) the areas under the observational and the theoretical spectra are 
0.22168834 and 0.22163162, respectively. That is, the difference between the areas is 
almost negligible. Nevertheless, since the brightest sources of scattered light 
(chromospheric or prominence material) have spectra which are quite different from 
those of F- or K-corona, we can now see that even a slight contamination of our coronal 
data by such scattered light can have a serious effect on the quantitative interpretation of 
our results. This was pointed out by Cram (1976) who suggested that any observations 
should avoid such bright sources. However in MACS twenty fibers were exposed to all 
around the sun and at different heights above the solar limb. As a result, it is almost 
inevitable that some of our fibers were unfortunately exposed to bright prominence 
material. The effects due to exposure to these bright sources could be avoided only 
through a scatter free instrument. In this regard every effort was made to keep 
instrumental scattering to a minimum. The present experiment, the very first of its kind, 
reveals the extent and effects of instrumental scattering and the need for more remedial 
measures in this regard. The same is true of other fibers that were exposed to the coronal 
light: each was contaminated by scattered light to a greater or lesser extent. 

Table ( 6 . 1 ) gives the wind and the temperature sensitive intensity ratios for fiber #04 
from figure ( 6 . 5 ). From table ( 6 . 1 ) it is apparent that the wind and temperature sensitive 
intensity ratios are well matched in spite of the fact that the observational spectrum 
deviates from the prediction. This may be somewhat fortuitous, since if the amount of 
scattered light had been different, the observed spectmm might have matched a totally 
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different theoretical spectrum. This matter is further compounded by the non-existence of 
any mechanisms to quantify the degree of contamination of instrumental stray light on 
the individual fiber, hi this regard it appears that in spite of the problems caused by the 
instrumental scattering the temperatures ratios were somehow maintained. Thus the 
accuracy of the temperature is to within the error estimates made in chapter- (4 .7). That 
is to within ~± 6 % (± 1.2 MK). As for the wind measurements, as discussed in chapter- 
(4.7), its measurement was deemed to have been a failure in the current experimental 
effort. However the high wind values are in line with the uncertainty of ± 200 km/sec 
for a measurement uncertainty of ± 1 % . 

Table (6.1). Comparison of the wind and the temperature sensitive 
intensity ratios between the observational and the theoretical K-coronal 
spectrum, for fiber #04. These values are based on figure (6.5). 


Fiber #04 

Temp: I(4100)/I(3850) 

Wind: I(4233)/I(3987) 

Observational 

1.531 

1.242 

Theoretical 

1.529 

1.242 


Figure (6.6) shows plots of the observed K-coronal spectrum of fiber #04, as 
shown in figure (6.5), corrected for the scattered light in a cmde manner. Here various 
fractions of the scattered light recorded by the fiber located in the center of the lunar 
shadow were subtracted from the spectrum recorded by fiber #04. It is evident from 
figure (6.6) that the observed K-coronal spectrum then begins to follow closely the shape 
of the theoretical K-coronal spectrum implying the influence of the scattered light. 
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Comparison between theory and observation for fiber #04- 


Figure (6.6). Plots showing the effect on the shape of the observed K- 
coronal spectrum of fiber # 04, as shown in figure (6.5), in comparison 
with the shape of the theoretical K-coronal spectrum by subtracting 
various fractions of the scattered spectrum from the observed K-coronal 
spectrum. A, B, C, D, E and F correspond to a subtraction of 0.0%, 
15.0%, 10.0%, 8.0 %, 5.0% and 3.0%, respectively, of the scattered light 
spectrum recorded by the fiber that was located in the center of the lunar 
shadow from the spectrum recorded by fiber #04. These plots show that 
the observed spectrum closely follows the shape of the theoretical 
spectrum with various fractional subtraction of the scattered spectrum. 
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A similar pattern as above follows for all the other fibers for which temperatures 
were determined in chapter- 4. Instrumental scattering can be minimized only by paying 
careful attention to the optical quality of the components of the instrument. The stricter 
the demands on optical quality, the higher the cost will be. The very first operation of 
MACS in conjunction with the total solar eclipse of 1 1 August 1999 has shed much light 
and understanding on the effects of instrumental scattering associated with such 
observations. Apart from better optics, other remedial measures envisaged for future 
operation of MACS are listed in section (6.3). 

6.3 Remedial measures 

The principle and the method described in this thesis is unique in its character that 
allows for the simultaneous and global measurements of the thermal electron 
temperatures and the solar wind velocities on the solar corona using the shape of the K- 
coronal spectrum. However the steps described in section (4.6) in isolating the K-coronal 
spectrum from a terrestrial observation in conjunction with an eclipse are beset with 
problems. The main problems that arise from terrestrial observation during an eclipse are 
time constraints and state of the sky. Apart from these, there are other difficulties that 
arise in connection with the data analysis. We summarize these here. 

1. In order to determine the wavelength sensitivity of our instrument, we have assumed 
that the sky spectrum recorded by MACS (corrected for Rayleigh scattering) should 
match the photospheric spectrum depicted in figure (2.1). The comparison between 
the two spectra shown in figure (4.21) demonstrates clearly the nature of the 
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difficulties we face in this regard, especially close to the Ca line at 3933.7 angstrom 
and regions at 4000.0 angstrom. Since this method depends on the shape of the K- 
coronal spectrum these inadequacies, although seem to be irrelevant, do affect the 
shape of the K-coronal spectrum in a significant manner. However a near perfect 
wavelength sensitivity curve could be obtained with the use of absolute wavelength 
calibration using sources of known strengths as commonly done for all advanced 
instruments with optical elements. This suggests that in this particular experiment the 
derivation of the wavelength sensitivity curve for MACS was limited in its accuracy 
by choice. 

2. Another problem arises from our assumption that the F-corona must be allowed for in 
our Rayleigh-corrected MACS sky spectrum. The intensity of the F-corona (which 
depends on the properties of interplanetary dust between the sun and the earth) has 
some wavelength dependence, perhaps proportional to yfk . However on the 
assumption that the F-corona is unpolarized while the K-corona is polarized (although 
only orthogonal scattering gives complete polarization as per equations (A.48) and 
(A.49) there exists an easy practical means to determine the contribution from F- 
corona. This is done by observing the corona through a polarized filter with the 
polarization axis at minimum of three known angular positions. A detailed 
explanation could be found in Golub and Pasachoff (1997). The only concern that 
hinders this easy maneuver in conjunction with an eclipse is the time constraints. This 
suggests that the ideal platform to use MACS be in tandem with a coronagraph. 


180 



3. The instrumental stray light could be reduced to a great extent only through superior 
optics. A fair estimate of the distribution of the instrumental stray light could also be 
estimated by placing many more fibers in the shadow of the moon and outside the 
field of view in the focal plane of the telescope as depicted in figure (6.7). 



Figure (6.7). The extra fibers placed in the shadow of 
the moon and outside the field of view in order to 
record the instrumental stray light distribution. 


4. Another important deficiency in MACS in its present state is in its inability to 
simultaneously visually image the focal plane during the observations. This prevented 
us from knowing whether the fibers stayed in focus at their intended locations. 
Deviations from the intended positions can arise from slight movements of the 
telescope due to wind, and also due to deficiencies in the tracking mechanism. This 
understanding could be easily satisfied by attaching a CCD-based video camera to the 
viewfinder of the telescope that in turn is co-aligned with the telescope. The live 
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video image of the sun through a sun filter projected on a TV screen, at least from one 
hour prior to the first contact, can be placed on a suitable 2_D grid and its course 
followed for tracking purposes. The bright pixels prior to the first contact should 
appear as dark between the second and the third contacts with the sun filter removed. 
In order to achieve this goal, the viewfinder and the telescope should be perfectly 
aligned. Although this enables us to monitor the tracking efficiency it does not prove 
the centering of the sun on the image plane of the telescope. To accomplish this the 
following maneuver could be followed. 

(a) . First center the sun’s image through the sun filter on a reticle eyepiece. The reticle 

pattern in the Meade brand is a double cross line with two concentric circles. At this 
stage if the viewfinder is perfectly aligned with the telescope then the sun’s image 
seen through a video recorder attached to the viewfinder should display the sun’s 
image at the center of the detector. Centering on the reticle eyepiece and making 
adjustments to the telescope’s orientations can be carried out alternatively until the 
image remains centered on the viewfinder and should satisfy the tracking efficiency. 

(b) . Then attach a camera with the lenses demounted for prime focus photography of 
the sun through a sun filter at the focal plane of the telescope. For the Meade brand of 
12-inch Schmidt_Cassegrain telescopes one end of the #62-T adapter connects to the 
F/6.3 focal reducer and the other end threads into a universal T-adapter with bayonet- 
mount unique to the brand of 35mm camera and focus the image. This should 
facilitate the focusing for the experiment since the fiber optic plate too will be located 
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at the prime focus of the telescope during observations. If step (a) was successfully 
performed then the image should be centered too. 

However the above features do not yet allow for the simultaneous imaging of the 
coronal image formed on the front surface of the fiber optic plate. The following optical 
feature as shown in figure (6.8) may be contemplated for MACS in order to obtain 
simultaneous imaging of the image formed on the focal plane. In figure (6.8) the fibers 
are embedded on an optical flat (well polished reflecting surface facing the telescope) 
along with a very thin high transmission glass surface (AB) inclined at 45° to the optical 
axis and with a semi-reflecting coating on the side facing the fiber optic plate. The 
reflected image from (AB) is a copy of the image formed on the fiber optic plate with 
dark dots marking the light transmitted through the fibers. Such an image will identify the 
actual location observed by each of the individual fibers on the corona. This is ideal to be 
used in attachment to a coronagarph. This is because any additional photons lost in the 
transmission through the additional optical component could be compensated by a longer 
integration time. 
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Figure (6.8). A schematic diagram showing the introduction of a very 
thin high transmission glass surface (AB) inclined at 45° to the optical 
axis with a semi-reflecting coating on the side facing the fiber optic plate 
along with the fibers embedded on an optical flat. This alteration will 
enable the simultaneous imaging of the image formed on the fiber optic 
plate. 


The proposal (NRA:99-OSS-01) to build a coronagraph to accommodate MACS 
at a cost of ~$ 500,000.0 has been approved by the NASA’s Office of Space Science. 
NASA’s Goddard Space Flight Center and The Catholic University of America will 
jointly build this instalment over the next three years. Possibilities for observing sites 
include Sacramento Peak, Flawaii and the Canaries. 
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Chapter 7 


CONCLUSION 


7.1 Summary of this thesis 

The methodology presented in this thesis is unique in its character: it allows for 
the simultaneous and global measurement of the temperature and the solar wind velocity 
in the solar corona. This line of thinking originated from Cram’s (1976) theoretical 
formulation for the formation of the K-corona and its properties. The physical properties 
that were taken into consideration in the formulation of the theory included the effects 
due to thermal Doppler broadening, limb darkening effects, angular dependence of the 
Thomson scattering and the cumulative contribution from different positions along the 
line of sight. The properties that were derived were the existence of temperature sensitive 
anti-nodes in the K-coronal spectrum and the remarkable independence of the 
temperature insensitive nodes with height above the solar limb. It was suggested that 
these properties could be exploited to measure the coronal temperature from the K- 
coronal intensity profile. Ichimoto et al. (1996) demonstrated the feasibility of Cram’s 
method by analyzing the K-coronal spectrum obtained using a slit spectrograph during 
the total solar eclipse of 3 November 1994 in Putre, Chile. 
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Ill this thesis Cram’s theory was extended by including the effects of coronal 
expansion in the formation of the K-corona. While confirming the properties derived by 
Cram, an additional property was discovered with the introduction of the solar wind. This 
allows for the simultaneous determination of both the temperature and the wind 
velocity from the same K-coronal spectral measurement. An overview of the theory and 
of the method can be found in chapters-1 and 2. The dependence of the theoretical model 
on the selection of the electron number density model, different temperature and wind 
profiles and the numerical procedures are discussed in chapter-5. 

In this thesis, instead of a slit spectrograph a fiber optic spectrograpgh was 
designed whose features allowed for global determination of the K-coronal spectra all 
around the solar limb and at different heights. The instrument is discussed in chapter-3. 
Combining the theoretical prediction on the feasibility to determine both the temperature 
and the wind velocity from a single K-coronal spectrum and the instruments ability to 
globally measure the K-coronal spectra, the instrument was named MACS for Multi 
Aperture Coronal Spectrometer. The very first effort of its kind was tested in conjunction 
with the total solar eclipse of 1 1 August 1999 in Elazig, Turkey. The observational results 
are discussed in chapter-4. We were successful in determining the electron temperature at 
various points in the inner corona. 

However, an important (and unwelcome) aspect of our observations concerns the 
degradation of the data due to instrumental scattering. Because of unexpectedly high 
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levels of scattering, we were not successful in obtaining reliable estimates of wind 
velocity. Although remedial measures were taken in expectation of possible instrumental 
scattering, the actual magnitude of the scattered light was not realized. Our efforts to 
reduce scattering were compounded by the inability to replicate the observing conditions 
in a laboratory and the lack of funds to obtain high quality optical components. One of 
the vital ingredient required to minimize instrumental scattering is careful selection of 
high quality optical components. Another important feature for the success of MACS 
would be to perform an absolute wavelength calibration in the laboratory. The extent of 
the instrumental scattering and some of the remedial measures are discussed in chapter-6. 

Despite the lack of success in determining the wind velocity, our success in 
determining the coronal temperature using MACS has already been recognized as a 
successful proof-of-concept by a funding agency. A proposal to mate MACS concept 
with a coronagraph was submitted to NASA and it received favourable review. As a 
result, funds have been awarded for an upgraded version of our experiment in which the 
scattered light problem should be much reduced. It is hoped that the upgraded version of 
MACS will be ready for the solar eclipse of 21 June 2001 in Africa. At this eclipse 
MACS is to be complimented with a filter based system, where filters would be centered 
at the temperature and wind sensitive wavelength positions. The concept of MACS has 
also been recommeded as an instrument in the proposal for the instrumental package that 
would fly aboard the Solar Probe. A brief review of the Solar Probe is given in section 
(7.2). 
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7.2 Using a spacecraft platform 


A space platform envisaged for implementing the theoretical concept of MACS is 
the unmanned Solar Probe Spacecraft that will approach the sun at four solar radii at its 
perihelion. In this proposal a series of filters centered on the temperature and the wind 
sensitive wavelength positions is proposed. This would allow for both temperature and 
wind velocity measurements. Figure (7.1) shows the parabolic path of Solar Probe. 



(r.Q t - polar coordioatts of the parabola w 
ft -angld P \Fi 
h -angler \ PI ■ 

V ■ angle! < ■ I* \i 


ib=angle(EFP| 


Y 


P - solar probe 

PCD - parabolic orbit of the probe 
PA - tangent at P 

F - center of the sun & focus of parabola 
FG - perpendicular to PG 
PB, PE - tangents to the sun 


Detector direction 


Figure (7.1). The path envisaged for the unmanned Solar Probe 
Spacecraft. At the perihelion (C) and the poles the spacecraft will be four 
and eight solar radii, respectively, from the center of the sun. The 
position of the spacecraft (P) is given by (r, 0 ) with respect to the center 
of the sun. In this picture the parabolic path, the detector plane and the 
center of the sun are in the same plane. 
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The equation of the parabolic orbit of the Solar Probe as depicted in figure (7.1) in polar 


coordinates is given by equation (7.1). 


r(0) = 


2.0 xFC 
1 - cos(0) 


8.0R 


1 - cos(0) 

where R is the solar radius & FC = 4R 


(7.1) 


The angle OC made by the tangent at P with the x-axis is given by equation (7.2). 


a = tan 


-i 


1 - cos(0) 
sin(0) 


(7.2) 


The derivation of equation (7.2) is as follows. 

y = r sin(0) —> dy = r cos(0)d0 + sin(0)dr 
x = rcos(0) — » dx = -rsin(0)d0 + eos(0)dr 
dy rcos(0)d0 + sin(0)dr _ reos(0) + sin(0)(dr/d0) 


dx - r sin(0)d0 + cos(0)dr -rsin(0) + eos(0)(dr/d0) 

dr 


from equation (1) 


2xFCxsin(0) „ . 2xFC 

& using r = 


(l - cos(0))^ 


l-cos(0) 


, , . dy l-cos(0) 

tan(a) = — = — 

dx sin(0) 


From the triangle APE the relationship between a, (3 and 0 is given by equation (7.3). 


0 = a + (3^(3 = 0- a 


(7.3) 
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Supposing the detector direction is 7 degrees from the tangent to the orbit at position (P) 


then the perpendicular distance p from the center of the sun to the line of sight of the 
detector is given by equation (7.4). 

p = rsin(P + y) (7-4) 

At this instance the distance between the spacecraft and the intersection between the line 
of sight and the normal from the center of the sun is given by equation (7.5). 

PG = r cos(p + 7 ) (7*5) 

In order to prevent the detector being pointed directly at the sun the angle 7 will have to 
satisfy the following condition given by equation (7.6). 


7 _ > -(P - P mi „ ) measured in the (-) direction 
7 + < 2 tc - (P + P Illin ) measured in the (+) direction | 

. n • -ifR' 

where p, nin =sin 


(7.6) 


In determining the theoretical K-coronal spectrum observed by the Solar Probe the only 
change occurs in the x-integration (along the line of sight) as shown in equation (7.7). 


Jdx — > Jdx 

-PG 


(7.7) 
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The above integration may be performed the following way as shown in equation (7.8). 
The first and the second integrations are performed using trapezoidal composite rule and 
Laguerre integration, respectively. 


-PG 


J f (x)dx = j f (x)dx + 1 f (x)dx 

-PG 0 

0 +°° 

= J f (x)dx + J e“ x (f (x)e +x )dx 


-PG 
( 20 


h 


\ 20 


2x£f(z.)-f(z 0 )-f(z 28 ) +^ w i e+x, f(x i ) 


i=0 


I where h = & z . = -PG + i x h 

20 


(7.8) 


Figure (7.2) and figure (7.3) show for solar wind velocities 0.0 and 400.0 knt/sec, 
respectively, the theoretical K-coronal intensity variation for various detector elevations 
with respect to the tangent when the spacecraft is positioned at the pole (P). Figure (7.4) 
shows the variation of the wind-sensitive intensity ratio against the detector tilt y for 
various wind velocities at a constant temperature of 1 .0 MK and the spacecraft positioned 
at the pole. Figure (7.5) shows the variation of the temperature- sensitive intensity ratio 
against the detector tilt y for various temepratures at a constant wind velocity of 400.0 
km/sec and the spacecraft positioned at the pole. Figures (7.6) and (7.7) show the 
variation of the temperature-sensitive intensity ratio against temperature for various wind 
velocities and the variation of the wind- sensitive intensity ratio against wind velocitiy for 
various temperatures, respectively, for a given detector tilt angle y and the spacecraft 
positioned at the pole. 
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Figure (7.2) The theoretical K-coronal intensity variation 
for various detector elevations for 0.0 km/sec wind velocity 
and the spacecraft positioned at the pole (P). 


Figure (7.3). The theoretical K-coronal intensity variation 
for various detector elevations for 400.0 kni/sec wind 
velocity and the spacecraft positioned at the pole (P). 
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( 41 uO}/i(385Gj VS Temperature for 50 angstrom bandpass 
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Figure (7.6). Temperature-sensitive intensity ratioVS 
Temperature for various wind velocities at y = 0.0 
degrees and the spacecraft positioned at the pole. 


I(4233)/I(3987) VS Wind for 50 angstrom banc 



+00 
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Figure (7.7). Wind-sensitive intensity ratioVS Wind 
velocity for various temperatures at 7 = 0.0 degrees and 
the spacecraft positioned at the pole. 
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These plots show the promise behind the theoretical concept of MACS being 
employed from a space platform deep inside the solar corona for the measurement of the 
coronal temperature and the solar wind velocity thus avoiding the F-corona and the use of 
a coronagraph. Even the spectrograph could be avoided with the use of optical filters 
centered at the temperature and wind sensitive wavelengths. The temperature and wind 
sensitive intensity ratios remain almost the same even at fifty angstrom bandpass. Figures 
(7.8) and (7.9) show the temperature and wind sensitive intensity ratios, respectively, for 
the use of fiters centered at those wavelengths with a fifty angstrom bandpass for line of 
sight integration from ± °° at 1.1 solar radii. This could be verified by the comparison of 
figures (7.8) and (7.9) with figures (2.3) and (2.20), respectively. 


i(41 00}/I(385G) VS Temperature far 50 angstrom bandpass 



Figure (7.8). Temperature-sensitive intensity ratioVS 
Temperature for various wind velocities with intensities 
measured through filters of fifty angstrom bandpass 
centered at those wavelengths for line of sight at 1.1 solar 
radii. 
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l(4233)/i(3987) VS Wind for 50 angstrom bandpass 



' ; 

~i | i r— 

J T , 

^ i 1 j 1 i 1 i 1 : 


I i 

i 

1 

1 

1 

Temperature MK / 


I 

1 

1 

/ i 


i 

i 

1 

1 

1 

i o.5gstioooo i I 

1,270 

— i 

1 

1 

1 ! — 


: | 

1 

i 

1 

1 

1 

1 

l" l 


i 

i 

1 

l 

1 

1 

1 

y' | ^ ; _ 

1.260 


l 

1 

S 1 i 



1 

1 

s' 

s’ | 

,t " 1 


- 

1 

1 

1 

1 

J" 1 

; 

1.250 

E- ! 

1 

1 


! J ^ 


: i 

y'\ 




- | 

s' ,J'- 

J,- 

i ! i 

1.240 

r | 


i 

! ; - 


Z yv 


l 

i 

1 1 


i 

i 

i 

l 

i 

j ! : 

1.230 

i 

i i 1 

1 

_i 1 i i 

, . ! . 1 . ... 4 


400 600 

Wind km/sec 


800 1000 
rim in solrjr radius 1 .1000000 


Figure (7.9). Wind-sensitive intensity ratioVS Wind 
velocities for various temperatures with intensities 
measured through filters of fifty angstrom bandpass 
centered at those wavelengths for line of sight at 1.1 solar 
radii. 


In conclusion, with the approval of funding for a coronagraph for MACS the 
future is very optimistic as a tool for simultaneous and global coronal temperature and 
wind measurements. There is also good reason to hope that the concept can be applied as 
an instrument on a space-based platform. 
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Appendix A 


THE THEORY OF FORMATION OF THE K-CORONA 


A.l The general theory on scattering of radiation by free electrons 



Figure (A.l). The general coordinates where f is a field 
/ 

point and is a source point. 
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Consider a system whose charges and currents are varying in time. There is no 
loss in generality by restricting our considerations to potentials, field and radiation from 
localized systems that vary sinusoidally in time. 


The time dependent charge and current densities could be written, respectively, as given 
by equation (A.l) and equation (A.2), based on coordinates shown in figure (A.l). 


p(C,t) = p(r) e KOt 


J(r',t) = J(r) e ' 


icot 


(A.l) 

(A.2) 


And the continuity equation is given by equation (A.3). 


V.J + ^- = V.J + (-icop) = 0 

at 


(A.3) 


The scalar and vector potentials for the charge distributions are given, respectively, by 
equation (A.4) and equation (A.5). 


V(r) = 1 

47T£ 


1 .p(r> ik fc-^ 


dV 


r-r 


p ftr m ,kr_r 

A(r) = — J d*X 


4 n 


r-r 


(A.4) 


(A.5) 
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Ill the far radiation zone, where r » □ and r » r', the denominator Ir- bl can be 
considered to be independent of b although the argument of the complex exponential is 
not. Thus, equation (A.5) could be written as equation (A.6). 


p Jkr ' , 

A(r) = — - — Ir) <,3^ 

4k r 


(A.6) 


Also if the source dimension b is small compared to the wavelength, then k b « 1 which 
justifies writing equation (A.6) as equation (A.7). 


A(r) = 


V ikl - ,/r.r'^ 1 


4n r | 1! 


-JJ(0 


v r j 


d 3 L 


(A.7) 


The above is tme since the electron diameter of 5x 10 -15 m is very small compared to the 
wavelength of visible light of 5x 10 07 m. Therefore the expansion of the exponential 
term in the integrand of equation (A.6) is justified by the summation given equation 

(A.8). 


ik r.r / r 



(A.8) 


The scattering of electromagnetic radiation by systems whose physical 
dimensions are small compared with the wavelength of the wave being scattered, it is 
reasonable to assume that the incident radiation as inducing electric and magnetic 
multipoles and these to oscillate in definite phase relationship with the incident wave and 
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also radiate energy in directions different from the direction of the incident wave. For the 
case where the wavelength of the incident radiation is very long compared to the size of 
the scatterer, only the lowest multipoles, usually the electric and magnetic dipoles ( 1 = 0 ) 
are important. Therefore for electric dipole radiation (1=0) equation (A.7) reduces to 
equation (A.9). 


p ikr 

A(r) = — JJ(L)d 3 L 

471 r 

Consider the following mathematical operation 


(A.9) 


V'.(x'J) = (VV)J + x'(V / .J) 

= J x + x'(V'j) 

•••Jx = V'.(x'J)-x'(V'.J) 

iJxd 3 ?; =|V'.(x'j)d 3 r -Jx'(V'.j)d 3 r 

= |x'J.ds -J x'(V'j)d 3 r 
Using equation(1.3) 

= |x'J.ds +Jx'— d 3 r 

at 

Using equation(l.l)and integrating over a large enough volume 
= -icojVpd 3 r 

and for all components of J 

JJ(r )d*V =-ico|r'p(r )d^r =-icop 

where dipole moment p =Jr'p(r )d^r 


(A.10) 
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Equation (A.ll) can be written from equation (A.9) and equation (A.10). 

—icon Jkr 

A(r) = 5-- — p (A.ll) 

4k r 

Then, using equation (A.ll) the magnetic field is given by equation (A.12). 


(A.12) 

From Maxwell’s equations for outside the source (J=0) the magnetic and the electric 
fields are related by equation (A.13). 

(A. 13) 


From equation (A.13) the electric field is given by equation (A.14). 

. 2 

E(r) = — VxB(r) (A.14) 

co 
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And using the following vector identities 


A x (B x C) = (A.C)B - (A.B)C 
V x (A x B) = A(V.B) - B(V x A) + (B.V)A - (A.V)B 

and the relations 

2 

V.r=— , V.r = 3 and p ^ p(r) 
r 

on equation (A.14) the expression for the electric field is given by equation (A.15). 



ckcou 

o 

47t 



1 


4718 

o 




(A.15) 


In the radiation zone kr»l where k = co/c the equations (A.15) and (A.12) reduce to 
equation (A. 16). 



(A.16) 
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The Poynting vector (average power radiated per unit area) is given by equation (A.17). 


* 

(s) = — — — Watts/m 2 (A.17) 

2p 0 

Using equation (A. 16) and the vector identity 

(A x B )x C = (A.C)B - (B.C)A and I).(I) x E) = 0 and k=co/c the Poynting vector 

given by equation (A.17) reduces to equation (A.18). 

4 

(0 LI 

O 

,,, 2 2 

32 k r c 

Then the total power radiated through an area da is given by equation (A.19). 


<s> 


L \ 2 
(rxp)' 


Watts/m' 


(A.18) 


dL = (s).da = (s).rr 2 dQ Watts 


(A.19) 


Substituting equation (A.18) in (A.19) gives equation (A.20). 


CflV / \ C0 4 p , , 

dL = 5_(rxpj 2 r.rr 2 dQ = 9_(rxp) 2 dQ Watts 

32Tt 2 r 2 c 32 k 2 c 


(A.20) 


Now for an oscillating electric dipole the induced moment p is given by equation (A.21). 


P = 


Bo 6 


— icot 


(A.21) 
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Differentiating equation (A.21) twice with respect to time gives equation (A.22). 


p = to 2 p 

Substituting equation (A.22) in (A.20) gives equation (A.23). 

dL = (rxp) 2 dQ Watts 
32ti 2 c W 

For a single electron the induced moment p is given by 
p(t) = er^(t) 

where is a source point with respect to its origin. 

From the equation of motion in the non-relativistic case 

mb = e(E + v x B) = eE 

where m is the mass of the electron. 

Substituting equation (A.25) in (A.24) gives equation (A.26). 

.. e 2 E 

P = — — 

-* m 


(A.22) 


(A.23) 


(A.24) 


(A.25) 


(A.26) 
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Again substituting equation (2.26) in (2.23) gives equation (A.27). 


( ,,2 \ 


dL 


V "1C , 


^(fxE) 2 dn 


(A.27) 


32ti 


Diving and multiplying equation (A.27) by 67ie o gives equation (A.28). 


dL = — o T e 0 c(rxE) 2 dQ 
16k 1 * 


(A.28) 


where O is the Thomson scattering cross section given by equation (A.29). 


( „2 \ 


0 rji — 


6k£ 


o 


v mc J 


(A.29) 


Therefore from equation (A.28) the power radiated per unit solid angle is given by 
equation (A.30). 


— = — o T (rxE) 2 Watts/solid angle 
dQ 1671 "* 


(A.30) 
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A.2 Total power radiated per solid angle for the cases where the electric vector is 
aligned parallel and perpendicular to the scattering plane 

Consider the case where the electric vector is parallel to the scattering plane. In 
figure (A.2) n in and Il out are unit vectors in the directions of the incident and the 
scattered radiation and they both lie on the scattering plane and 0 is the scattering angle. 

/V Sc 

E /; is the unit vector of the component of the electric vector parallel to the scattering 

* Sc 

plane and E is the unit vector of the component of the electric vector perpendicular to 
the scattering plane. 


7 

^SC 



Figure (A.2). The general coordinates of the incident and scattered 
radiation and the orientations of the parallel and the perpendicular 
components of the electric vectors with respect to the scattering plane. 
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A A 

The direction cosines of n in and Il out are given by equation (A.31) and equation (A.32), 


respectively. 


f 

f 

K N 


( 

K ' 

\ 

-sin 

0- 

2 , 

,-cos 

0- 

2 , 

,0 

V 

v 

J 


V 

J 

/ 


Sout - (IjOjO) 


(cos 0,-sin0,O) 


(A.31) 

(A.32) 


From figure (A.2) the dot product ofn. and n gives equation (A.33) 


Sin'Sout _ Sin Soutl cos 


«(e) 


(A.33) 


and the square of its cross product gives equation (A.34). 


(Sin x Uout) — Inin I ISoutl s ' n (®) 


(A.34) 


a Sc A Sc 

The direction cosines for the electric vectors E and E are, respectively, given by 
equation (A.35) and equation (A.36). 


A £ c 

E /; = (-sin(7t-0),cos(7t-0),O)= (-sinO,-cos0,O) 


(A.35) 


E* C = (0,0,1) 


(A.36) 
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The sum of the cross products of equation (A.32) with (A.34) and (A.35) is gives 
equation (A.37). 



1 j i 


l J i 

/v ^ Sc ^ ^ Sc 




n ,xE„ +n ,xE = 

-►out — ►// -*• out — 

1 0 0 

+ 

10 0 


- sin 0 - cos 0 0 


0 0 1 


= -cos0k + j 


(A.37) 


/\ /\ 

From figure (A.2) using the relation T — H out and equation (A.37) the following 
relations are obtained as given by equation (A.38) and equation (A.39). 


. ^ Sc o 

rxE.Ef; 


= (Ef; ) 2 sin 2 (90 - 0) = (Ef ) 2 cos 2 0 


(A.38) 


^rxE^E? J =(E* C ) 2 ^ A - 39 ^ 

Now substituting the equation (A.38) and equation (A.39) in equation (A.30) give 
equation (A.40) and equation (A.41), respectively. 


dQ 


\Sc 

J H 



(8 o c(E^) 2 )cos 2 0 


(A.40) 


dL 

dQ 


xSc 

A 



(e 0 c(E?) 2 ) 


(A.41) 
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But for isotropic radiation equation (A.42) is true. 




(A.42) 


Substituting equation (A.42) in (A.40) and (A.41) give equation (A.43) and (A.44), 
respectively. 


dL 

dQ 


\Sc 

jH 


16tt 


ice 0 E 2 

V Z J 


cos 2 0 


(A.43) 


dL 

dQ 


.Sc 

J - L 




O 


E 2 


\ 

/ 


(A.44) 
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Using the Poynting vector S = — ExB the energy flux density transported by the plane 

monochromatic electromagnetic wave depicted in figure (A.3) is given by equation 

(A.45). 

S = ce 0 E„ cos 2 (kx-o)t + S)i (A.45) 
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and averaging equation (A.45) over a cycle gives equation (A.46). 



(A.46) 


From equation (A.46) the average power per unit area transported by an electromagnetic 
wave is called the intensity and is given by equation (A.47). 


C = (S) = -ce 0 E 0 2 


(A.47) 


Substituting equation (A.47) in equation (A.43) and equation (A.44) give equation (A.48) 
and equation (A.49), respectively. 


I Sc 

1 // = 


1 dL ^ Sc 


dQ 


Jn 


16 k 


o T I 0 cos 2 0 Joules/sec.steradian. A .electron 


(A.48) 


Sc 


j Sc _ 


dL v 


3 0 

o T I„ Joules/sec.steradian. A .electron 


1671 


(A.49) 


where I Q is the incident radiation. 
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Consider a radial plane inclined at an angle a to the scattering plane, as shown in figure 
(A.4). The variation of intensity through an angle (3 is given by equation (A.50). 

I = I M cos 2 (3 (A.50) 

where I M is the maximum intensity. 



Figure (A.4). Coordinates for the rotational transformation of the 
scattering plane through an angle a about the line of sight. 
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From figure (A.4) and using equation (A.50) equation (A.51) and equation (A.52) are 
obtained for the intensity parallel and perpendicular to the radial plane, respectively. 


If = I „ cos a + I f cos ' 


\ 


71 

— a | 

v 2 J 


= lf c cos 2 a + lf c sin 2 a 


(A.51) 


If = I/ c cos^ a + 1,; cos 


V 

= lf, c sin 2 oc + lf c cos 2 a 


71 

— a 
2 


J 


(A.52) 


Substituting equation (A.48) and equation (A.49) in equation (A.51) and equation (A.52) 
give equation (A.53) and equation (A.54), respectively. 


| Ra 
1 // “ 

f dL 1 

[dQj 

Ra ^ 

= - 1 — o T I 0 (sin 2 a + cos 2 0cos 2 a) 
u 16jc 

= lf c sin 2 a + if/ cos 2 a 

= Qf (a,0)l o Joules/sec.steradian.A°. electron 


NH 

1- 53 
P 

III 

( dL ) 

[dQj 

Ra ^ 

= — : — o T I 0 (cos 2 a + cos 2 ©sin 2 a) 

j_ 16 tc 

= lf c cos 2 a + lf c sin 2 a 

= Qf (a,0)l o Joules/scc.steradian. A ".electron 


(A.53) 


(A.54) 
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Using the scattering theories derived above the same could be applied for the 
scattering of photospheric radiation by the coronal electrons. Figure (A.5) shows the 
corresponding geometrical configurations depicted by figures (A.2) and (A.4) for the case 
of photospheric radiation being scattered by the coronal electrons to an observer on 
earth. 



E electron 
C center of Sun 
S source point 
O observer 
SEO scattering plane 
CEO radial plane 


Figure (A.5). Geometrical configuration for the scattering of photospheric 
radiation by coronal electrons. 
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Figure (A.6). Geometrical configuration for the formation of the K- Corona. 
Photospheric radiation emitted from a point S on the Sun is scattered from an 
electron E towards an observer O. The solar wind on the electron is radial and 
blows in the direction CE. 
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A.3 Contribution due to Compton scattering 


In the rest frame of the electron the scattering of the photospheric light incident on 
the electron is coherent. The red shift due to the Compton effect is very negligible. To 
quantify the contribution due to Compton effect, consider radiation of wavelength X 
incident on a stationary electron, as shown in figure (A.7). Let the radiation scattered off 
the electron be of wavelength X' and the velocity of the electron now be v. The 
scattering angles for the radiation and the electron are, 0 and (p , respectively. 


Scattered radiation 



Figure (A.7). Scattering of radiation off an electron also known as Compton 
scattering. 
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From conservation of energy in figure (A.7) equation (A.55) is obtained. 


f 


\ 


he 

T 


he 

— + me 

X 1 


2 



1 


J 


(A.55) 


From figure (A.7) using conservation of momentum in the X the Y directions, 
respectively, equation (A.56) and equation (A.57) are obtained. 


h h 

— = — COS0 + 

X X' 


111V 



cos 9 


0 = — — sin 0 + 

u 


mv 



sintp 


where 


h = Planck cons tan t = 6.63 x 10“ 34 J.s 
m = electron mass = 9.11x10 31 Kg 
c = speed of light = 3 x10 s ms 1 


(A.56) 


(A.57) 


Eliminating cp and v from equation (A.55), (A.56) and (A.57), the wavelength shift is 
given by equation (A.58). 


AX = X'-X = -^-(l-cos0) (A.58) 

me 

= 0.0243A* (l - cosG) < 0.0486A 0 
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A.4 Expression for the total K-Coronal scattered intensity 


From figure (A.6) for the calculation of the total K-Coronal scattered intensity in the rest 
frame of the electron, first, it is necessary to consider an electron velocity distribution over a 
volume element at P in order to detennine the number density. Consider a Maxwellian velocity 
distribution for the coronal electrons. Then the number density at the point P in the velocity 
interval (u, U + du ) is given by equation (A.59). 


f,(u)=N,(p) 




where 


q(P) = mean electron thermal velocity 



k = Boltzmann constant = 1.38 x 10 23 JK 1 
m 3 = electron mass = 9.11 x 10“ 31 kg 
T e in 10" Kelvins 


5508^/l7 kms 1 


(A.59) 
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However, in the rest frame of the observer, the scattered radiation off a moving 
electron will be altered in wavelength from the monochromatic radiation incident on that 
electron. Consider a coronal electron with velocity u subjected to radial solar wind 

velocity w in a coordinate system where the x-axis bisects the supplement of the 
scattering angle, as shown in figure (A.8). 


Jf 

s 



Figure (A.8). Construction to obtain an expression for the scattered 
intensity in the rest frame of the observer. Consider a coronal electron with 
velocity u subjected to radial solar wind velocity w in a coordinate system 
where the x-axis bisects the supplement of the scattering angle. 
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From figure (A.8) the direction cosines of the unit vectors n jn and n out are then given by 
equation (A.60) and equation (A.61). 


n. n = (-cosy, sin y,0) 

(A.60) 

Sout = (cosy, sin y,0) 

(A.61) 


And the net velocity vector V is given by equation (A.62). 
V = (u x + w x ,u y + w ,u, + wj 


(A.62) 


From Doppler effect, in the rest frame of the observer, the observed radiation has 
wavelength given by equation (A.63). 


1 = f \ 

* observed * scattered 



n 

.V 

1- 

out — ► 

c 


n 

.V 

1 + 

-♦ out — ► 


(A.63) 


The relationship between X observed and ^ photosphere is given by equation (A.64). 


X =X 

* scattered photosphere 


1 + 


n .V 


n .V 


1 - 


(A.64) 
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From equation (A.63) and equation (A.64) eliminating X scattered gives equation (A.65). 


( n . .V ^ f 


\ — \ 

* observed w photosphere 


w n. .V ^ 


1 


► out 


n ,.V 

J _|_ -►out — ► 




1 + 


n .V 


(A.65) 


Since the speed of light c» V equation (A.65) can be reduced to equation (A.66). 


— X, 

photosphere 

f (n. -n ).V 

j 1 Wm -►out / — ► 

c 

(n. -n ).V 


j Win -►out / — ► 


l c 


^2 


(A.66) 


By Taylor expansion of equation (A.66) it could be reduced to equation (A.67). 


X 


observed 


{ 1 + 




i 





photosphere 




v 

( 


= X 


photosphere 


1 + 


v 



j 


\ 


c 


J 


(A.67) 
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From equation (A.60), (A.61) and (A.62) equation (A.68) could be obtained. 


(a 


in 



2cosy^ 


'"u + YV A 




(A.68) 


From equation (A.67) and equation (A.68) and using the notations used in figure (A.8) 
equation (A.69) could be obtained, which relates between X obseved and X photosphere . 


x=r 


1 - 2cosy 


u + w 




where 


observed 




photosphere 


(A.69) 


Equation (A.79) is obtained by rearranging equation (A.69). 


( 2cosyA,' N 

U x - 

r 

2cosy ' 

1 w x 

-X 

= 0 

l c J 

r 

_ 

l c ' J 

. 



(A.70) 


Equation (A.70) satisfies the condition that scattering is coherent in the rest frame of the 
electron. 
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Now, from figure (A.8) and equation (A.53), (A.54), (A.59), (A.70) the intensity 
scattered from an electron at point P is given by equation (A.71). 


If (A, nj+ If (A, n ou t ) = l(A', n o ut )x (q«- (a, 0) + Qf (a, 0))x ND 

where 


(M) -- LI c 


A'[ i_2cosy w A. ] du x du du z 
c 


Tlie expression for ND in equation (A.71) reduces to equation (A.72). 


ND = ^ije ^uje’-du, 

IvTtqJ - 


2bA'u. 


f 2bA'u x ^ (.,( 2bw x ^ 

O - A 1 -A 

c c 




= -jNAix(V^(|)x(VS'q)xexp ‘ j 

(vTtq) f 2 bqA ^ 


N e (P) 

j2rJnAb] 


■exp - 


A-Afl- 2 ^ 

L 1_ 

(2Ab) 


, * qA r _ _ (k-Q 

where A = and b = cosy = cos 

c 2 






From equation (A.71) and (A.72) the scattered intensity from point P is given by equation 


(A.73). 


lS a (^n oiit ) = l(r,n in )xQ“ a (a,0)x^^Iexp| 


where O = (//, _L) 


f ( 2bw. 


x-r 


i- 


2Ab 


(A.73) 


To obtain an expression for the total observed scattered intensity, from figure (A.6), the 
equation (A.73) needs to be integrated over the following parameters. 


1 . All wavelengths from each point on the photosphere 


JdV 


(A.74) 


2. From all points on the photosphere 


271 0 ) 271 1 

Jd(pJsincodco= Jdcp Jdcosco 


(A. 75) 


3. From all points along the line of sight 


Jdx (A.76) 
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From figure (A.6) and figure (A.8) the radial solar wind velocity has components given 
by equation (A.77). 

w = (— cos cocos y, cos cosin y, sin co) = (w x , w y , w z ) (A.77) 

Using equation (A.74), (A.76), (A.76) and (A.77) in (A.73) gives the expression for the 
total observed intensity for a given observed wavelength X at a given line of sight 
distance p from the center of the Sun (see figure (A.6)), as given by equation (A.78). 


+oo2n 1 o° 

Iq 11 (X, p) = J J J J dA,' d cos co dcp dx x 

— 00 0 cos co* 0 


N,(r)xQ"*(a,0)x 


X-X'( 


2y[nAb 


l(V,n. )exp 


, , 2b 2 cos cow(r) radia , 


2Ab 


L V J 

where O = (//,!) 

// parallel to the radial plane 
± perpendicular to the radial plane 


(A.78) 
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A.5 Expressions for the dependent variables in equation (A.78) in terms of the 
independent variables 

Figure (A. 9) is a highlight from figure (A.6) of the triangles formed by SEW 
(scattering plane), CEF (radial plane) and CSE. 



Figure (A.9). Highlighted map of the radial and the scattering planes of figure 
(A.6). 

From figure (A.9) equation (A.79) and equation (A.80) can be obtained. 

ES = |ES|(sino)coscp, sin cosin <p, cos co) (A.79) 

EF = |EF|(0, sin x, cos x) (A.80) 
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And the dot product of equation (A.79) and equation (A.80) gives equation (A.81). 


ES.EF = E S EF (sin co sin cp sin y + cos cocosy) 
= |ES||EF|cos(7r-0) 


(A.81) 


From equation (A.81) the expression for 0 is given by equation (A.82). 


0 = k - cos 1 (sin cosin cpsin y + cos cocos %) 


(A.82) 


Also from figure (A.9) the expression for angle y is given by equation (A.83). 


x 

cos x = 

r 


(A.83) 


From figure (A.6) and measuring distances in solar radius the expression for angle co is 
given by equation (A.84). 


. . HC 1 

sin co = = — 

EC r 


(A.84) 
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From figure (A.6) consider the spherical triangle formed by EDJI, as shown in figure 


(A.10). 


I 



Figure (A.10). Highlight of the spherical triangle formed by EDJI in figure 
(A.6). 


From figure (A.10) the relationship between the angles is given by equation (A.85). 


sin a _ sincp 
sinco sin(7i - ©) 


(A.85) 
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From equation (A.85) the expression for angle a is given by equation (A.86). 


. sin co sin cp 

a = sin — 4- 

sin(jt - 0) 


Figure (A.ll) shows the triangle ECS from figure (A.6). 



Figure (A.ll). Highlight of the triangle ESC of figure (A.6). 

From figure (A.ll) the relationship between the angles 0 and co is given by equation 

(A.87). 

1 r 

sin co sin(rc-0) (A.87) 

.*. sin0 = rsinco 
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A.6 Expression for the incident solar intensity on the coronal electrons 

In order to evaluate equation (A.78) the incident solar intensity l(A/,n. n ) need to 
be known. The following expressions are from Astrophysical Quantities by Allen. Let 
I (A/,0) be the intensity of the solar continuum at an angle 0 from the center of the disk 
and l(A/,0) the continuum intensity at the center of the disk. The results may be fitted by 
the expression given by equation (A.88). 


i(v,e) 

l(V,oj 


= l-u 2 - v 2 + u 2 cos0 + v 2 cos 2 0 


where u 2 and v 2 are limb darkening constants 


(A.88) 


Or less accurately equation (A.88) can be written as equation (A.89). 

i(r,e) , 

— 7 1 — r = 1 — U , + U, COS0 

i(r,o) 1 1 

where cos0 is the heliocentric angle shown in figure (A.6). 


(A.89) 
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And for determining the value of u, it is preferable to make a fit at cos0 = 0.5 , which 


gives the expression given by equation (A.90). 


cos 2 0 = (3 cos 0 - 1) at cos0 = 0.5 and 

I(V oj = 1-112 ~ V 2 + u 2 cos0+v 2 |(3cos0-1) 


1- 

( 3 1 

Ui +-v 2 

+ 

( 3 ) 

Ui +-v 2 


l 2 J 


2 

'v J 


^I-Uj+UjCosO where 


( 


»1 = 


3 




u 2 + -v 2 

V 1 J 


(A.90) 


The ratio of mean to central intensity is given by equation (A.91). 



(A.91) 


From equation (A.90) and (A.91) the ratio of l(X , ,0)/F(X / ) is given by equation (A.92). 
l(X',0) _ 1 - Uj + u t COS0 

W) ~ j.i,, (A.92) 

3 1 
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The emittance of solar surface per unit wavelength range is given by equation (A.93). 


E(V) = jcF(V) 


(A.93) 


And the solar flux outside the Earth atmosphere per unit area and wavelength range is 
given by equation (A.94). 


f(V) = E(V) 


R 


solar 


AU 



v solar 

Alf 


1 F(V) 


(A.94) 


From equation (A.90) and equation (A.92) the incident intensity on the coronal electrons 
is given by equation (A.95). 



(A.95) 


In equation (A.94) and equation (A.95) AU is the Sun-Earth distance. The wavelength 
dependent limb-darkening coefficient Uj can be obtained from Astrophysical Quantities 
by Allen 
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A.7 Final expression for the observed intensity 


From figure (A.6) and equation (A.78) all the following expressions are given in 
terms of solar radii, as shown in equation (A.96). 


r — > rR 



(A.96) 


And the observed intensity is given by equation (A.97). 


+o°2n 1 


C(XpR S oiar)= JJ J JdVdtpdcosco d(xR solar ) 


X 


5 0 COS CD* 0 

N,(rR„JxQ?(a,e)x 


— jL — I (A/, co,x)exp 
2VttAb 


x-r 


i + 


2b 2 coscow(rR solar ) 


33 2 


radial 


2Ab 


L V J 

where O = (//, ±) 

// parallel to the radial plane 
L perpendicular to the radial plane 


(A.97) 


The expressions for physical parameters in equation (A.97) in terms of independent 
variables are given in equation (A.98). Equation (A.99) gives the parameters in equation 
(A.97) for which suitable models or actual measurements need to be used. 
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l--u t 
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I© = n - cos -1 (sin oo sin cp sin % + cos cocos x) 


a = sin 


-i 


sincosintp 
v sin(7i- ©) 


X = cos 


ko =sin 


v r y 

fi \ 


-i 


v r / 


0 = sin ^rsinco) 


r = x 2 + p" 


(A.98) 


(X 7 ) = limb darkening coefficient 
f (V) = extraterrestrial solar irradiance 
N e (rR so , ar ) = electron density model 
T(rR sol . ir ) = coronal temperature model 
W(rR SO |ar ) = solar wind model 


(A.99) 
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Appendix B 


CODE 


B.l The basic equation to be solved 

The basic equation to be solved is equation (A.97) as developed in Appendix A 
and reproduced below. The observed K-Coronal intensity is given by equation (B.l). The 
expressions for physical parameters in equation (B.l) in terms of independent variables 
are given in equation (B.2). Equation (B.3) gives the parameters in equation (B.l) for 
which suitable models or actual measurements need to be used. 


+<*>271 1 


Io a <XpR S oiJ= J J J J dV d cos co dcp d(xR solar ) 


X 


: 0 cos co 0 

N,(rR«,JxQ?(a,e)x 


— jL — l(V, co,x)exp 

2VrcAb 


x-v 


1 + 


2b 2 coscow(rR solar ) radia| 


33 2 


2Ab 


L V J 

where O = (//, J_) 

// parallel to the radial plane 
_L perpendicular to the radial plane 


(B.l) 
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u^V) = limb darkening coefficient 
f (V) = extraterrestrial solar irradiance 
N e (rR solar ) = electron density model 
T(rR so i ar ) = coronal temperature model 
W(rR SO ] ar ) = solar wind model 


(B.2) 


(B.3) 
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B.2 Approximations to the J d?C integral 


Consider the exponent in equation (B.l) be represented by y as shown in equation 


(B.4). 


X-rfl + 2 b 2 cosco — 


2Ab 


where A 


q?i' 


(B.4) 


The first approximation to be made on equation (B.4) is given in equation (B.5). 

. qiC qA. 

A = — > — 


(B.5) 


This approximation is justified on the grounds that the speed of light c (3 x 10 5 ) 
km/sec is much larger than the thermal electron velocity q (5500 kin/sec) so that a 
difference of about 10 Angstroms between A, and X' would hardly have any effect on the 
value of A . 

Then differentiating equation (B.4) with respect to X' together with the approximation 
given by equation (B.5) gives equation (B.6). 


w 


, 1 + 2b 2 cos co— 

dy _ c 


dr 


(B.6) 


2Ab 
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The second approximation to be made on the integration limits in equation (B.4) is given 


by equation (B.7). 



Using the limits derived in equation (B.7), the exponent in equation (B.l) together with 
all other X' dependent parameters reduces the integral over J dX' to equation (B.8). 
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B.3 b,Q“ a and Q“ a in terms of the independent variables co, q> and x 

From equation (B.2) the expression for b is given by equation (B.9). 

(B.9) 

Using the trigonometric identity cos(2y) = 2cos 2 (y)-l on equation (B.9) together with 
expressions for cos(tt-0), siny, cosy andr from equation (B.2) gives the following 
expression for b, as shown in equation (B.10). 
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From equation (B.2) the expression for Q“ a is given by equation (B.ll). 



Now, using the relationships sin(2y) = sin 0 and sin a = sin 0)8111 ^ f rom equation (B.2) 

sin 0 

and the trigonometric identities sin(2y) = 2sin ycosy and sin 2y +cos 2 y = l on equation 
(B.ll) reduces equation (B.ll) to equation (B.12). 


Q// a = — — c T (l + sin 2 cosin 2 9-4(1 -cos 2 y)cos 2 y) 

16ji 


(B.12) 


The expression for cosy as a function of (00, 9, x) is given by equation (B.10). 
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Again, from equation (B.2) the expression for Q Ra is given by equation (B.13). 



* , • , , • • sincosuup . ... 

And using the relation sin a = — — from equation (B.2) in equation (B.13) gives 

sin0 

equation (B.14). 

(B.14) 

B.4 Rewriting the basic equation with all parameters expressed in terms of the 
independent variables (co,cp,x,y) 

With change of variables introduced in section B.2 and the relationships 
developed in section B.3 the expression for the observed K-Coronal intensity is given by 
equation (B.15). The expressions for the physical parameters in equation (B.15) in terms 
of the independent variables are given in equation (B.16). Equation (B.17) gives the 
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parameters in equation (B.15) for which suitable models or actual measurements need to 
be used in terms of the given distance p and observed wavelength X . 


(B.15) 
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u j (%) = limb darkening coefficient 

l'(^) = extraterrestrial solar irradiance 

( „ \ 


5 = - 


X\ 1 + 2— by 
c 

i— 

1 + 2b 2 cosco— 


NeW x 2 + p 2 R solar )= electron density model 
T (yjx 2 + p 2 R solar ) = coronal temperature model 
W (y/x 2 + p 2 R SO | ar )= solar wind model 


OB. 17) 


B.5 Changing the independent variable cosco to cos0 


Consider the integral in equation (B.15) as shown in equation (B.18). 

i 

Jdcosco 

COS 0)* 


(B.18) 


From equation (B.2) the expression for cosco* and 0 are given by equation (B.19) and 


equation (B.20), respectively. 


* L . 2 * 

cosco = VI -sin co 

= i 

tr* 

x +p“ 



(B.19) 


(B.20) 
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Substituting cos0 = p equation (B.20) can be written as shown in equation (B.21). 


COS (0 = cos 


f 

< n : — 

• -1 

l-p- 

sin 

y x 2 + p 2 

V 

l )) 


(B.21) 


Differentiating equation (B.21) with respect to |i gives equation (B.22). 



(B.22) 


i 

And the limits on the integration Jdcosoo will change to the following form. 

COSO) 


From equation (B.16) the expression for p can be written as equation (B.23). 


p = cos 0 
= Vl-sin 2 0 

=VT 


2 • 2 

r sin co 


: -y/l - (x 2 + p 2 \l ~ COS 2 CL)) 


when coso) = 1 then p = 1 and 


when coso) = coso) = 1 


—T -- — r then p = 0 
x +p" 


(B.23) 
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Incorp orating the change of variable from co to 0 the basic equation (B.15) changes to 
equation (B.24). ). The expressions for physical parameters in equation (B.24) in terms of 
independent variables are given in equation (B.25). Equation (B.26) gives the parameters 
in equation (B.24) for which suitable models or actual measurements need to be used in 
terms of distance p and observed wavelength X . 


+°°2n 1 +oo 

C(^PK 

solar ) J* J* J* J" d(xR so i ar ) x 

-oo 0 0 

N e (Vx 2 +p 2 R solar )xQ ( R ;(p, 9 ,x)x 

tJx 2 + p 2 ^x 2 + p 2 + p 2 -1 

-^l(y,p,(p,x)- CXp L y | 

1 + 2b 2 (3— 
c 

where O = (//, JL) 

// parallel to the radial plane 
-L perpendicular to the radial plane 


(B.24) 
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(B.25) 


(B.26) 


246 





B.6 Expression for the extraterrestrial solar irradiance !'(£,) 



3000 3500 4000 4500 5000 5500 

ang £3 tor o m s 

Figure (B.l). The plot of extraterrestrial solar irradiance VS 
wavelength obtained with a Fourier Transform Spectrometer at the 
McMath/Pierce Solar Telescope on Kitt Peak, Arizona. 


The spectrum of extraterrestrial solar irradiance shown in figure (B.l) was 
downloaded from ftp.noao.edu/fts/fluxatl. The spectrum was then recreated in intervals 
of 0.0025 Angstroms through linear interpolation. The reasons for this change will be 
apparent in section (B.10). Here f(^) is the wavelength dependent extraterrestrial solar 
irradiance with wavelength q measured in Angstroms. 
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limb darkening coefficient 


B.7 Expressions for the limb-darkening coefficient Uj (^) 


limb darkening coefficient VS wavelength 



wavelength in angstroms f 


Figure (B.2). The plot of limb-darkening coefficient versus wavelength 
with linear approximation super imposed. 
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Using the wavelength dependent limb-darkening coefficient from Astrophysical 


Quantities by Allen the following linear fits were made for the different wavelength 
regions. 


(1) 3000 - 3200 Angstroms 

^ 538.0 0.11 x^ 

u i(s) = - 

1 200.0 200.0 


(B.27) 


(2) 3200 - 3700 Angstroms 


u,ft) = 0.922 


(5- 3220.0) x 0.120 
480.0 


(3) 3700 - 5000 Angstroms 


(B.28) 


u 1 (4) = 0.862- 


fe- 3700.0) x 0.212 
1300.0 


(4) 5000 - 5500 Angstroms 


(B.29) 


“i(y 


625.0 0.06 x £, 

500.0 500.0 


(B.30) 


In the above equations £ is the wavelength measured in Angstroms. 


249 



B.8 Expression for the coronal electron density model N e (-y/x 2 + p 2 R solar ) 

If the eclipse takes place at sunspot maximum, the corona has an approximate 
circular form. This means that the coronal rays are most probably equally distributed 
between the polar and the equatorial regions of the Sun. Since the Sun is expected to 
reach sunspot maximum in the year 2000 the solar corona at the August 11, 1999 total 
solar eclipse will assumed to be of circular form. This assumption justifies using a radial 
dependent expression for the electron number density. 

From the analysis of the photometric data on ten eclipses from 1905 to 1929 
Baumbach (1937) deduced the following expression for the coronal electron density. 


N e (r) = 10 s (o.036r _1 ‘ 5 + 1.55r“ 6 + 2.99r 16 )cm 3 

(B.31) 

where r = yx 2 +p 2 


Cram (1976) used the expression given in equation (B.32). 

N e (r) = 1.67 x jo (4+4 - #4/r) cm" 3 
where r = -y/x 2 + p 2 


(B.32) 


Equation (B.32) agrees within 2% in the region 1.5 < r < 2.0 to the minimum equator 
model given by Van de Hulst (1950). In this code the electron density model given by 
equation (B.31) was used. Figure (B.3) shows the electron density plots given by 
equation (B.31) and equation (B.32). 
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Electron number density VS radius 



Figure (B.3). The plot of electron number density 
versus the radial distance from the solar surface for 
the two models given by equation (B.31) and 
equation (B.32). 


B.9 Expressions for the coronal temperature model T {Jx 2 + p 2 R solar ) and the 
solar wind model W {<Jx 2 + p 2 R solar ) 

In this code the corona is considered to be isothermal which implies that that the 
temperature T is a constant. As regards the solar wind W, it is considered to be radial, 
isotropic and constant. The idea is to construct multiple models for different 
combinations of temperature T and solar wind velocity W. Nevertheless, the code does 
not restrict from applying a coordinate dependent temperature and wind profiles that are 
measured independently. 
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+o°2K 1 +oo 


I * a (?i,p) = constant x J J J Jdy dcpdpdxx 


f j 1.5 , \-6 I 16 \ 

10 8 x 0.036xL/x 2 +p 2 + 1.55xU/x 2 +p 2 + 2.99xL/x 2 +p 2 ) x 


Vx 2 + p 2 Vx 2 +p 2 +p 2 -: 

(l-(l-p 2 )sin 2 cp)x 

l-u t fe) + ^Uife) f ^ x 

expf - y2 l where 


l + 2b 2 [3 
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3xq t (AU) 2 
167r s/2 R.„ 1 


. i L I i — (3 2 . p P 
S{ 1* 2 +p 2 V777 

f ( 1 1- 2 Yl 

P = cos sin” 1 — 

\x 2 +p 2 


X 1 + 2 — by 
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l + 2b P 


20 W 
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o T =6.677 xl0“ 25 cm 2 
R solar = 6.96 x 10 10 cm 
AU = 1.496 xlO 13 cm 


From equation (B.33) and equation (B.34) the total observed intensity is given by 
equation (B.36). 

It„,(vp)=i;-(vp)+ii - (vp) 


fOO 

1. Evaluation of the J" A(x)dx integral 


The above integral is evaluated using a 20-point Hermite polynomial expansion as 
shown in equation (B.37). 


jA(x)dx= j(A(x)e x2 )r x2 dx 
= f A(h,>^w, 


where h t and w, are the Flennite polynomials and the associated weights, respectively. 
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271 

2. Evaluation of the J B(cp)d(p integral 

o 

This integral is evaluated using 20-point trapezoidal composite rule as shown in 
equation (B.38). 

(B.38) 

3. Evaluation of the J C(p)dp integral 

0 

This integral is evaluated using 10-point trapezoidal composite rule as shown in 
equation (B.39). 

(B.39) 

+oo 

4. Evaluation of the | D(y)e ' dy integral 

This integral is evaluated using a 12-point Hermite polynomial expansion as shown in 
equation (B.40). 

|D(y)e“ y2 dy = £ D (h,)wi (B. 40 ) 

1=1 

where h, and w, are the Elermite polynomials and the associated weights, respectively. 
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T" 

Consider the integral J D(y)e” v dy its y dependent parameters and the expansion outlines 


in equation (B.40) in equation (B.41). 


J, Xy ) e V dys jizr4±l^ ffe ) e -r dy 


1 


where s = ^(y) 


g £ i-»,te,) + n»,te,) ffe>| 


1=1 




=1 — “ - -- ffe.V, + di u ; fel) ffe.v, 

l=1 l--Uifei) 1=1 1 - — Ulfel ) 


AJ 


6 q ^ 

1 + 2 — bh, 


where = 


1 + 2h 2 (3 


(B.41) 


Due to the rapid intensity variations associated with the Fraunhofer lines, as 
evident in figure (B.l), this effect is included by taking the mean value for l'(q, ) over 
each interval in the Flennite polynomial expansion coefficients h, . The mean is then 
given by equation (B.42). 
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In order to evaluate the integral in equation (B.42) the following procedure was 
followed in order to expedite the calculation. 


1. The extraterrestrial solar irradiance spectrum was recreated in intervals of 0.0025 
Angstroms (r|'). 
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2. For a given set of values for (x,cp,[i) from the integrals ( J d x , J d(p, J d Lt ) the values 
(n upper » n lower >B R ) were calculated. 

3 - (n upper » n lower ) were assigned new values (n' pper , bLer ) which are divisible by 0.0025 
and closest of (V) to (b uppe r Slower )• 

4. The values (f(r| upper ),f (blower )) now corresponded to (b' upper iVL* ) • 

5. Tliis procedure allowed the integral to be expanded in intervals of 0.0025 Angstroms. 

6. Within the range of 0.0025 Angstroms the limb-darkening coefficient U, 
extraterrestrial solar irradiance f could be assumed to be a constant. This allows for 


writing equation (B.42) as shown in equation (B.43). 


^ R Slower 1 - — U , (b) ^ R Tll< ™" 1 - — U| (b) 


1 

b R 


T1 lower +0-0025 A” . ( \ T| lower +0-0050 A” , / \ 

J 2^ahf(n)in+ J -f^-f(n)an 

^1 lower I — 11 (Vl 1 +1 lower +0-0025 A 0 I II . (n) 


1- 3 U i( t >) 


Tltower +0-0025 A 0 1 - -- |J , (b) 


,+ 


+ 


Br 


T 

T1 upper - 0*0025 A 0 l-i U[ (b) 

Slower +0-0025 A° ( \ +0.0050 A 0 ( \ 

I — — f (n)dri + I -^_f(p)dp+ 

Slower 1 U | (b) Slower +0-0025 A “ ] U | (b ) 


,+ 


Y u,(n) 

J 1 / S 

-0.0025 A 0 1 — — U , (b ) 


f(b)db 


(B.43) 
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Now to elaborate on the method, by which the integration is expedited, consider the 


following table. Columns 1 and 4 are from section B.6 and columns 2 and 3 are from 
section B.7. Here T is divisible by 0.00025. 


Table (B.l). The dependence of the extraterrestrial solar irradiance and the 
limb-darkening functions with the wavelength.. The columns 1 and 4 are from 
section (B.6) and columns 2 and 3 are from section (B.7). 


Wavelength r\ 
in Nanometers 
(nm) 


«i(n) 

Extraterrestrial 

Solar 

irradiance f(r|) 



T 

A1 

B1 

Cl 

T+0.00025 

A2 

B2 

C2 

T+0.00050 

A3 

B3 

C3 

T+0.00075 

A4 

B4 

C4 

T+O.OOIOO 

A5 

B5 

C5 

T+0.00125 

A6 

B6 

C6 

T+0.00150 

A7 

B7 

C7 

T+0.00175 

A8 

B8 

C8 

T+0.00200 

A9 

B9 

C9 

continued 

continued 

continued 

continued 


In the code T is 300.00025 nni and then increases in intervals of 0.00025 nm to 
550.00000 nm. 


259 








Now create a new table with each row reflecting the cumulative sum of table (B.l) for the 


columns 2, 3 and 4. 


Table (B.2). Cumulative values of the columns 2,3 and 4 of table (B.l). 


T 

A1*C1 

B1*C1 

Cl 

T+0.00025 

A2*C2 

B2*C2 

C2 

T+0.00050 

A3*C3 

B3*C3 

C3 

continued 

continued 

continued 

continued 


In the code the following format is used, where columns 2 and 3 are multiplied by 
column 4 of table (B.2). 


Table (B.3). Multiplying column 4 of table (B.2) with column 2 and 3 of table (B.2). 
This is the format of the table used by the code. 


T 

A1*C1 

B1*C1 

T+0.00025 

( A1 *C 1 )+( A2*C2) 

(B1*C1)+(B2*C2) 

T+0.00050 

(A1*C1)+(A2*C2)+ 

(A3*C3) 

(B 1 *C 1 )+(B2*C2)+ 
(B3*C3) 

T+0.00075 

( A1 *C 1 )+(A2*C2)+ 
(A3*C3)+ (A4*C4) 

(B1*C1)+(B2*C2)+ 
(B3*C3)+ (B4*C4) 

continued 

continued 

continued 
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Going back to equation (B.43) the integrals could be written as a summation, as 
shown in equation (B.44). Figure (B.4) illustrates the evaluation of the first integral of 
equation (B.44) where the integration under each interval is approximated as the area 
under the square. This procedure is justified on ground of the very small interval. 


upp 

J 


1 — U, (r|) 


f(n)dTi + — J 

n J 


«i(b) 


-f(Tl)dTl 


J 1 \ * / * J 1 

Hlower 1 - — U, (r|) ^ R ^'ower 1 - — U, (r|) 

s Ati_ ^ 1 ItAii n £' _u i (n)_f( T1 ) 

n=nL„ i_i_ n (n) n=n[„„er i _i_ n fn) 

in intervals ~ w 1 \ 1 1 / in intervals A ^ u 1 V 1 1 / 

•J nf 0.0025 A 0 * 


of 0.0025 A u 


where Ar| = 0.0025 A® 


Tlu =B 


upper 


I) 


lower 


(B.44) 



Figure (B.4). Blustrating the integration procedure of the first 
integral in equation (B.44). 
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Taking the cumulative value of the parameter reflected on the y-axis of figure (B.4) 


would give the following, as shown in figure (B.5). 



Figure (B.5). Plot of the cumulative y-axis values of figure (B.4) against the 
wavelength t|'. 


The same procedure as above could be performed for the second integral in equation 
(B.44) too. 
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Ill terms of expediting the calculation of the integral presented in equation (B.42), 
equation (B.44) allows the code do the calculation in two easy steps in conjunction with 
table (B.3) and illustrated by figure (B.5) for the first integral, as follows. 


1. Find the corresponding row numbers, say, rows (R,P), respectively, which equal the 
values (r|' pp(;r , rq ower ) in column 1 of table (B.3) and as depicted in figure (B.5). 


2. Find the values in columns 2 and 3 corresponding to the rows (R,P) in table (B.3) and 
as shown by the arrows cutting across the y-axis in figure (B.5). 


3. Now calculate [(2(R) - 2(P)) + p(3(R) - 3(P))]x — where the numbers correspond 

% 

to the column numbers. The first part 2(R) - 2(P) corresponds to the first integral in 
equation (B.44) and the second part 3(R) - 3(P) corresponds to the second integral 


in equation (B.44). 


In figure (B.5), 


(2(R)-2(P)) 

Fr 


evaluates the first integral in 


equation (B.44). 
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B.ll The symbols used in the code 

The following symbols have been used in the code. 

x — > x 

<p^y 

(X — » z 

r| — » cetaU 

upper 

rj — > cetaL 

lower 

ri' —> a(0,u) 

upper 

r\' -^a(0,t) (B.45) 

lower 

h, — > hermite_12(l) 
w, — » weights _12(l) 
hj — > hermite_20(i) 

Wj — > weights _20(i) 
b-^bb 
p — » rho 

X -> lambda 

In the code table (B.3) is given by 3000_5500_cum_irradiance.dat 
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From equation (B.2) the following relationships could be derived. 



(B.46) 


(B.47) 


(B.48) 


angle _ 4 = sin(%) 

= sin(cos (angle _ 3)) 


(B.49) 


angle _ 5 = cos(7t - 0) 

= cos(co)cos(%) + sin(co)sin(%)sin (9) 

= angle _ 2 x angle _ 3 + angle _ 1 x angle _ 4 x sin (y ) 
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(B.52) 


(B.53) 


(B.54) 
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B.12 The IDL code 


The following is the code written in IDL where (;) is a comment symbol and ($) a 


command continuation symbol. 

;README 

;DECIDE ON THE RADIAL WIND VELOCITY IN KM\SEC 
;DECIDE ON THE ISOTHERMAL CORONAL TEMPERATURE IN MK 
;DECIDE ON THE VALUES FOR RHO IN SOLAR RADIUS 
; 000.0 0.50 1.10 
; 400.0 2.00 1.30 
; 800.0 1.50 1.40 

;ENTER THE ABOVE INFORMATION IN BATCH FILE.PRO 

;COMPILE ECLIPSE SIMULATION.PRO 
;RUN ECLIPSE SIMULATION.PRO 

;THE DATA WILL BE STORED AS 
;’ES wOOO.O T0.50 rl.10.dat’ 

; ’ESw400.0 T2.00rl.30.dat’ 

; ’ES_w800.0_T1.50_rl.40.dat’ 


;CALULATING AVERAGE INTENSITY IN THE REGION [cetaU,cetaL] 

function KURUCZ_TABLE, a, z,cetaU,cetaL 


u=long(round((doiible(cetaU)-300.00025d0)/0.00025d0)) ;a(0,0)=300.00025 
t=Iong(roiind((double(cetaL)-300.00025d0)/0.00025d0)) ;a(0,0)=300.00025 
interval =a(0,u)-a(0,t) 

iiradiance=((a(l,u)-a(l,t))+z*(a(2,u)-a(2,t)))*0.00025d0/iiiterva 

; microwatts/cm squared/nm 

return, irradiance ; microwatts/cm squared/nm 

end 

CALCULATE FOURTH INTEGRAL PARAMETRS AND THE VALUE 

function calculate_fourth_integral_value ,a, lambda, rho,hermite_l 2, weights...! 2, x,y,z,w,T 
q =5508.()*sqrt(T) ;thermal electron velocity in km/sec 
c =3.0d5 ; velocity of light in km/sec 

final_sum_4==dblarr(2, 1 ) 

angle.. 1 =sin (acosf z))/sq rt(x A 2+ rh o A 2) ;sin(omega) 

angle_2=cos(asm(angle„l)) ;cos(omega) 
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;cos(chi) 
; sin (chi) 


angle 3=x/sqrt(x A 2+rho A 2) 

angle 4=sm(acosfaiide 3)) 

angle_5=angle_2*angle_3 + angle.. J *angle_4*sin(y) 
angle_6=acos(angle„5) 

r=sqrt(x A 2+rho A 2) ;radius 

aa=l .0d8*(0.036/(r A l .5) + 1 .55/(r A 6) + 2,55/(r A 16)) :Baumbach model 

bb=cos(angle_6/2.0) ;cos(gamma) 

cc_t= 1 .Q-(sin(y) A 2)*(angle 1 A 2) ;tangential 

ec....r=ld)+(sm(y) A 2)*(ang]e....l A 2)-4,0*(l ,0-bb A 2)*(bb A 2);radial 
dd=(l .0*z)/(sqrt(x A 2+rho A 2)*sqit(x A 2+rho A 2+z A 2- 1 .())) 

;changing d(cos(omega))to d(mu) 

ee=l .0+2.0*bb A 2*angle„2*(w/c) ;wind component 


N=1 1 ; Nth order hermite polynomial expansion is used 

sum J... tangenti al =0.0 
sumJLradial =0.0 


for 1=0, N do begin; first Joop 


if (i eq 0 ) then begin 

cetaL=(l.ambda-lambda*(q/c)*bb*hemiite_12(i))/ee 

;nanometers 

end else if (i ne 0) then begin 

cetal^(iambda-lambda*(q/c)*bb*(herniite_12(i)+hermite_12(i-l)))/ee 

manometers 

endif 


if (i eq 11) then begin 

cetaU=(lambda-lambda*(q/c)*bb*hennite__12(i))/ee 

;nanometers 

end else if (i ne 11) then begin 

cetaU=(lambda-lambda*(q./c)*bb*(liermite_12(i)+hermite__12(i+l)))/ee 

manometers 

endif 


if (((cetaU or cetaL) It 300.00025d0) or ((cetall or cetaL) gt 550.00000d0)) then begin 
tinal_sum_4(0,0)=999999999999.9999 
final_siim_4( 1 ,0)=999999999999.9999 
goto, terminated 
endif 
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CALCULATING THE AVERAGE INTENSITY CORRESPONDING TO 
;[cetaU,cetaL] 


irradiance=KURUCZ__TABLE(a,z,cetaU,cetaL) ;microwatts/cm squared/nm 
FC=irradiance* LQd-7 ;W/cm squared/ Angstroms 


sum_i_tangential=sum_i_tangential+weights_12(i)*(FC) 
sum_i_radial =sum_i_radial +weights_12(i)*(FC) 


endfor ;loop_i 


final_sum_4(0,0)=sum_i_tangential*aa*cc_t*dd*(1.0/ee) 
final_sum_4( 1 ,0)= sum_i_radial*aa*cc_r*dd*( 1 ,0/ee) 
term in ate... 4: 
return, fmal_sum_4 
end 


CALCULATE THIRD INTEGRAL PARAMETRS AND THE VALUE 

function calculate_third_integral_value ,a,lambda,rho,hemiite_12,weights_12,x,y,w,T 

; THIRD .INTEGRAL BOUNDS 

a3=0.0 

b3=i.O 

n3=10.0 

h3=(b3-a3)/n3 

f3=fix(n3) 

sum__z_t=dbl arr(t3 + 1 } 
sum_z_r=dblarr(f3 + 1 ) 
total_sum_z_t=0.0 
total__sum__z__r=0.0 
final_sum_3=dblarr(2, 1 ) 

for r=0,£3 do begin 

z=a3+r*h3 ;mu=cos(theta) 

final. ...su m_4=calculate__fourth„integra Lvalue $ 

(a, lambda, riioJh.ermite_12,weights_12,x,y,z,w,T) 
if ( final_sum_4(0,0) eq 999999999999,9999) then begin 
final_sum_3=final_sum_4 
goto, terminated 
endif 

sum_z__t(r)=final_.sum„4(0,0) 

total sum z t=to la I sum z t+i'i naLsu m... 4f0,0) 
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sum_z_r(r)=final_sum_4(l,0) 

totaLstim__z__r==totaLsum__z„r+fma!__sum 4(1,0) 

endfor 


final_sum__3(0,0)=(h3/2,0)*(2.0*total_suni„z_t-sum_z„t(0)-sura„z_t(n3)) 

fmal__sum_.3(l,0)=(li3/2,0)*(2.0*total_sum__z_r-sum__z__r(0)-sum__z__r(ii3)) 

terminate__3: 

return, final . sum 3 

end 


CALCULATE SECOND INTEGRAL PARAMETERS AND THE VALUE 

function calculate, second,integral„yalue ,a,lambda,rho,hermite_J2,weights_.12,x,w,T 


; SECOND INTEGRAL BOUNDS 

a2=0.0 

b2=2.0*!pi 

n2=20.0 

h2=(b2-a2)/n2 

f2=fix(n2) 

sum_y_t=dblarr(f2+ 1 ) 
sum_y_r=dbl arr(f2+ 1 ) 
total_smn_y_t“0.0 
total_sum_y_r=0.0 
final,, sum,. 2=dblarr(2, 1 ) 

for v=0,f2 do begin 
y=a2+v*h2 ; in radians 
.fin.al_sum_3=calcuiate_third_integral_value $ 

(a, lambda, rho,hermite„l 2, weights,. 1 2,x,y,w,T) 
if ( final, sura, 3(0,0) eq 999999999999,9999) then begin 
fmal_sum_2=final_sum_3 
goto, terminate_2 
en d if 

sum_y_t(v)=final_sura_3(0,0) 

total_sum_y_t=total_sum_y_t+finaLsum_3(0,0) 

sum y, r(v)==final sum 3( 1 ,0) 

tot al_sum„y„r=total„sum_y , _r+final_su m_3 ( 1 ,0) 


endfor 


tinal_sum_2(0,0)=(h2/2.0)*(2,0*total_sura_y_t-sum_y_t(0)-sum_y_t(n2)) 

final,.stirn,2(1 ,())=( h2/2.())*(2.()*totaLsuniy„r-suTn y.,r(0)-sum y r(n2)) 

terminate, .2: 
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return, final_sum_2 
end 

CALCULATE FIRST INTEGRAL PARAMETERS AND THE VALUE 

functioncalculate_first_integral_value $ 

, a, lambda, rho,hermite_ 1 2,weights_ 1 2,hermite_20,weights_20,w,T 

;FDRST INTEGRAL, USING HERM1TE POLYNOMIALS 
n 1=19 

sum_x_t=dblarr(iil+l ) 
total_sum_x_t=0.0 
s um_x_r=dbl arr(n 1 + 1 ) 

total sum x_..r=0.0 

intensity=d blarr (2, 1 ) 

for h= 0,nl do begin 
x=hermite_20(h) 

final_sum_2=calculate_second_integral.„value $ 

(a, lambda, rho,hermite... 12, weights. ...1 2,x,w,T) 
if (fmal_sum_2(0,0) eq 999999999999.9999) then begin 
intensity=final_sum_2 
goto, terminate....! 
endif 

sum„x_t(h)=final_siim_2(0,0)*weights„20(h)*exp(x*x) 

total_sum_x_t==total...sum._x_t+snni x.j(h) 

sum_x_r(h)=fmal_sum_2(l,0)*weights...20(h)*exp(x*x) 

total_sum_x_i^total_smm_x_r+sum_x_r(h) 

end for 

sigma=6.677E-25 ;Thomson Scattering Cross Section cm squared 

AU = 1 .49597870E+ 1 3 ;cm-symbol D in CRAM 

solar... radius=6.96E+ 10 ;cm- symbol R in CRAM 
constant=(3.0*sigma*(AU) A 2.0)/(16.0*solar_radius*(!pi) A (2.5)) 

.i.ntens.i.ty(0,0)=constant*total_sum_x_t ;Watts/cm squared/Angstroms/steradians 
intensity! 1 ,0)=constant*totai„su m_x_r ; Watts/cm squared/Angstroms/steradians 

terminate... I : 
return, intensity 
end 


print, ’starting time systime(O) 
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;FILES TO BE READ 

; READING THE KURUCZ EXTRATERRESTRIAL IRRADIANCE AT 0.00025A 
INTERVALS FROM 3000A TO 5500A 


a=dblarr(3, 1 0000001) 

opem;2, 3OOO_550O_cirniJhxadiance.dat’ 

readf,2,a 

close, 2 

; READING THE ROOTS OF 12TH ORDER HERMITE POLYNOMIALS 
;HANDBOOK OF TABLES OF MATHEMATICS: SELBY 3RD EDITION PAGE 
;843 


liennite. ... 1 2=dblarr( 1,12) 
openr, 3, , hermite_ 1 2_roots. dat ’ 
readf,3 ,he:nmite__12 
close, 3 


; READING THE WEIGHTS OF 12TH ORDER HERMITE POLYNOMIALS 
;HANDBOOK OF TABLES OF MATHEMATICSrSELBY 3RD EDITION PAGE 
;843 


weights.....! 2=dblarr( 1 ,12) 

opem;4,lieri:riite 12_wei.ghts.daf 

readf,4,weights....l2 
close, 4 


; READING THE ROOTS OF 20TH ORDER HERMITE POLYNOMIALS 
; HANDBOOK OF MATHEMATIC FUNCTIONS EDITED BY M.ABRAMOWITZ 
;AND L.A. STEGUN PAGE 924 


hermite_20=dblarr( 1 ,20) 

openr, 5, ’hermite 20.roots.dat ’ 

readf,5,hermite_2Q 
close, 5 
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; READING THE WEIGHTS OF 20TH ORDER HERMITE POLYNOMIALS 
; HANDBOOK OF MATHEMATIC FUNCTIONS EDITED BY M.ABRAMOWITZ 
;AND L.A. STEGUN PAGE 924 

wei ghts 20=dbl arr( 1 , 20) 

opem*,6,Tieraiite__20_.weights.dat ! 

readf,6,weights__20 

close, 6 


; READING THE BATCH FILE.PRO 

;THIS FILE CONTAINS THE WIND, TEMPERATURE AND RHO VALUES 

data=read__ascii( batch__fi.le.pro ’,$ 
c omment_symbol = *xxx.x ’,d ata_start= 1 ) 
data_elements=n_ elements(dataiield 1 ) 


counter=(data_elements/3)-l 


for j=0, counter do begin 


wind= data.fieldl (3*j+0);wind in kmVsec 
temp= data.fieldl (3*j+l);Temperature in MK 
rhoo= data.fieldl(3*j+2);rho in solar radius 

str_w=strmg(wind) 
parts=str_sep(str_w, ’) 
wl=parts[0] 
w2=strm.id(parts[ 1 ] ,0, 1 ) 

str_t=string(temp) 
parts=str_sep(str_t, ) 
tl=parts[0j 

t2=strmid(parts[ 1 ] ,0,2) 

str_r=string(rhoo) 
parts=str _ septstr _ r, 3 
rl=parts[0] 

r2=strmid (parts) 1 ] ,0,2) 


w=double(wind) 

T=double(temp) 

rho=double(rhoo) 

print,' ’wind = IwdkmXsec temp = ’,temp,’MK$ 
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; RECORDING THE OUTPUT 

open w, 9, ES_w’+trira(w 1 , ’(i3.3)’)+’, 5 +trim (w2, ’(i 1.1) )+’.„TS 
+trim(t 1, \i 1.1) )+ ’+trim(t2, ’(12.2) ’)+’„r'+trira(rl , ’(il . 1 ) ’)+ $ 
+trim(r2 , ’(12.2) ’)+ ’.dat ’ 

for r=0,l(X)0 do begin ;370.0 nm to 470.0 nm 
lambda=370.0d0 + r*0.1d0 ;nano meters 


iiitensity=calculate._first„integral__valiie $ 

(a, lambda, rho,liemiite„12,w'eights„12,hemiite__20,weights__20,w,T) 
if (intensity(0,0) eq 999999999999.9999) then begin 
print, lntensity(0,0) = 999999999999.9999 for $ 
lambda = lambda 
goto, inten sity„out„of_range 
endif 


totai_mtensity=intensity(0,0)-f-intensity( 1,0);' Watts/cm squared/ Angstroms/steradians 

print, lambda- lambda* 10.0 Angstroms’,’ $ 
intensity = total Jntensity* 1 .QE+7 ,$ 

’ ergs/sec/cm squared/Angs/steradians’ 

;print,lime = ’,systime(0) 

printf, 9, format= ’(4E15.5)’, lambda* 10.0,total_intensity,intensity(0,0),$ 
intensity(l.O) 
flush, 9 

intensity_out_of_range: 
endfor ; the r loop for intensity 

close, 9 

endfor ; the j loop for batch_file 

print, ’ending time ’,systime(0) 
end 



All example of the input data information given in BATCHFI LE .PRO 


Table (B.4). The format of the input information on the wind velocity in 
km/sec, isothermal temperature in MK and the distance to the line of 
sight in SR. 


Wind(km/sec) 

Temperature! MK) 

RHO(SR) 

400.0 

1.50 

1.10 

300.0 

0.50 

1.50 

continued 

continued 

continued 

xxx.x 

XX. X 

XX.X 
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Appendix C 


A STREAMER MODEL 
C.l The effects due to streamers 

Due to the structural nature of the solar corona it is important to consider the 
effects on the K-corona due to coronal inhomogeneities. In this regard consider a 
streamer imbedded in a spherically symmetric isothermal corona. The model streamer 
considered has the following structure as depicted in figure (C.l). The streamer is of 
thickness d and the extremes of the streamer, namely, RT and UV are equidistant from 
the solar center C and are inclined at an angle 8 to the plane of the solar limb. Within the 
streamer the electron number density is enhanced by a factor f and the thermal electron 

s 

temperature is T s . From figure (C.l) the expressions for Xi (=FV) and X 2 (=FT) are 
given by equation (C.l). 


pxtan(e)- 
pxtan(e) + 


d/2 

cos(e) 

d/2 

cos(e) 


and 


(C.l) 
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Figure (C.l). Model streamer of thickness d in front of the solar 
limb with an electron density enhancement by a factor f s and 
thermal electron temperature T s . 


And in the streamer model considered the electron density is enhanced by a factor f s 
in the region between x 2 and x, along the line of sight. The incorporation of the effect 
due to the streamer in the code, as described in section (B.12), is as follows. 

1. For specified values for d, p and 8 calculate x 2 and x 2 . 

2. In the x-integration (integration along the line of sight) if x lies between and including 
x, and x 2 then multiply the electron number density by a factor f s and substitute T s 

for T where f s and T s are specified. 
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3. Perform the x-integration by trapezoidal composite rule instead of the Hermite 
polynomial expansion. This was found to be necessary in order to pick up sufficient 
x-values in the region [ x, and x 2 1. Since the Hermite polynomial expansion 
considered in the code is of twentieth order there may not be sometimes any x-values 
in the region [ x 2 and x 2 ] . 

4 . The x-integration is substituted by the following commands 

CALCULATE FIRST INTEGRAL PARAMETERS AND THE VALLE 

function calculate_first_integral_value $ ,a,lambda,rho,hermite_12,weights_12,w,T 

; FIRST INTEGRAL BOUNDS 

al=-6.0 ;lower bound at - 6.0 solar radii 

h 1 =+6.0 ;upper bound at +6.0 solar radii 

nl=120.0 ; region divided into 120 equal intervals 

hl=(bl-al)/nl 

fl=fix(nl) 

sum_x_t=dblarr(f 1 + 1 ) 
sum_x_r=dblarr(f 1 + 1 ) 
total_sum_x_t=0.0 
total_sum_x_r=0.0 
intensity=dblarr(2, 1 ) 
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for h=0.fl do begin 
x=al+h*hl 

final_sum_2=calcu]ate_second_mregral_value $ 

(a,lambda,rho,hermite_.12,weiglits„12,x,w,T) 

if (final....sum 2(0,0) eq 999999999999,9999) then begin 

mtensity=fiiiaLsum„2 
goto, terminate,...! 
endif 

sum...x.....t(!i)=iinal,_sum__2(0,0) 

total... sum...x.„. t=to tai suit i_x_t+su n.i__x t (h) 

sum_x_r(h)=final_sum_2(l,0) 

totaI_sum_x_r=total_sum_x_r+sitm. x....r(li) 

endfor 

Although x-integration need to be considered in the interval [ - °o,+oo ] it is 
sufficient to restrict the integration to [-6.0, +6.0] solar radii due to the rapid falloff in the 
electron number density. 
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Consider a streamer model as depicted by figure (C.l) with a streamer width (d) 
of 0.5 solar radii, density enhancement factor (f s ) of 4.0 within the streamer and 

assuming the thermal electron temperature (T s ) within the streamer to be equal to 
isothermal coronal temperature (T) for the rest of the corona. Figure (C.2) shows the 
effect of streamers, at angles ( e ) - 45.0 degrees (front of the limb) and + 45.0 degrees 
(behind the limb) to the limb of the solar corona, on the K-coronal spectrum. 


observed total intensity VS observed wavelength 



Figure (C.2). Plot to show the effect of streamers on the K- 
coronal spectrum with d=0.5 SR, f s =4.0, T S =T=1.0 MK and at 

angles 8 =-45.0 and +45.0 degrees. 
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Although the two streamers considered in figure (C.2) are symmetric about the 
plane of the solar limb the asymmetry in the two theoretical spectra is due to the fact that 
the streamers scatter through different angles in reaching the observer. Figure (C.3) 
illustrates this point. 



Figure (C.4) shows the effect of streamers, at series of angles (8) in front and 
behind the plane of the solar limb, on the K-coronal spectrum. As per the streamer model 
considered the intensity of the K-coronal spectrum decreases as the location of the 
streamer moves away from the plane of the solar limb. 
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observed total intensity VS observed wavelength 



The following figures (C.5) to (C.6) show the effect on the wind and temperature 
sensitive intensity ratios I(4233)/I(3987) and I(4100)/I(3850), respectively, by a streamer 
of the shape depicted by figure (C.l) and with a streamer width (d) of 2.0 solar radii, 
density enhancement factor (f s ) of 4.0 within the streamer and assuming the thermal 

electron temperature (T s ) within the streamer to be equal to isothermal coronal 
temperature (T). The horizontal lines show the intensity ratio for a given solar wind 
speed, isothermal coronal temperature and line of sight in the absence of streamers. The 
curves are the corresponding ratios with a streamer of the shape depicted in figure (C.l) 
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Figure (C.5). Effects by streamers on the ratio I(4233)/I(3987) for 
W=400 km/sec and rho=l.l SR. The horizontal lines show the 
intensity ratio for a given solar wind velocity of 400.0 km/sec and 
line of sight at 1.1 solar radii, and for different isothermal coronal 
temperature in the absence of streamers. 


From figure (C.5) it is apparent that even for a simple streamer structure shown in 
figure (C.l) the wind sensitive intensity ratio can differ from the ratios predicted by 
models that do not consider streamer structures. The horizontal lines in figure (C.5) are 
from models that do not consider any streamer structures. A similar occurrence could also 
be seen for the temperature sensitive intensity ratio as shown in figure (C.6). 
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Figure (C.6). Effects by streamers on the ratio I(4100)/I(3850) for 
W=400 km/sec and rho=l.l SR. The horizontal lines show the 
intensity ratio for a given solar wind velocity of 400.0 km/sec and 
line of sight at 1.1 solar radii, and for different isothermal coronal 
temperature in the absence of streamers. 


Unfortunately it is not an easy task to incorporate realistic streamer models simply 
by observing an eclipse. Nevertheless figure (C.5) and figure (C.6) are examples that 
show the effect on the wind and the temperature sensitive intensity ratios by streamers. 
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